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SECTION  1 
INTRODUCTION 

This  report  is  a  user's  guide  for  the  utilization  of 
several  interdependent  computer  programs  which  are  designed  to 
provide  the  survivability/vulnerability  engineer  with  a  tool 
for  the  utilization  of  finite  element  models  (FEM's)  in  the 
simulation  of  structural  flight  loading  of  wings  and  wing 
components.  This  Structural  Analysis  Computer  Program  User's 
Manual,  in  conjunction  with  References  1  and  3,  is  designed  to 
completely  describe  the  utilization  and  interaction  of  four 
computer  programs  for  predicting  the  static  response  of 
undamaged  and  ballistically  damaged  wing-like  structures.  The 
objectives  of  utilizing  a  simplified  structural  analysis  tool 
are  detailed  in  Reference  3  but  may  be  summarized  as  follows: 

(a)  to  provide  a  check  on  experimentally  obtained  data, 

(b)  to  compute  the  internal  stresses  in  an  undamaged 
structure, 

(c)  to  predict  the  stress  redistribution  in 
ballistically  damaged  structures,  and 

(d)  to  estimate  the  threshold  and  residual  strength  of 
ballistically  damaged  aircraft  wing  and  empennaqe 
components . 

These  computer  programs  employ  the  finite  element 
modeling  (FEM)  technique  to  discretize  a  wing  structure  into 
its  various  components  in  a  fashion  that  provides  an  accurate 
mathematical  representation  of  the  wing.  This  model  is  then 
subjected  to  simulated  experimental  loading  to  determine  the 
information  desired  in  items  a-d  above. 

This  section  will  provide  a  brief  introduction  to  finite 
element  modeling  (FEM)  and  will  then  discuss  briefly  each  of 
the  component  programs  in  the  FEM  (finite  element  model)  approach 
to  simulated  flight  loading.  The  remainder  of  this  report 
addresses  each  of  the  components* of  FEM  in  turn:  preprocessing 


(WINGEN  in  Section  2);  structural  analysis  (MAGNA  in  Section  3); 
and  postprocessing  (PLOTBOB  and  CONTOUR  in  Section  4).  Included 
at  the  end  of  this  report  are  several  appendices  illustrating 
procedures  for  executing  the  programs  covered  in  this  report. 

A  glossary  of  terms  is  also  provided  with  cross-referencing 
to  help  the  user  with  terminology. 

1.1  FINITE  ELEMENT  MODELING 

The  finite  element  method  is  a  systematic  process  of 
reducing  a  structure,  such  as  the  one  illustrated  in  Figure 
1.1.1,  into  a  set  of  interrelated  nodes  and  elements.  A  node 
is  a  three  coordinate  point  in  space.  A  series  of  nodes 
are  connected  together  to  form  mathematical  (structural) 
components  or  finite  elements.  These  components  are  then 
assembled  to  construct  the  FEM  (finite  element  model) .  Figure 
1.1.2  illustrates  the  assignment  of  nodes  to  the  model 
illustrated  in  Figure  1.1.1.  In  this  case  a  node  has  been 
defined  for  each  alteration  of  the  wing  planform  or  for  the 
intersection  of  spars,  ribs,  stringers  or  skins.  Figure  1.1.3 
illustrates  the  finite  element  components  of  the  same  wing 
model  after  the  nodes  of  Figure  1.1.2  have  been  connected 
appropriately. 

The  2-D  elements  (plate  elements)  numbered  1-16  of 
Figure  1.1.3  form  the  upper  skin  while  the  1-D  elements  (bar 
elements)  numbered  129-136  are  stringers.  The  2-D  plate 
elements  33-52  comprise  the  five  ribs  defined  for  this  model 
with  the  rib  caps  (1-D  bar  elements)  being  numbered  65-104. 

The  spars  are  represented  as  2-D  plate  elements  numbered  53-64 
with  the  spar  caps  (1-D  bar  elements)  numbered  105-128. 

Finally,  the  lower  skin  is  comprised  of  2-D  plate  elements 
numbered  17-32  with  stringers  (1-D  bar  elements)  numbered 
137-144.  Vertical  posts  are  defined  to  provide  the  model  with 
structural  integrity  and  to  prevent  collapsing  under  analytical 
conditions.  These  are  numbered  145-159.  These  159  elements 
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Critical 


combine  to  completely  define  a  finite  element  model  suitable 
for  simulated  loading  analysis  by  an  appropriate  finite  element 
model  computer  analysis  program. 

The  process  of  developing  an  accurate  model  capable  of 
efficient  analysis  is  one  that  takes  experience  and  intuition 
to  be  done  properly.  Due  to  the  requirements  of  accuracy  and 
detail  necessary  for  the  construction  of  a  valid  FEM  model, 
computer  programs  have  been  developed  recently  that  alleviate 
the  user  from  the  tedium  of  model  definition  in  cases  where 
the  model  falls  into  a  well-defined  series  of  geometric  shapes. 
The  wing  structure  for  which  this  report  is  directed  does  fall 
inpo  such  a  well-defined  class  of  geometric  structures  and  is, 
therefore,  readily  capable  of  having  models  generated  by 
computer  programs  with  the  engineer  required  to  define  only 
basic  node  points  and  the  geometry  of  the  wing  desired.  A 
computer  program  can  take  this  basic  geometry  data  and  generate 
all  the  nodes,  elements,  model  constraints  and  loading  values 
required.  This  file,  containing  all  the  model  information,  can 
be  analyzed  with  the  appropriate  computer  program.  The  output 
file  from  the  analysis  program  must  be  examined,  generally  with 
a  post-analysis  program  to  extract  information  which  is  desired 
and  present  it  in  a  form,  such  as  graphs,  which  is  readily 
comprehendable.  Figure  1.1.4  illustrates  the  process  of  model 
definition  from  the  engineer's  conceptualization  through  the 
analysis  process.  Figure  1.1.5,  discussed  below,  more  fully 
illustrates  the  process  of  computer  utilization  in  FEM 
development,  analysis  and  postanalysis.  The  various  items 
shown  in  Figure  1.1.5  are  discussed  briefly  here  to  aid  the 
reader  in  forming  a  clearer  general  understanding  of  the 
structural  analysis  tool;  more  detailed  discussions  on  these 
areas  are  contained  in  subsequent  sections  and  chapters  of  this 
report. 

(a)  Input  Data  —  This  item  represents  the  data  which  the 
survivability/vulnerability  engineer  must  supply  to  the 
preprocessor  program  in  order  to  define  completely  the  physical 
structure  to  be  simulated  analytically. 
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Figure  1.1.4a.  Relationship  of  structural  analysis  technique 
with  experimental  technique. 


Figure  1.1.5. 


Finite  Element  Model  Structural  Analysis 
Technique. 


(b)  Preprocessor  —  Assuming  that  the  problem  to  be 
considered  falls  into  one  of  a  number  of  defined  classes  of 
problems,  the  primary  functions  of  the  preprocessor  are  to 
accept  input  data  in  an  abbreviated  format  and  to  generate 
an  expanded  set  of  data  which  is  accepted  by  the  finite 
element  program.  Another  function  of  the  preprocessor 

is  to  display  information  concerning  the  generated  finite 
element  model  to  the  survivability/vulnerability  engineer  in 
the  form  of  printed  output  and  undeformed  geometry  plots. 

(c)  Finite  Element  Program  Data  —  This  item  denotes  the 
data  which  is  generated  and  stored  on  auxiliary  storage  by 

the  preprocessor  and  is  accepted  by  the  finite  element  program. 

(d)  Finite  Element  Program  —  The  finite  element  program 
is  the  backbone  of  the  structural  analysis  technique.  This 
program  takes  the  finite  element  program  data  (load  deck)  and 
performs  analysis  simulating  experimental  conditions. 

(e)  Finite  Element  Program  Output  —  This  item  refers 
to  the  output  which  is  printed  by  the  finite  element  program. 

For  the  current  program,  this  information  is  directed  to 
auxiliary  storage  in  the  form  of  a  WINGMPOST  file  for  use  by  the 
postprocessing  programs  in  addition  to  a  printed  tab  of  the 
data  generated. 

(f)  Postprocessor  —  The  primary  function  of  the 
postprocessor  is  to  edit  the  output  from  the  finite  element 
program  to  eliminate  any  information  which  is  not  pertinent 
to  the  survivability/vulnerability  engineer. 

(g)  Output  Data  —  This  item  represents  the  edited 
output  data  from  the  postprocessor.  The  output  data  available 
to  the  survivability/vulnerability  engineer  is  generated  in  the 
form  of  plotted  displacement  and  stress/strain  results,  deformed 
geometry  plots  and  other  graphs  to  illustrate  the  results  of 
the  analysis  for  quick  and  accurate  interpretation. 
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The  reader  can  see  that  the  pre-  and  postprocessors 
essentially  "straddle"  the  finite  element  program  in  such  a 
manner  that  the  composite  analysis  tool  is  convenient  for  the 
survivability/vulnerability  engineer  to  use  effectively. 

The  reader  may  feel  that  these  programs  and  their 
interaction  are  very  complex  and  provide  very  little  effective 
results.  However,  once  familiar  with  the  processes  involved 
and  the  mechanics  explained  in  this  report  to  model  and 
analyze  wing  specimens,  the  user  will  find  that  considerable 
savings  in  time,  effort,  and  money  will  be  realized  apart 
from  the  savings  in  materials  once  these  analytical  procedures 
are  able  to  fully  supplement  regular  structural  testing. 


1.2  PREPROCESSOR  PROGRAM  WINGEN 

The  nature  of  finite  element  modeling  is  such  that  much  of 
the  tedious  work  of  generating  a  model  may  be  performed  very 
effectively  by  a  computer.  The  preprocessor  program  WINGEN  was 
designed  specifically  to  alleviate  the  tedium  of  finite  element 
modeling  by  providing  for  the  generation  of  nodes  and  element 
connectivities  for  several  classes  of  wing  structures.  The 
program  will  accept  as  input  an  abbreviated  data  file  containing 
the  basic  defining  parameters  for  the  class  of  wing  the  user 
wishes  to  develop  and  will  generate  all  the  nodes  and  finite 
elements  necessary  to  fully  define  the  model.  Since  the 
abbreviated  data  file  is  generally  no  more  than  three  or  four 
dozen  data  cards  and  the  output  file  generated  can  contain,  at 
the  least  several  hundred  and  very  easily  several  thousand 
cards,  the  savings  in  time  and  energy  is  immediatley  apparent, 
especially  when  alterations  to  damage  specifications  and 
loading  conditions  are  made  to  the  model.  WINGEN  has  been 
developed  to  accommodate  a  batch  format  (card-image)  data  deck 
or  can  be  used  in  a  fully  interactive  mode.  The  same  input 
information  is  required  for  either  case  and  the  output  will 
be  identical.  If  the  user  selects  to  utilize  the  interactive 
input  method  (no  cards  required)  the  program  will  create  a 
model  definition  data  deck  for  future  use  with  the  program  to 
prevent  the  user  from  having  to  interactively  input  the  basic 
model  data  more  than  once.  Alterations  and  damages  to  the 
model  or  new  loads  may  be  specified  in  the  batch  format  (model 
definition)  data  deck  and  the  finite  element  data  deck  can  be 
regenerated  with  these  changes  very  quickly.  Figure  1.2.1 
illustrates  how  WINGEN  operates  with  the  various  files  and 
user  interaction. 

The  flow  of  information  diagram  (Figure  1.2.1)  for  WINGEN 
illustrates  the  various  input  and  output  features  of  the  program. 
Two  input  options  are  available: 


* 


(a)  Card  Deck  input  for  batch  mode  execution.  This  option 
requires  the  user  to  place  the  model  data  on  cards 

in  the  card  deck  or  to  attach  a  permanent  computer 
resident  disk  file  as  a  local  file,  TAPE4 ,  which 
contains  the  model  definition  data  detailed  in 
Section  2.8. 

(b)  Interactive  input  allows  the  user  to  select  fully 
interactive  input  where  all  model  definition  data 
is  entered  directly  through  a  terminal  (causing  a 
USRDATAFILE  to  be  created)  or  semi-interactive  input 
where  all  model  definition  data  is  placed  on  a  local 
file  (TAPE4)  for  the  program  to  use. 

Several  output  options  are  available  from  WINGEN: 

(a)  A  Batch  mode  output  listing  is  generated  if  the 
card  deck  input  option  is  selected. 

(b)  Model  Geometry  Plots  are  generated  on  an  interactive 
graphics  device,  if  requested. 

(c)  A  Model  Definition  Data  file  is  created  and  made  a 
permanent  disk  file  with  the  name  USRDATAFILE  if 
the  user  selects  the  fully  interactive  input  option. 

(d)  A  Model  Load  Deck  file  which  is  used  to  initiate 
an  analysis  of  the  FEM  by  MAGNA  is  generated  if 
requested  by  the  user.  This  file,  TAPEll,  must  be 
catalogued  by  the  user  if  a  permanent  disk  file  is 
desired. 

WINGEN  has  been  designed  primarily  to  accommodate  a  user 
who  wishes  to  execute  the  program  interactively  to  generate  a 
USRDATAFILE.  Any  modifications  made  to  the  Model  Definition 
Data  stored  on  the  USRDATAFILE  can  be  made  and  the  USRDATAFILE 
attached  as  a  local  file  (TAPE4 )  to  be  used  for  a  semi-interactive 
execution  of  WINGEN.  This  means  the  user  should  not  need  to 
utilize  the  full  interactive  execution  option  more  than  once  for 
each  model.  Section  2  discusses  WINGEN  in  detail. 
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Figure  1.2.1.  WINGEN  operating  environment.  I  =  interactive 

mode;  SI  =  semi-interactive  mode;  B  =  batch  mode. 
These  operating  modes  are  explained  in  Chapter  2. 
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1.3  ANALYSIS  PROGRAM  MAGNA 

Once  a  finite  element  model  is  defined  it  may  be  analyzed 
utilizing  a  finite  element  analysis  program.  MAGNA  is  a 
special  purpose  computer  program  developed  specifically  for  the 
analysis  of  large  displacements  of  finite  element  structures 
with  up  to  several  thousand  degrees  of  freedom.  It  can  handle 
material  nonlinearities  as  well  as  geometric  nonlinearities. 
WINGEN  provides  the  user  with  all  the  necessary  parameters  to 
invoke  the  appropriate  analysis  types  when  the  finite  element 
model  load  deck  is  created.  Since  MAGNA  has  a  wide  variety  of 
applications  this  allows  the  user  to  utilize  MAGNA  without 
the  initial  delay  while  becoming  highly  familiar  with  it.  MAGNA 
will  perform  the  analysis  of  the  model  utilizing  the  data  file 
supplied  by  WINGEN.  This  will  result  in  the  production  of  an 
output  or  -MPOST  file  (WINGMPOST)  which  can  be  utilized  by  the 
postprocessing  programs  PLOTBOB  and  CONTOUR.  Figure  1.3.1 
illustrates  the  MAGNA  operating  environment. 

This  flow  of  information  figure  illustrates  how  the  Load 
Deck  created  by  WINGEN  as  TAPE11  can  be  used  as  input  to  execute 
a  MAGNA  analysis.  TAPE11  can  be  either  saved  as  a  permanent 
disk  file  and  accessed  via  card  input  or  through  interactive 
terminals  or  it  can  be  punched  as  a  card  deck  and  submitted 
through  card  reader  to  initiate  a  MAGNA  analysis.  Since  MAGNA 
must  execute  as  a  Batch  job  there  can  be  no  user  interaction 
with  the  analysis.  There  is  capability  for  two  output  files 
from  MAGNA.  An  execution  listing  of  the  analysis  containing 
results  in  a  readable  format  is  always  generated.  In  addition 
the  user  will  find  a  permanent  disk  file  with  the  name  WINGMPOST 
containing  the  same  information  as  is  contained  on  the  output 
listing  but  in  a  condensed  format  for  use  by  the  postprocessing 
programs  CONTOUR  and  PLOTBOB  (see  Section  4).  Reference  1 
contains  complete  information  on  MAGNA  and  how  to  use  it. 

Section  3  provides  an  overview  of  MAGNA  features  and  its 
interaction  with  the  pre-  and  post-processors  as  well  as  a 
sample  output  listing. 
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.3.1.  MAGNA  operating  environment.  MAGNA  can  only  be 
executed  in  a  batch  mode.  Refer  to  Chapter  3 
for  further  information. 


1.4  POSTPROCESSOR  PROGRAMS  CONTOUR  AND  PLOTBOB 

Postprocessing  becomes  an  important  tool  in  finite  element 
modeling  and  analysis  due  to  the  large  volumes  of  output 
necessarily  generated.  Postprocessing  capability  described  in 
this  report  relates  primarily  to  graphical  representation  of 
MAGNA  analysis  results.  The  results  are  displayed  as  either 
deformed  structure  plots  or  as  contour  plots  with  stress  or 
strain  values  illustrated  by  contour  lines  drawn  on  the  model. 
This  type  of  data  presentation  is  very  effective  in  allowing 
the  user  to  immediately  comprehend  the  distribution  of  forces 
throughout  a  structure  and  thereby  help  to  more  rapidly  isolate 
where  and  when  the  structure  might  fail.  Since  graphical 
presentation  of  data  is  considerably  easier  to  understand  them 
large  computer  output  tabs  and  because  much  data  presented  in 
output  tabs  is  useless  unless  the  user  can  know  precisely  how 
the  data  relates  to  the  model  in  question,  plotted  results 
answer  a  critical  need  by  bridging  the  theoretical  and  the 
conceptual  ideas  of  finite  element  modeling. 

The  interaction  of  WINGEN,  MAGNA,  PLOTBOB,  CONTOUR,  and 
the  user  is  illustrated  in  Figure  1.4.1.  These  four  programs 
provide  a  complete  set  of  tools  for  the  model  development, 
analysis,  and  post-analysis  data  reduction  required  in  finite 
element  model  analysis. 

The  flow  of  information  illustrated  in  Figure  1.4.1  is  a 
combining  of  the  information  from  Figure  1.2.1  and  1.3.1  with 
the  addition  of  the  postprocessing  program. s  It  becomes 
apparent  from  this  figure  that  PLOTBOB  has  the  capability  of 
making  geometry  plots  illustrating  the  structure  to  be  analyzed 
from  the  same  Load  Deck  data  file  that  initiates  a  MAGNA 
analysis.  PLOTBOB  will  also  generate  geometry  plots  of  the 
WINGMPOST  file  that  is  produced  by  a  MAGNA  analysis.  CONTOUR 
will  generate  contour,  displacement  and  geometry  plots  of  the 
WINGMPOST  data  file  only.  Section  4  details  use  of  these  two 
programs.  In  addition  Reference  4  contains  some  additional 
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information  about  these  two  postprocessing  programs.  In  Reference 
4  CONTOUR  has  been  renamed  to  CPLOT  and  PLOTBOB  has  the  new 
name  GPLOT. 
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Figure  1.4.1. 


MODEL  development  and  analysis  for  wing  damage 
flight  loads  simulation.  This  illustrates  the 
interaction  of  WINGEN,  MAGNA,  PLOTBOB,  and  CONTOUR 
programs  in  model  development,  analysis  and 
postprocessing. 


SECTION  2 

FINITE  ELEMENT  MODEL  PREPROCESSOR 

2 . 1  INTRODUCTION 

WINGEN  is  a  convenient  and  flexible  preprocessor  computer 
program  designed  specifically  for  the  generation  of  finite 
element  nodes  and  connectivities  input  data  for  damaged  and 
undamaged  wing  models  for  the  finite  element  analysis  program 
MAGNA.  The  program  can  be  run  in  either  the  interactive  or 
batch  modes.  Interactively  the  model  designer  can  initiate  the 
program  and  generate  an  entire  data  file  with  simple  geometries 
and  basic  wing  planform  dimensions.  For  batch  mode  processing 
an  input  data  file  containing  this  same  basic  information  is 
required. 

In  either  case  of  utilizing  batch  or  interactive  mode  a 
load  deck  is  created  which  fully  defines  the  geometry  of  any 
one  of  four  classes  of  wings  (see  MODEL  DEFINITION)  with 
necessary  linear  constraints,  load  conditions,  finite  element 
boundary  conditions,  damage  specifications,  and  model  mesh 
refinements  over  selected  areas.  In  addition,  plots 
illustrating  the  basic  geometry  of  the  wing  with  labeling  of 
the  elements  and/or  nodes  may  be  generated  interactively  on  a 
graphics  terminal  or  chart  plotter.  Capability  is  currently 
present  that  generates  8-node  shell  elements  for  upper  and  lower 
wing  skins  and  necessary  constraints  for  analysis  of  a  refined 
mesh  model.  Figure  2.1.1  illustrates  the  position  WINGEN  holds 
in  the  process  of  wing  model  development. 


2.2  WINGEN 


PROGRAM  ORGANIZATION 


The  WINGEN  preprocessor  program  is  designed  to  be  executed 
in  a  batch,  interactive  or  semi-interactive  mode.  Batch  mode 
use  is  completely  non-interactive  and  is  initiated  with  cards 
or  a  card-image  file.  Figure  2.2.1  illustrates  a  sample  batch 
input  deck.  The  important  process  to  be  aware  of  in  batch  job 
execution  is  that  all  data  must  be  established  prior  to 
executing  the  program  and  be  placed  on  a  permanent  disk  file 
accessible  to  the  job  or  on  cards  (Figure  2.2.2).  The  user  must 
catalog  the  output  load  deck  file  for  future  use.  If  plots  of 
the  basic  structure  are  desired,  program  PLOTBOB  may  be 
executed  with  the  MAGNA  input  file  (load  deck) .  The  advantage 
to  using  batch  mode  is  that  no  user  interaction  is  required. 

Semi-interactive  mode  use  is  accomplished  by  using  a 
data  file  established  prior  to  the  execution  of  the  program 
with  the  information  defining  the  wing  formatted  as  in  batch 
mode  input.  In  this  case,  however,  the  program  control 
directives  are  displayed  on  the  computer  terminal  and  you  must 
answer  them  as  they  appear.  In  addition  a  number  of  questions 
concerning  the  type  of  analysis  to  be  done  will  be  requested. 

The  advantage  to  this  type  of  execution  is  that  it  requires 
considerably  less  time  for  the  operator  to  get  a  load  deck 
than  a  fully  interactive  run  and  one  may  generate  plots  of  the 
structure  without  needing  to  execute  PLOTBOB  with  the  output 
load  deck  file,  as  is  required  when  using  batch  mode. 

Full,  interactive  use  of  WINGEN  requires  the  individual 
to  have  before  him  all  pertinent  data  to  define  the  wing.  The 
program  will  ask  all  questions  and  the  user  must  answer  with 
the  appropriate  piece  of  information.  The  program  will  create 
and  catalog  (make  a  permanent  disk  file  of)  the  data  input  by 
the  operator  in  model  definition  data  format.  This  file  is 
given  the  name  'USRDATAFILE'  and  will  be  cataloged  under  the 
ID  which  the  user  logged  into  the  computer  with.  It  must  be 
noted  that  on  the  CDC  there  is  a  limit  of  five  permanent  files 
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XXXX , T200 , 10600 , CM70000 , STANY .  DXXXXXX , DOE , BLDG4 5,2556666. 
ATTACH, F,WINGEN 
ATTACH , TAPE 5 , WINGDATA . 
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Figure  2.2.1.  A  sample  batch  input  data  deck.  File  WINGDATA 
must  contain  the  wing  data  discussed  in 
Section  2.8. 
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Job  Control  Cards 


A.  Example  of  batch  job  deck  for  use  of  a  permanent  file 
WINGEN  input  model  deck  (batch  deck) .  This  deck  is 
listed  in  Figure  2.2.1. 


B.  Example  of  batch  job  deck  for  use  of  a  card  deck  of 
WINGEN  input  model  data.  This  deck  is  listed  in 
Figure  2.2.3. 


Figure  2.2.2.  Two  Examples  of  Batch  Execution  for  WINGEN. 


w 


XXXX ,T2 00, 10600, CM70000 , STAN Y .  DXXXXXX , DOE , BLDG4  5,2556666. 
ATTACH, F,WINGEN. 

COPYCR , INPUT , TAPE 5 . 

REWIND, TAPE5. 

REQUEST, TAPE 11 , *PF . 

F. 

CATALOG, TAPE11,WINGINPUT,RP=300. 

7/8/9  (EOR) 


Pf. PLICA  TEST  SPECIMEN 


PROFILE  PI  Cl 

SI 

PLhN  S3. 

DEPTH 

S8  74.  S 

17. 7S 

SKIN  1 

RIBS  S 

1  .25000 

.25000 

a  i 

1  0.00 

9.30 

53.88 

9.30 

.1000 

.1500 

2  1 

1  0.00 

19.60 

53.88 

19.60 

.1000 

.1500 

2  1 

1  0.00 

31.20 

53.88 

31.20 

.1090 

.1500 

2  1 

1  0.00 

43.90 

53.88 

43.90 

.1000 

.1500 

2  1 

1  0.00 

59.69 

53.88 

59.60 

.1000 

.1500 

SPARS  4 

2  1 

1  0.00 

0.00 

0.00 

74.50 

.1875 

2.6250 

2  1 

1  17.97 

0.00 

17.97 

74.50 

.1250 

1.8750 

2  1 

1  35.91 

0.00 

35.91 

74.50 

.1250 

1 .8750 

2  1 

1  53.88 

0.00 

53.88 

74.50 

.1875 

2-6250 

MODIFY 

POSTS  -1 

2  0 

0.0000 

.  3Q00 

REFINE  1 

THICK 

DAMAGE 

l 

7/8/9 

(EOR) 

N 


Y 
N 

Y 

Y 
N 

XXXX 

DXXXXXX , DOE , BLDG4 5,255-6666. 

1 

1 

1 

1 

7/8/9  (EOR) 

6/7/8/9  (EOF) 


Figure  2.2.3.  Sample  Batch  Input  Deck  for  Figure  2.2.2b 


(cycles)  with  the  same  name,  so  it  is  suggested  that  the  user 
rename  the  file  created  to  be  something  more  representative 
of  the  project  being  done.  (If  the  program  aborts  the  user  may 
examine  local  tape  file  7  (TAPE7 )  for  correct  data  and  catalog 
it  as  a  WINGEN  input  file  for  batch  or  semi-interactive  use. 

This  may  save  another  lengthy  session  of  interactive  use  of 
the  program.)  The  ' USRDATAFILE '  generated  can  then  be  used 
as  input  to  WINGEN  to  create  a  load  deck  for  input  for  a  MAGNA 
analysis  run.  Figure  1.2.1  illustrates  a  flow  chart  of  the 
preprocessor  operation.  Note  that  a  'USRDATAFILE',  which 
actually  formats  and  echoes  the  input,  a  MAGNA  input  deck 
(TAPEll )  and  plots  can  all  be  generated  in  a  single  execution 
of  WINGEN. 

During  execution,  the  program  has  four  processing  phases: 
Control  Directive  Input;  Node  and  Element  Generation;  MAGNA 
Input  Deck  Creation;  and  Plotting. 

Control  Directive  Input  is  requested  by  the  program  to 
determine  if  the  run  is  interactive  or  batch,  whether  or  not  to 
create  a  MAGNA  input  deck,  if  shell  elements  should  be  used  for 
the  upper  and  lower  skins,  etc.  The  input  will  look  as  follows: 

INTERACTIVE  RUN?  (Y,N) . : 

CREATE  A  LOAD  DECK?  (Y ,N) . . . : 

GENERATE  A  GRAPH?  (Y,N) . : 

LIST  NODES  AND  COORDINATES?  (Y,N)  . . . .  : 

LIST  ELEMENT  CONNECTIVITIES?  (Y,N) . . . : 

USE  SHELL  ELEMENTS?  (Y,N) . : 

Once  all  such  questions  have  been  answered  the  word  'START' 
will  appear  and  the  program  will  commence  with  the  generation 
of  nodes  and  element  connectivities. 

Nodes  and  Element  Connectivities  are  generated  as  the 
primary  function  of  WINGEN.  The  user  may  request  to  have  an 
output  listing  of  this  information  or  may  just  wish  to  see  the 
MAGNA  input  deck  that  is  created  at  the  conclusion  of  the  run. 


For  complex  models  the  output  listing  can  be  very  lengthy  and 
the  user  may  find  it  unnecessary.  Options  are  provided 
to  allow  suppression  of  this  output.  This  phase  of  WINGEN 
is  also  responsible  for  making  modifications  to  the  basic 
model  and  implementing  damage  specifications,  refinement 
directives  and  necessary  linear  constraints.  It  is  for  this 
phase  that  the  model  definition  card  input  data  file  is  most 
useful.  (Refer  to  Sections  2.7  and  2.8.) 

MAGNA  Input  Data  Deck  creation  follows  the  defining 
of  the  nodes  and  elements.  The  program  requests  several  pieces 
of  information  concerning  the  type  of  analysis,  solution 
parameters,  types  and  magnitudes  of  loads  applied  and  other 
necessary  information  concerning  the  execution  of  the  analysis 
run.  If  the  upper  and  lower  skin  meshes  of  the  wing  model  are 
refined  it  is  necessary  for  the  user  to  request  the  use  of 
three-dimensional  shell  elements  in  lieu  of  the  two-dimensional 
membrane  elements  used  for  the  basic  model.  Failure  to  do  this 
will  result  in  problems  arising  during  the  analysis.  WINGEN 
will  create  a  card-image  deck  containing  all  the  information 
necessary  to  initiate  a  MAGNA  analysis  run  including  all  job 
control  cards.  This  data  will  be  placed  on  local  file  TAPEll 
which  the  user  must  catalog  if  he  wishes  to  save  it.  Refer  to 
Section  2.9  'WINGEN  -  LOAD  DECK'  for  further  information. 

Plotting  capabilities  of  WINGEN  are  designed  to  provide 
the  user  with  a  basic  definition  of  the  model  just  generated. 
The  user  will  find  the  program  prompting  for  eye  position, 
node  and/or  element  labeling  and  type  of  view  (orthogonal  or 
perspective) .  These  capabilities  provide  the  necessary 
graphics  capabilities  for  this  type  of  model  development.  The 
program  will  plot  all  the  elements  each  time,  labeling  the 
nodes  or  elements  as  requested.  This  allows  rapid  visual 
verification  of  the  model  geometry  and  the  element  and  node 
generation.  For  further  reference  see  Section  2.10  'PLOTTING 
CAPABILITIES'.  Should  the  user  desire  more  selective  plotting 
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and  labeling  of  a  model  (e.g.,  for  a  very  complex  model)  he 
may  take  the  MAGNA  input  data  deck  created  by  this  program 
(on  TAPEll )  and  execute  the  PLOTBOB  plotting  program.  PLOTBOB 
will  provide  the  user  with  full  interactive  graphics 
capabilities  including  selective  element  plotting,  zoom  on 
model  parts,  labeling  of  components,  etc.  For  more  information 
on  PLOTBOB  see  Section  4.2  -  'INTRODUCTION  TO  PLOTBOB'. 

Should  the  user  select  shell  elements  for  the  upper  and 
lower  skins  and  wish  to  see  plots  of  the  finished  model  he 
will  have  to  execute  PLOTBOB  .  WINGEN  plots  done  with  the 
shell  elements  selected  will  only  illustrate  the  basic  model 
with  membrane  elements.  This  will  allow  the  user  an  opportunity 
to  verify  the  model  before  the  more  complicated  shell  elements 
are  generated.  The  shell  elements  add  such  a  degree  of 
complexity  to  model  plotting  that  the  requirements  for  plotting 
can  best  be  met  with  the  PLOTBOB  program. 
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2 . 3  WINGEN  -  MODEL  DEFINITION 


WINGEN  is  a  special-purpose  preprocessor  for  the  generation 
of  necessary  geometric  components  of  typical  wings  in  a  format 
acceptable  to  the  general  purpose  finite  element  analysis 
program  MAGNA.  Certain  common  components  have  been  identified 
and  linked  together  to  allow  the  rapid  generation  of  a  model 
suitable  for  efficient  and  accurate  analysis.  The  basic 
components  of  the  finite  element  model  (fern)  are  two-dimensional 
plane  stress  (membrane)  and  shear  panel  elements,  one-dimensional 
truss  elements,  and  three-dimensional  shell  elements  (MAGNA 
element  types  3,  4  and  5,  respectively). 

The  plane  stress  element  type  (element  subtype  1  of  MAGNA 
element  type  3)  is  used  for  the  upper  and  lower  skins  of  the 
wing  model.  These  must  be  converted  to  shell  elements  (MAGNA 
element  type  5)  if  the  skin  areas  are  refined  into  a  smaller 
mesh  as  there  will  be  no  constraints  or  other  element 
connectivities  to  prevent  the  nodes  from  achieving  an  undesired 
infinite  displacement  during  analysis.  The  3-D  nature  of 
shell  elements  prevents  this  from  occurring.  The  shear  panel 
element  type  (element  subtype  3  of  MAGNA  element  type  3)  is 
used  for  definition  of  the  ribs  and  spars.  Upon  refinement 
necessary  linear  constraints  are  generated  to  prevent  unreal 
displacements  of  newly  generated  nodes.  All  type  3  elements 
require  an  element  number,  an  element  subtype  (1  or  3) ,  a 
material  property  code  (1  or  2),  four  nodes  and  a  thickness 
for  complete  definition.  The  program  predefines  two  material 
property  codes  as  these  are  the  most  commonly  used  materials 
in  wing  construction. 


Material 


Elastic  Modulus  lb/in^ 


Poisson's  Ratio 


Mass  Density  ib  sec2/in4 


EQ  Stress  at  1st  Yield 
lb/in^ 


j  Strain  Hardening  Curve 


|  Uniaxial  Stress/Strain 
Data  Curve 


LINEAR  AND  NONLINEAR  ANALYSIS 

Material  1 

] 

Material  2 

aluminum 

steel 

.10*108 

.30*10* 

.  30 

.33 

.  259*1 0~3 

.  725*10~3 

.io*io5 

. 30*105 

NONLINEAR  / 

ANALYSIS  ONLY 

The  nodal  connectivity  will  be  generated  by  the  program 
for  each  element,  starting  at  one  corner  and  working  completely 
around : 


The  program  will  automatically  generate  all  nodes  required  to 
define  the  geometry  and  connect  these  nodes  to  produce  the 
proper  elements  for  the  model  definitions.  As  finite  element 
modeling  may  be  somewhat  restrictive  a  'MODIFY'  command  has 
been  included  to  alter  the  basic  structure  generated  to 
produce  desired  non-uniform  changes.  The  thickness  parameter 
requested  for  the  skins  is  the  average  thickness  over  the  entire 
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skin.  If  thickness  of  the  wing  skin  changes  dramatically  it 
may  be  necessary  to  alter  the  thickness  value  for  the 
appropriate  elements  on  the  MAGNA  input  file  after  the 
preprocessor  run  is  completed.  Please  note  that  when  a  model 
is  refined  in  the  upper  and  lower  skins,  shell  elements  are  used 
to  redefine  the  skins.  In  the  use  of  shell  elements  new  nodes 
must  be  generated.  To  accomplish  the  generation  of  new  nodes 
the  program  utilizes  the  upper  and  lower  skin  thickness  values 
and  generates  the  new  nodes  at  that  displacement  from  the  old 
nodes,  yielding  a  three-dimensional  element  with  the  third 
dimension  being  the  skin  thickness.  If  the  skin  thicknesses 
vary  considerably  it  is  advised  that  the  model  be  generated  and 
the  appropriate  component  (usually  the  z-component)  be  altered 
as  needed  in  the  nodal  coordinates  list  of  the  WINGEN  load  deck 
output  data  file  created  to  accommodate  the  variations. 

The  truss  element  type  (MAGNA  element  type  4)  is  used  to 
define  the  rib  and  spar  caps  and  the  vertical  posts.  The  truss 
element  is  a  one-dimensional  element  requiring  two  nodes  for  its 
definition  along  with  a  material  property  code  and  the  cross- 
sectional  area  of  the  structure  being  represented.  Listed  below 
are  the  two  material  sets  predefined  by  the  program  for 
aluminum  and  steel : 

LINEAR  AND  NONLINEAR  ANALYSIS 


Material 

— 

Elastic  Modulus  lb/in 

Mass  Density  2  4 

lb  sec  /in4 

EQ  Stress  at  1st  Yield 
lb/in2 


Material  1 
aluminum 


.10*10c 


,259*10 


.10*10' 


Material  2 
steel 

. 30*108 
.725*10~3 


.30*10' 
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NONLINEAR  ANALYSIS  ONLY 


Strain  Hardening 
Curve 


Uniaxial  Stress/ 
Strain  Data  Curve 


As  with  other  element  types,  additional  material  property  codes 
may  be  added  to  the  WINGEN  load  deck  output  file  for  additional 
material  used  in  the  wing.  For  information  on  how  to  make  this 
and  other  modifications  to  the  load  deck  output  file  see 
Section  2.9.  The  truss  element  requires  only  two  nodes  for 
its  spatial  definition: 


The  cross-sectional  area  of  the  element  is  obtained  by  computing 
the  cross-sectional  area  of  the  component  being  represented  by 
the  truss  element.  This  is  illustrated  in  Figure  2.3.1. 


Spar  Cap 


!  * - '> 


Spar 


V 


14 —  X 


Lower  Wing 
i  Skin 


Figure  2.3.1.  The  cross-sectional  area  of  a  truss  element  (spar 

and  rib  caps) :  (y*w)  +  (x*z)  =  area  of  cross-section. 
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The  3-D  thin  shell  element  type  (MAGNA  element  type  5)  is 
optionally  generated  to  provide  a  stable  model  where  the 
basic  model  has  been  refined.  This  is  necessary  when  any  new 
skin  nodes  generated  during  refinement  are  not  simultaneously 
incorporated  into  rib,  spar  and  skin  elements.  If  shell 
elements  or  linear  constraints  are  not  utilized  in  this  special 
case  then  the  newly  generated  nodes  will  undergo  infinite 
displacements  during  the  analysis,  which  will  yield  undesirable 
results.  The  definition  of  the  3-D  thin  shell  element  requires 
the  following  information:  a)  an  element  number;  b)  a  material 
property  code;  and  c)  eight  nodes.  The  element  numbers  are 
supplied  by  the  program.  The  program  automatically  substitutes 
the  previous  material  codes  established  for  the  original  2-D 
membrane  element  for  those  codes  required  in  the  3-D  element  type: 


LINEAR  AND  NONLINEAR  ANALYSIS 

Material  1 

Material  2 

Material 

aluminum 

steel 

2 

Elastic  Modulus  lb/in 

.10*108 

.30*108 

Poisson's  Ratio 

.30 

.33 

2  4 

Mass  Density  lb  sec  /in 

.259*10~3 

.  725*10~3 

Additional  material  properties  may  be  added  to  the  WINGEN  output 
file.  Further  information  on  this  and  other  changes  to  the  load 
deck  output  file  may  be  referenced  in  Section  2.9  -  'WINGEN  - 
LOAD  DECK ' . 

The  nodal  connectivity  for  the  thin  shell  element  is 
illustrated  below: 
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An  additional  advantage  with  the  3-D  thin  shell  element  is 
the  lack  of  a  need  to  maintain  an  even  thickness  throughout  the 
element.  This  is  the  case  in  the  2-D  plane  stress  element  type 
where  there  is  a  uniform  thickness  across  each  individual 
element.  The  3-D  shell  element  allows  the  user  the  capability 
of  altering  the  nodal  coordinates  to  arrange  for  any  systematic 
changes  in  skin  thicknesses. 

WINGEN  has  been  designed  to  accommodate  three  common 
profile  classes  of  wings:  rectangular,  full  swept  and  piecewise 
linear  (see  Figure  2.3.2)  as  well  as  a  general  class  where  the 
user  must  identify  all  nodal  coordinates  (see  P4-C4-S4  profile 
type  in  Section  2.8  'WINGEN  -  MODEL  DEFINITION  DATA  FORMAT').  Ir 
the  three  common  profile  classes  (P1-P3)  the  preprocessor 
program  will  generate  all  nodes  required  to  define  the  ribs, 
spars  and  skins  and  derive  all  necessary  nodal  connectivities 
for  the  basic  wing  model  elements.  The  general  class  (P4)  model 
will  have  all  elements  and  connectivities  generated  once  the 
nodal  coordinates  have  been  supplied.  The  nodal  coordinates 
are  produced  by  starting  at  the  initially  given  x  and  y  values 
and  alternating  with  the  generation  of  a  lower  skin  node  first 
then  an  upper  skin  node  working  chordwise  fore  to  aft,  moving 
to  the  next  rib  (in  a  spanwise  direction)  and  again  working 
chordwise  from  fore  to  aft  alternating  between  lower  and  upper 
skins  generating  all  nodes  for  the  basic  model. 

The  program  then  processes  any  refinement  directives, 
adding  nodes  as  required.  In  the  refinement  process  the 
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originally  numbered  nodes  are  resequenced  to  allow  the  newly 
generated  nodes  to  be  incorporated  sequentially.  This  provides 
for  a  more  comprehensible  model  when  all  modifications  are 
concluded.  The  refinements  are  generated  much  the  same  way  as 
the  model  was  originally  defined,  with  refinement  specifications 
being  acted  on  one  at  a  time.  The  refinement  nodes  are  added 
by  dividing  the  difference  in  the  associated  coordinates  by 
the  number  of  refinements  to  be  added  over  that  difference. 

This  gives  an  incremental  value  which  is  added  to  the  lower 
value  to  successively  generate  the  new  nodes.  The  refine 
command  will  always  completely  traverse  the  chordwise,  spanwise 
or  depthwise  bays  being  refined  to  ensure  element  compatibility. 
(See  Figure  2.8.14). 

Once  the  refine  directive  has  been  executed  all  the  nodal 
coordinates  will  have  been  defined.  The  program  then  processes 
the  damage/modify  directives  while  generating  the  element 
connectivities.  The  damaoe/modi f y  directives  essentially  tell 
the  program  where  to  omit  elements  not  relevant  to  the  model. 

The  program  establishes  a  table  of  damage  and  modification 
parameters  and  searches  it  to  determine  if  the  current  element 
under  consideration  should  be  eliminated.  If  the  element  meets 
the  criteria  established  for  removal  it  is  not  included  in  the 
model  element  list.  The  modify  directive  is  a  special  case  of 
the  damaqe  directive  for  use  when  the  model  must  be  altered 
after  using  the  FEM  generator  to  actually  achieve  a  true 
representation  of  the  wing.  Figure  2.3.3  illustrates  this  type 
of  situation.  There  is  no  processing  change  in  how  modify  and 
damage  directives  are  handled  by  the  program. 

The  generation  of  tne  elements  follows  the  same  format  as 
the  generation  of  the  nodes.  The  skin  elements  are  generated 
first  starting  with  the  lower  skin  and  alternating  between  lower 
and  upper  skins  chordwise  (fore  to  aft)  then  advancing  to  the 
fore  location  on  the  next  spanwise  (root  to  tip)  rib  and 
continuing  until  all  skin  elements  have  been  defined. 
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Following  the  skin  element  definition  all  the  ribs  are 
next  defined  starting  at  the  root  end  and  proceeding  fore  to 
aft  chordwise  then  moving  to  the  next  rib  spanwise  (root  to  tip) 
until  all  ribs  have  been  defined.  The  spar  elements  are  then 
generated  starting  at  the  foremost  rib  and  working  spanwise 
(from  root  to  tip)  then  moving  chordwise  (fore  to  aft)  to  the 
next  spar  and  continuing  until  all  spars  have  been  defined. 

Once  all  type  3  elements  have  been  defined  the  program 
progresses  to  identifying  the  type  4  elements.  This  process 
parallels  the  defining  of  the  element  type  3  ribs  and  spars 
described  above:  first  the  rib  caps  are  defined  starting  on 
the  lower  skin  and  alternating  between  lower  and  upper  skins, 
moving  fore  to  aft  and  root  to  tip  (chordwise  and  spanwise) . 

The  spar  caps  definition  begin  on  the  lower  skin  and  alternate 
between  lower  and  upper  skins;  moving  root  to  tip  and  fore  to 
aft  (spanwise  and  chordwise).  The  last  component  to  be  defined 
are  the  posts.  Posts  are  vertical  elements  included  to  give 
the  structure  stability  and  prevent  collapsing  of  the  structure 
under  analytical  conditions.  These  are  defined  originating  at 
the  wing  coordinate  origin  and  proceeding  chordwise  then 
spanwise  until  all  posts  have  been  defined. 
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2.4  PROCEDURE  FOR  EXECUTING  WINGEN 

To  utilize  WINGEN  one  must  follow  the  procedure 

a)  LOGIN... 

b)  ATTACH, TAPE4 ,  (Model  Definition  data  file)  ] 

c)  ATTACH , F , WINGEN 

d)  REWIND, F[,TAPE4] 

e)  ATTACH , LIB, PLOT3D, ID=KING,SN=AFFDL  J 

f)  ATTACH , LIB1 , TEKLIB , ID=LIBRARY , SN=ASD  ] 

or  ATTACH , LIB1 , HPPLOT21 , ID=KING , SN=AFFDL  ] 

g)  LIBRARY, LIB, LIB1  ] 

h)  F. 

The  above  procedure  executes  a  segment  loaded  binary  version 
of  WINGEN.  TAPE4  is  required  for  a  batch  or  semi-interactive 
execution  of  the  program.  Data  formats  and  definitions  for  TAPE4 
are  given  in  Sections  2.3  -  MODEL  DEFINITION  and  2.8  -  MODEL 
DEFINITION  DATA  FORMAT.  If  no  plots  of  the  model  geometry  are 
desired  then  steps  e-g  above  do  not  have  to  be  executed.  These 
steps  serve  to  attach  the  plotting  routines.  Step  f  above  must 
have  the  user  attach  the  proper  library  of  plotting  routines 
dependent  on  the  plotting  device  being  utilized.  If  plots  are  to 
be  done  on  a  Tektronix  graphics  terminal  or  PLOT10  emulator  then 
TEKLIB  must  be  attached,  otherwise  the  user  is  constrained  to 
plotting  on  the  Hewlett-Packard  7221  plotter  and  must  attach 
HPPLOT21  instead  of  TEKLIB.  Refer  to  Section  2.10  -  WINGEN 
PLOTTING  CAPABILITIES  for  more  information  on  plotting.  Appendix 
A  illustrates  in  detail  how  to  execute  WINGEN.  An  alternative 
method  is  as  follows: 

a)  LOGIN,.... 

b)  ATTACH, TAPE4 ,USRDATAFILE 

c)  ATTACH , PROC , WINGEN , ID= BRUNER , SN=AFFDL 

d)  BEGIN, WINGEN, PROC.  (for  Tektronix  plots) 

or  BEGIN, WINGEN, PR0C,H,  (for  Hewlett-Packard  plots) 

This  procedure  will  execute  WINGEN  and  catalog  TAPEll  as  WINGINPUT. 


listed  below: 


not  required  for 
interactive  mode 


only  required 
if  model 
geometry 
plots  are 
desired 
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2.5  WINGEN  OPERATION  -  PROGRAM  INITIATION 

Once  the  program  WINGEN  has  been  initiated  the  user  will 
be  prompted  by  several  questions  regarding  the  nature  of  the 
processing : 

INTERACTIVE  RUN?  (Y,N) . : 

CREATE  A  LOAD  DECK?  (Y,N) . : 

GENERATE  A  GRAPH?  (Y,N) . : 

LIST  OF  NODES  AND  COORDINATES?  (Y,N) . : 

LIST  ELEMENT  CONNECTIVITIES?  (Y,N) . : 

USE  SHELL  ELEMENTS?  ( Y , N) . : 

These  questions  serve  as  the  basic  controllers  of  the 
program  and  must  be  answered  ' Y'  or  'N'  for  yes  or  no. 

INTERACTIVE  RUN?  (Y ,N)  . : 

Requests  if  the  user  wishes  to  input  all  data  into  the 
program  dependent  on  prompting  questions  from  the  program 
(response  =  'Y').  Otherwise  the  program  looks  for  a  data  file 
on  TAPE4  set  up  as  described  later  in  this  report  under 
'WINGEN  -  MODEL  DEFINITION  DATA  FORMAT'  (response  =  'N'). 

CREATE  A  LOAD  DECK?  ( Y , N ) . : 

Should  the  user  response  ' Y'  to  this  question  an  input 
file  fully  formatted  with  CDC  job  control  cards  and  all  data 
necessary  for  a  MAGNA  FEM  (finite  element  model)  analysis  will 
be  created.  TAPEll  will  contain  this  information  at  the 
conclusion  of  the  run.  It  is  up  to  the  user  to  save  the  TAPEll 
file  as  a  permanent  file  for  later  use.  An  'N'  response 
prevents  a  load  deck  from  being  created.  Please  note:  if 
shell  elements  are  selected,  a  load  deck  will  be  created 
reg.  -dless  of  the  response  to  this  question  and  the  user  will 
need  to  utilize  PLOTBOB  postprocessing  program  to  obtain  model 
geometry  plots  of  the  structure  defined  on  TAPEll. 

GENERATE  A  GRAPH?  (Y,N) . : 

WINGEN  has  limited  plotting  capability  which  can  generate 
the  complete  structure  with  any  viewing  position  and  with  or 


without  nodes  and  elements  labeled.  The  user  must  have  access 
to  a  Tektronix  or  other  device  utilizing  PL0T1J8  graphics 
software  or  a  Hewlett-Packard  utilizing  HPPL0T21  library 
subroutines  to  make  use  of  the  graphics  capability  of  the 
program.  Further  information  is  detailed  in  Section  2.10 
’WINGEN  -  PLOTTING  CAPABILITIES'. 

USE  SHELL  ELEMENTS?  (Y,N) . : 

Shell  elements  (MAGNA  element  type  5)  are  an  optional 
substitution  for  the  membrane  plate  elements  (MAGNA  element 
type  3  subtype  1)  used  for  the  upper  and  lower  skins  of  the 
wing.  These  are  used  in  the  event  the  user  is  interested  in 
refining  the  element  mesh  to  be  finer  than  the  basic  geometry 
defined.  The  use  of  shell  elements  eliminates  the  need  for 
linear  constraints  to  be  generated  for  all  skin  nodes  not  bound 
by  geometric  restrictions  imposed  by  being  situated  at  the 
junction  of  rib,  spar,  and  skin  elements.  If  the  model  is 
refined  such  that  more  than  two  nodes  are  required  to  define 
the  thickness  through  any  part  of  the  structure  the  program 
will  generate  linear  constraints  to  prevent  unrealistic 
displacements  to  occur  during  analysis.  Shell  elements  should 
be  requested  for  damaged  models  as  well  as  refined  models  for 
more  accurate  analysis  results. 

Two  options  control  the  suppression  of  output  listing  of 
nodal  components  and  element  connectivities.  This  is  for  the 
case  where  the  terminal  operator  is  not  concerned  with  getting 
an  immediate  listing  of  the  node  and  element  data  and  may  elect 
instead  to  examine  the  FEM  input  data  stored  on  TAPEll  following 
the  WINGEN  run. 

LIST  NODES  AND  COORDINATES?  (Y,N) . : 

LIST  ELEMENT  CONNECTIVITIES?  (Y ,N) . . . . : 

An  'N'  response  to  either  of  these  questions  will  inhibit 
the  printing  of  the  appropriate  material  at  the  terminal. 


2.6  WINGEN  -  INTERACTIVE  MODE 

WINGEN  has  the  capability  for  fully  interactive  wing  model 
generation.  This  interactive  mode  requires  no  predefined  data 
sets  although  the  user  must  have  all  data  necessary  to  define 
the  model  available  for  the  terminal  session.  The  data  required 
is  defined  in  detail  in  Section  2.8  -  MODEL  DEFINITION  DATA 
FORMAT.  WINGEN  will  prompt  the  user  with  questions  such  as  those 
illustrated  in  Figure  2.6.1  where  the  underlined  values 
represent  responses  to  the  questions  prompted  by  the  program. 

All  information  requested  by  the  program  follows  the  same 
sequence  of  input  as  described  for  the  Model  Definition  data 
format.  The  questions  were  designed  to  guide  the  user  as  much  as 
possible  in  understanding  what  information  is  being  requested. 

A  glossary  is  provided  for  terms  which  are  utilized  in  this 
report  or  are  contained  in  questions  prompted  by  the  program  in 
the  event  the  questions  are  insufficient  for  the  user  to 
understand  what  is  being  requested.  Once  all  the  questions 
pertaining  to  the  model  definition  have  been  answered,  the 
program  will  build  a  WINGEN  input  data  file  from  the  parameters 
input  by  the  user.  If  a  load  deck  was  requested,  the  program 
will  then  initiate  further  questions  as  described  in  Section  2.9 
to  define  the  loaded  nodes,  fixed  nodes  and  analysis  parameters 
for  MAGNA.  Finally,  if  plots  were  requested,  the  program  will 
prompt  the  user  with  the  necessary  commands  to  generate  the 
model  just  defined.  The  plotting  capabilities  of  WINGEN  are 
fully  described  in  Section  2.10. 

Interactive  mode  is  initiated  as  described  in  Section  2.5. 
The  program  will  request  several  control  directives  to  be  input 
by  the  user  such  as: 

INTERACTIVE  RUN?  (Y,N) . : 

CREATE  A  LOAD  DECK?  (Y,N) . : 

GENERATE  A  GRAPH?  (Y , N) . : 

DO  YOU  WANT  A  LISTING  OF  NODAL  COORDINATES?  (Y,N) . : 

DO  YOU  WANT  A  LISTING  OF  ELEMENT  CONNECTIVITY?  (Y ,N) . . . : 

USE  SHELL  ELEMENTS?  ( Y , N ) . ' . : 


START 
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Figure  2.6.1. 


Sample  of  Interactive 


Mode  Questions. 
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The  user  must  respond  'Y'  to  the  question  -  Interactive  Run? 

This  will  direct  questions  to  the  user's  terminal  requesting 
input  of  necessary  data  for  defining  the  model.  The  remaining 
control  directives  may  be  answered  depending  on  the  user's 
wishes.  Once  the  program  issues  a  'START'  line,  the  interactive 
mode  has  been  engaged  and  questions  will  be  displayed  requesting 
input  data. 

The  input  data  required  by  the  preprocessor  falls  into 
twelve  categories: 

1.  Program  title 

2.  Profile  definition 

3.  Planform  description 

4.  Chordwise  description 

5.  Spanwise  description 

6.  Rib  definitions 

7.  Span  definitions 

8.  Modifications  to  basic  model 

9.  Post  definitions 

10.  Refinement  directives 

11.  Damage  specifications 

12.  Loads  to  be  applied 

These  12  categories  are  followed  sequentially  while  the  program 
requests  data  and  defines  the  model.  Each  of  these  categories 
is  fully  explained  and  all  values  are  defined  in  Section  2.8. 

The  user  should  refer  to  that  section  for  more  detailed 
information  concerning  the  data  requested  by  the  program.  Once 
the  data  has  been  fully  input,  the  program  will  catalog  the 
data  under  the  name ' USRDATAFILE ' .  This  file  may  be  attached 
as  a  preprocessor  input  data  file  (TAPE4)  and  used  to  generate 
the  same  model  in  later  executions  of  WINGEN  or  modified  and 
used  to  generate  a  new  model.  This  file  is  illustrated  in 
Figure  2.6.3.  The  user  should  be  aware  that  a  maximum  of  five 
permanent  files  with  the  same  name  (cycles)  are  allowed  on  the 
CDC .  If  repeated  use  of  the  interactive  mode  is  made,  the  user 
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will  have  to  ensure  there  is  a  cycle  available  on  which  to 
store  the  data  file.  Should  there  be  any  problem  (a  PF  abort, 
for  instance) ,  all  the  model  data  for  WINGEN  input  is  contained 
on  TAPE7  at  the  conclusion  of  the  program  execution.  The  user 
may  catalog  TAPE7  with  a  different  permanent  file  name  at  the 
conclusion  of  the  program  run.  Figure  2.6.2  illustrates  the 
T-38  wing.  The  appendices  provide  examples  of  attaching  and 
executing  WINGEN  to  make  use  of  the  various  processing  options. 
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WINGEN  SEMI-INTERACTIVE  MODE 


Capability  has  been  provided  for  a  "semi-interactive"  mode 
of  utilization  of  the  preprocessor  program  WINGEN.  This  mode 
differentiates  between  interactive  and  batch  modes  in  that  it 
allows  the  user  the  advantages  of  building  a  batch  input  data 
file  for  model  definition  and  permits  the  user  the  advantage  of 
graphic  verification  of  the  model  geometry  without  a  lengthy 
terminal  session.  The  semi-interactive  mode  is  utilized  by 
first  creating  a  model  data  file  as  defined  in  Section  2.8  and 
saving  the  file  or  attaching  a  cataloged  file  as  a  local  file  - 
TAPE4 .  The  user  then  follows  the  procedure  listed  in  Section 
2.5  to  initiate  WINGEN.  The  user  must  respond  'N'  to  the 
question  "Interactive  Run?"  All  required  model  definition  data 
will  be  taken  from  TAPE4 .  Once  the  model  has  been  defined, 
the  program  will  request  information  for  the  creation  of  a  load 
deck,  if  one  was  requested,  then  it  will  request  information 
for  plotting  parameters,  if  graphs  were  requested.  The  user 
may  refer  to  Sections  2.9  and  2.10,  respectively,  for  further 
information  on  load  decks  and  plotting  capabilities  of  WINGEN. 
Figure  2.7.1  illustrates  a  semi -interactive  mode  execution  of 
WINGEN  utilizing  the  model  definition  (batch  input)  data  format 
file  illustrated  in  Figure  2.7.2  to  generate  the  model  illustrated 
in  Figure  2.7.3. 


RIB  AMD  RIB-CAP  PARAMETERS  l 

RIB  NO  RATI  CAP  NATL  SECT  DIMENSIONS 


Figure  2.7.1.  (continued) 
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Figure  2.7.1.  (continued). 
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Figure  2.7.1.  (continued). 
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Figure  2.7.1.  (continued). 
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Figure  2.7.3.  Model  generated  with  the  input  batch  mode  data  format  illustrated 
in  Figure  2.7.2  during  the  semi-interactive  execution  of  the 
program  listed  in  Figure  2.7.1. 


2.8  WINGEN  -  MODEL  DEFINITION  DATA  FORMAT 

The  following  section  describes  the  Batch  mode  or 
Semi-interactive  mode  card  input  data  for  the  WINGEN  wing  model 
preprocessor  program.  The  definition  procedure  for  the  wing 
can  be  broken  down  into  12  sections  as  follows: 

1.  problem  title 

2.  wing  profile  definition 

3.  planform  description 

4.  wing  depth  description 

5.  wing  skin  properties  definition 

6.  ribs  definition 

7.  spars  definition 

8.  modify  directives 

9.  posts  definition 

10.  refinement  directives 

11.  damage  specifications 

12.  loads  definition 

For  each  item  of  input  a  corresponding  FORTRAN  variable 
name  is  listed.  Unless  otherwise  noted,  the  type  of  the  input 
variable  corresponds  to  the  standard  FORTRAN  naming  conventions 
(names  beginning  with  the  letters  I  -  N  are  integer;  all  others 
are  floating  point  numbers) .  All  floating  point  data  may  be 
input  with  or  without  exponents  in  the  data  field  provided. 
Integers  and  exponents  must  be  right- justified  in  the  data 
field.  The  terms  utilized  extensively  in  this  section  are 
illustrated  in  Figure  2.8.1  and  a  glossary  is  provided  to  assist 
the  user  with  term  definition. 

The  user  must  include  one  card  for  each  of  the  twelve 
categories  as  explained  below  in  each  section.  When  the  file  is 
completed  the  user  should  save  it  as  a  permanent  disk  file 
and  attach  it  as  TAPE 4  for  use  with  Batch  or  Semi -interactive 
WINGEN  execution  modes.  This  section  also  describes  the  data 
required  for  full  interactive  execution  of  WINGEN,  In  this 
case  the  program  will  prompt  the  user  for  the  appropriate  data. 
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2.8.1  Problem  Title 


(required) 


CARD 

COLUMN 

DATA 

DESCRIPTION 

NOTES 

1 

1-80 

TITLE 

Alphanumeric  problem  title 

1 

NOTES : 

(1)  This  title  is  utilized  in  the  postprocessor 
program  CONTOUR  for  labeling  the  plot  and  as 
a  header  on  all  pages  of  the  MAGNA  analysis 
run  and  WINGEN  output  listing. 


s 


2.8.2  Wing  Profile  Parameters 
(required) 


CARD 

COLUMN 

DATA 

DESCRIPTION 

NOTES 

1 

1-7 

PROFILE 

Literal  'PROFILE1 

8 

(blank) 

- 

9-10 

NP 

Planform  shape  code  literal: 

1 

'Pi'  =  rectangular 
'P2*  =  full  swept 
'P3*  =  piecewise  linear 
'P4*  =  general 

11 

(blank ) 

12-13 

NC 

Chordwise  depth  distribution 
code  literal: 

2 

'Cl'  =  constant  chord  depth 
'C2*  =  linear  chord  depth 
'C3'  =  piecewise  linear 

chord  depth 

'C4'  =  general 

14 

(blank) 

15-16 

NS 

Spanwise  depth  distribution 
code  literal: 

3 

' SI '  =  constant  depth  span 
'S2'  =  linear  depth  span 
'S3'  =  piecewise  linear 

depth  span 

'S4'  =  general 

NOTES : 

(1)  Diagramed  in  Figure  2.8.2  are  the  three  program- 
defined  planform  shapes  corresponding  to  Pi,  P2  and 
P3.  P4  is  a  general  class  and  should  only  be  used 
with  chord  and  spanwise  depth  distributions  C4  and  S4 

(2)  The  three  program-defined  chordwise  depth 
distributions  for  Cl  -  C3  codes  are  illustrated  in 
Figure  2.8.3.  Chordwise  depth  code  C4  is  a  general 
class  and  should  only  be  used  with  planform  and 
spanwise  codes  P4  and  S4 . 
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(3)  The  three  program-defined  spanwise  depth  distributions 
for  SI  -  S3  codes  are  illustrated  in  Figure  2.8.4. 
Spanwise  depth  code  S4  should  only  be  used  in 
conjunction  with  planform  and  chordwise  depth  codes 
P4  and  C4 . 


PROFILE  CONFIGURATION 


PLANFORM 


DESCRIPTION 


RECTANGULAR 


FULL  SWEPT 


PIECEWISE  SWEPT 


Figure  2.8.2.  Three  Planform  Classes. 


SKETCH 


PROFILE 

CONFIGURATION 

CODE 

DESCRIPTION 

CHORDWISE 

C  1 

CONSTANT  DEPTH 

C  2 

LINEAR  DEPTH 

C  3 

PIECEWISE  LINEAR  DEPTH 

SKETCH 


Figure  2.8.3.  Three  Chordwise  Depth  Classes 
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CODE 

DESCRIPTION 

SPANWISE 

S  1 

CONSTANT  DEPTH 

S  2 

LINEAR  DEPTH 

S3 

PIECEWISE  LINEAR  DEPTH 

SKETCH 


Figure  2.8.4.  Three  Spanwise  Depth  Classes 


2.8.3  Planform  Description 
(required) 


A.  Profile  Type  PI 


CARD 

- - 

COLUMN 

DATA 

DESCRIPTION 

NOTES 

1 

1-4 

PLAN 

Literal  'PLAN' 

_ 

5 

(blank) 

- 

6-15 

ROOTC 

Root  chord  dimension 

16-25 

SPAND 

Span  dimension 

2 

NOTES : 

(1)  The  root  chord  dimension  is  measured  from  fore  to 
aft  (leading  edge  to  trailing  edge  of  wing)  along 
the  root  edge  of  the  wing,  as  illustrated  in 
Figure  2.8.5. 

(2)  The  span  dimension  is  the  distance  from  the  wing 
root  to  the  wing  tip  measured  perpendicular  to 
the  root  chord.  This  is  illustrated  in 

Figure  2.8.5. 


I 


Figure  2.8.5.  Root  chord  dimension  and  span  dimension. 


B.  Profile  Type  P2 


CARD 

COLUMN 

DATA 

DESCRIPTION 

NOTES 

1 

1-4 

PLAN 

Literal  ’PLAN* 

_ 

5 

(blank) 

6-15 

ROOTC 

Root  chord  dimension 

1 

16-25 

SPAND 

Span  dimension 

2 

26-35 

TIPC 

Tip  chord  dimension 

3 

36-45 

SWEEP 

Leading  edge  sweep  angle 
(degrees ) 

4 

NOTES : 

(1)  The  root  chord  dimension  is  measured  from  fore  to 
aft  (leading  edge  to  trailing  edge)  along  the  root 
edge  of  the  wing,  as  illustrated  in  Figure  2.8.5. 

(2)  The  span  dimension  is  the  distance  from  the  wing 
root  to  the  wing  tip  measured  perpendicular  to 

the  root  chord.  This  is  illustrated  in  Figure  2.8.5. 

(3)  The  tip  chord  dimension  is  measured  from  fore  to 
aft  (leading  edge  to  trailing  edge)  along  the  tip 
edge  of  the  wing.  This  is  illustrated  in  Figure  2.8.5 

(4)  The  leading  edge  sweep  angle  is  the  angle  in  degrees 
that  the  leading  edge  of  the  wing  deviates  from  the 
perpendicular  intersection  of  the  root  chord.  This 
is  illustrated  in  Figure  2.8.6. 
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C.  Profile  type  P3 


COLUMN 

DATA 

DESCRIPTION 

1-4 

PLAN 

Literal  ’PLAN' 

5 

(blank) 

6-10 

NSEGS 

Number  of  spanwise 
stations  to  be  used  in 
defining  planform  shape 

1-10 

XLE 

Leading  edge  x-coordinate 

11-20 

YLE 

Leading  edge  y-coordinate 

21-30 

XTE 

1  Trailing  edge  x-coordinate 

31-40 

YTE 

Trailing  edge  y-coordinate 

NOTES 


NOTES : 


A  spanwise  station  is  defined  for  this  program  as 
being  the  point  where  two  wing  sections  with 
differing  sweep  angles  intersect.  For  N  spanwise 
sections  there  must  be  N+l  spanwise  stations 
defined  as  illustrated  in  Figure  2.8.7. 

Each  spanwise  section  is  to  be  defined  beginning 
at  the  leading  edge  root  segment  and  moving 
outboard.  Repeat  card  2  NSEGS  times  to  define  all 
spanwise  segments.  The  origin  of  coordinates  is 
the  leading  edge  root  point;  x  is  measured  positive 
aft  and  y  is  measured  positive  outboard. 
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(outboard) 
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Figure  2.8.7.  A  new  spanwise  depth  station  must  be  defined 

for  each  change  in  wing  thickness.  This  model 
has  four  depth  stations  defined  (NSD.-NSD. ).  For 
comparison  the  seven  spanwise  nodal  stations 
are  also  labeled  (NSN.-NSN-),  one  for  each  rib 
in  the  model. 
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D.  Profile  type  P4 


This  section  defines  the  Direct  Nodal  Point  Input 
for  the  general  class  of  wings  designated  as  P4-C4-S4  in  the 
planform  definition. 


16-20 

21-25 

26-30 

31-35 

36-40 

41-45 


6-10 

11-15 

# 

76-80 


DATA 

DESCRIPTION 

NSD 

Number  of  spanwise  nodal 
stations  to  be  input 

NODES 

Total  number  of  nodal 
points 

ISW 

Number  of  pairs  of  nodal 
coordinates  to  be  switches 

ISW1 (1) 

Index  of  1st  coordinate  in 

1st  pair  switched 

ISW2 (1) 

Index  of  2nd  coordinate  in 

1st  pair  switched 

ISW1 (2) 

Index  of  1st  coordinate  in 
2nd  pair  switched 

ISW2 (2) 

Index  of  2nd  coordinate  in 
2nd  pair  switched 

ISWl (3) 

Index  of  1st  coordinate  in 

3rd  pair  switched 

ISW2 (3) 

Index  of  2nd  coordinate  in 

3rd  pair  switched 

FACT (1 ) 

Scale  factor  for  (switched) 
x-coordinate 

FACT (2) 

Scale  factor  for  (switched) 
y-coordinate 

FACT { 3 ) 

Scale  factor  for  (switched) 
z-coordinate 

FORMAT 

FORMAT  (including 
parenthesis)  for  nodal 
coordinate  input  to  follow 

NL  (1 ) 

Spanwise  nodal  stations 
input 

NL  (2 ) 

NL  { 3 ) 

• 

• 

NL (16 ) 

up  to  16  per  card 

NOTES 


CARD 


COLUMN 


DATA 


DESCRIPTION 


NOTES 


5 


NODE 

X 

Y 

Z 


Node  number 
x-coordinate  of  node 
y-coordinate  of  node 
z -coordinate  of  node 


8 

8 

8 

8 


NOTES : 

(1)  NSD  must  be  greater  than  0.  This  parameter  defines 
the  number  of  spanwise  nodal  stations  to  be  input 
on  card  4.  A  spanwise  nodal  station  required  for 
this  input  is  the  last  node  defined  for  each  row 

of  nodes  located  between  the  leading  and  trailing 
edges  (spanwise) .  This  is  illustrated  in  Figure 
2.8.8  where  A,  B,  C,  D,  E  and  F  denote  spanwise 
nodal  stations. 

(2)  The  P4-C4-S4  planform  type  requires  the  user  to 
input  all  nodes  to  WINGEN  in  x-,  y-  and 
z-coordinate  form.  NODES  is  the  total  number  of 
nodes  required  to  define  the  model.  The  program 
will  request  format  specifications  on  card  3  for 
input  of  NODES  nodal  points  from  card  5. 

(3)  A  special  option  for  translating  nodal  coordinates 
from  one  axis  to  another  (e.g.,  from  x  to  y  or  z  to 
x  and  vice  versa)  is  provided.  This  is  for  the 
case  where  a  previously  defined  model  perhaps 
developed  for  an  alternative  analysis  program  is 
being  input  for  analysis  by  MAGNA.  The  switched 
nodal  coordinates  allow  for  the  translation  of  the 
model  in  space  so  it  will  conform  to  the  x,  y,  z 
coordinate  axes  established  by  this  program  where 
the  wing  model  coordinate  origin  is  at  the  junction 
of  the  leading  edge,  the  root  chord  and  centered  at 


Il _ _  1 


.  1  _  X.  .  J 


•T*  r- 


the  midpoint  of  the  wing  depth  with  the  z-axis 
vertical,  the  x-axis  extending  chordwise,  fore  to 
aft,  and  the  y-axis  extending  spanwise  root  to  tip 
(outboard) . 

TSW  is  the  number  of  coordinate  groups  to  be 
switched  where 

0  =  no  coordinate  transformations 
1=1  pair  of  coordinate  transformations 
2  =  2  pair  of  coordinate  transformations 
3=3  pair  of  coordinate  transformations. 

The  maximum  number  of  transformations  is  3. 

(4)  Coordinate  transformations  must  be  entered  in  pairs 
where  1  =  x,  2  =  y  and  3  =  z.  The  transformation 
should  result  in  having  x  at  1  y  at  2  and  z  at  3. 

For  example,  if  data  being  input  is  as  follows:  node, 
z,  x,  y.  Then  two  pairs  of  coordinate  transformations 
are  necessary: 

1st  pair  of  switches  is  1  and  2  which  leaves 
x  in  the  1st  array  position  and  z  in  the  2nd  array 
position.  2nd  pair  of  switches  is  2  and  3  which 
leaves  y  in  the  2nd  array  position  and  z  in  the 
3rd  array  position. 

The  value  range  for  ISWl  and  ISW2  is  from  1  to  3. 

(5)  FACT ( 1 )  -  FACTO)  correspond  to  coordinate 
multipliers  for  altering  the  model  origin  of 

the  nodes  being  input.  All  nodes  will  be  multiplied 
by  these  factors  therefore  a  factor  of  1.0  must  be 
specified  for  all  three  factors  if  no  coordinate 
shifting  is  to  take  place.  FACT(l),  FACT(2)  and 
FACT (3)  represent  the  (transformed)  x-,  y-  and 
z-coordinate  multipliers,  respectively. 
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(6)  FORMAT  must  contain  the  entire  format  including 
parenthesis  according  to  CDC  FORMAT  statement 
syntax.  All  data  on  card  5  will  be  input  utilizing 
this  format.  Each  read  using  FORMAT  will  look  for 
NODE,  X,  Y,  Z  values. 

(7)  Repeat  card  4  enough  times  to  define  NSD  spanwise 
nodal  stations  (e.g.,  NSD/16  =  #  card  4  required). 
Each  card  4  may  have  up  to  16  nodal  stations.  See 
notes  1  for  more  information  on  spanwise  nodal 
stations. 

(8)  Repeat  card  5  until  NODES  nodes  have  been  input. 

The  format  to  input  these  values  was  defined  on 
card  3.  NODE  should  begin  at  1  and  increment 
sequentially  by  1.  The  x-,  y-  and  z-coordinates 
will  all  be  multiplied  by  the  appropriate  FACT(l), 
FACT (2)  or  FACT (3)  as  defined  on  card  2. 
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2.8.4  Wing  Depth  Distribution 

(omit  for  P4-C4-S4  profile  type) 

A.  Constant  Depth  -  for  span  type  Si;  chord  type 
Cl  only.  (Refer  to  Figure  2.8.4.) 


CARD 

COLUMN 

DATA 

DESCRIPTION 

NOTES 

1 

1-5 

1-10 

DEPTH 

WDEPTH 

Literal  'DEPTH' 

Total  wing  depth 

1 

NOTES; 


(1)  This  is  for  a  wing  with  constant  wing  depth 
(thickness)  throughout. 
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B.  Chordwise  Variation  Only  -  for  span  type  SI;  chord 
types  C2  and  C3  only.  (Refer  to  Figure  2.8.4.) 


CARD 

COLUMN 

DATA 

DESCRIPTION 

NOTES 

1 

1-5 

DEPTH 

Literal  'DEPTH' 

_ 

6-10 

(blank) 

- 

11-15 

NCD 

Number  of  chordwise  stations 

at  which  depths  will  be 
specified 

1 

2 

1-10 

X 

x-coordinate  for  depth 
stations 

2 

11-20 

WDEPTH 

Wing  depth  at  this  station 

3 

NOTES : 

(1)  Repeat  card  2  NCD  times  to  define  all  chordwise 
depth  stations.  Begin  at  the  leading  edge  and 
proceed  to  the  trailing  edge  along  the  root  chord 
to  define  NCD  depth  stations.  There  must  be  a 
minimum  of  2  such  stations  specified.  See 
figure  2.8.9. 

(2)  The  x-coordinate  requested  is  in  the  units  system 
supplied  by  the  user  in  defining  the  model. 

(3)  The  wing  depth  should  be  measured  from  the  top  of 
the  lower  skin  to  the  bottom  of  the  upper  skin 
for  most  accurate  results. 
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Figure  2.8.9.  A  new  chordwise  depth  station  must  be  defined 

for  every  change  in  wing  thickness.  This  model 
has  four  depth  stations  defined  (NCD 
for  comparison  five  chordwise  nodal 
(NCN.-NCN-) ;  one  for  each  spar  defined  in  the 
modet. 


.  -NCDg )  and 
stations 
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C.  Spanwise  Depth  Variation  -  for  span  types  S2  and  S3j 

chord  types  Cl,  C2  and  C3  only.  (Refer  to  Figure  2.8.4.) 


CARD 

COLUMN 

DATA 

DESCRIPTION 

NOTES 

1 

1-5 

DEPTH 

Literal  'DEPTH* 

6-10 

NSD 

Number  of  spanwise  depth 
stations  for  which  data 
will  be  specified. 

1 

11-15 

NCD 

Number  of  chordwise  depth 
stations  for  which  data 
will  be  specified 

1 

2 

1-10 

X 

Depth  station  x  coordinate 

2 

11-20 

Y 

Depth  station  y  coordinate 

2 

21-30 

WDEPTH 

Wing  depth  at  this  coordinate 

3 

NOTES : 

(1)  Repeat  card  2  for  each  NCD  and  NSD  station  to  be 
defined;  a  total  of  NSD*NCD  cards  will  be  required. 
A  spanwise  depth  station  is  that  point  where 

the  depth  of  the  wing  changes  as  one  travels  from 
root  to  tip  along  the  wing.  A  chordwise  depth 
station  is  that  point  where  the  depth  of  the  wing 
changes  in  moving  from  leading  edge  to 
trailing  edge. 

(2)  The  user  should  define  the  depths  starting  at  the 
wing  coordinate  origin  and  moving  along  the  root 
chord  (fore  to  aft)  then  spanwise  to  the  next 
depth  station  along  the  leading  edge  and  working 
to  the  trailing  edge,  etc.,  until  all  depth 
stations  have  been  defined. 


(3)  Wing  depth  is  measured  from  the  top  of  the  lower 
skin  to  the  bottom  of  the  upper  skin. 


2.8.5  Wing  Skin  Properties 
(required) 


COLUMN 

DATA 

DESCRIPTION 

1-4 

SKIN 

Literal  'SKIN' 

5 

(blank) 

6-10 

NMATU 

Upper  skin  material 
property  code 

11-15 

. 

NMATL 

Lower  skin  material 
property  code 

16-25 

TUPPER 

Upper  skin  thickness 

26-35 

TLOWER 

Lower  skin  thickness 

NOTES 


Material 


Elastic  Modulus  lb/in' 


Poisson's  Ratio 


EQ  Stress  at  1st  Yield 
lb/in2 


LINEAR  AND  NONLINEAR  ANALYSIS 

_ _ _ — _ _ i 

Material  1 

Material  2 

aluminum 

steel 

.10*10® 

.30*10® 

.  30 

.33 

. 259*10~3 

.725*10~3 

.10*105 

. 30*105 

NONLINEAR  i 

ANALYSIS  ONLY 

Strain  Hardening  Curve 


Uniaxial  Stress/Strain 
Data  Curve 
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The  user  may  select  either  one  or  add  additional 
ones  as  required. 

(2)  Skin  thicknesses  should  approximate  as  closely 

as  possible  the  true  skin  thicknesses  of  the  wing. 
If  considerable  variation  exists  alterations  may 
have  to  be  made  to  the  WINGEN  output  file  (TAPEll) 
prior  to  an  analysis  run.  Note  that  in  the 
generation  of  shell  elements,  the  upper  and 
lower  skin  thicknesses  are  utilized  to  generate 
new  nodes  in  the  z -direction  at  TUPPER  and  TLOWER 
displacements  from  the  original  nodes  defining 
the  membrane  elements.  The  method  of  shell 
element  generation  is  illustrated  in  Figure  2.8.10. 
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skin  thickness  =  TUPPER 


B 


16  20  24 


skin  thickness  =  TLOWER 


Two  chord  bays  of  a  span  bay  are  used  here  to 
illustrate  the  thickness  expansion  for  the 
upper  and  lower  skin  elements  in  converting  from 
plate  elements  to  shell  elements.  A)  Nodes  of 
basic  model  before  shell  element  conversion  B) 
Nodes  of  same  elements  after  conversion  to  shell 
elements . 


2.8.6  Rib  Definition 
(required) 


CARD 

COLUMN 

DATA 

DESCRIPTION 

NOTES 

1 

1-4 

RIBS 

Literal  'RIBS' 

5 

(blank) 

6-10 

NRIB 

Number  of  ribs  to  be 
defined 

1 

2 

1-5 

MAT 

Material  property  code  for 
rib 

2 

6-10 

MATC 

Material  property  code  for 
rib  cap 

2 

16-25 

XI 

Leading  edge  x-coordinate 
of  rib 

3 

26-35 

YI 

Leading  edge  y-coordinate 
of  rib 

3 

36-45 

XO 

Trailing  edge  x-coordinate 
of  rib 

3 

46-55 

YO 

Trailing  edge  y-coordinate 
of  rib 

3 

3 

1-10 

RPROP(l) 

Rib  thickness 

11-20 

RPROP (2) 

Rib  cap  cross-sectional  area 

4 

NOTES : 

(1)  Repeat  cards  2  and  3  to  define  NRIB  ribs.  If 
NRIB  =  0  or  blank,  card  1  completes  this  item. 

(2)  WINGEN  predefines  two  material  property  codes 
as  follows: 


LINEAR  AND  NONLINEAR  ANALYSIS 

Material  1 

Material  2 

Material 

aluminum 

steel 

2 

Elastic  Modulus  lb/in 

.10*108 

.30*10® 

Poisson's  Ratio 

.30 

.33 
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Mass  Density 

..  2,.  4 

lb  sec  /xn 

.259*10~3 

. 725*10~3 

EQ  Stress  at  1st  Yield 
lb/in^ 

.10*105 

. 30*105 

NONLINEAR  ANALYSIS  ONLY 

Strain  Hardening  Curve 

1 

1 

Uniaxial  Stress/Strain 
Data  Curve 

2 

2 

Additional  codes  may  be  added  but  the  user  must 
be  careful  to  define  these  material  property  codes 
in  the  WINGEN  output  file  (TAPEll)  prior  to 
executing  an  analysis.  If  MAT  =  0  or  blank  no 
rib  web  is  defined  for  the  rib  definition.  If 
MATC  =  0  or  blank  no  rib  caps  are  defined  for  the 
rib  definition. 

(3)  Define  each  rib  in  terms  of  leading  and  trailing 
edge  x- ,  and  y-components;  z-components  will  be 
supplied  based  on  wing  depth  distributions. 

To  define  each  rib  begin  with  the  rib  closest  to 
the  wing  coordinate  origin  and  wing  root  and  work 
from  the  wing  root  to  the  wing  tip  defining  NRIB 
ribs  as  diagrammed  in  Figure  2.8.11. 

(4)  Each  rib  web  will  have  associated  with  it  an  upper 
and  lower  rib  cap  if  MATC  >  0.  The  user  must  define 
the  cross-sectional  area  of  these  rib  caps  (note: 
this  is  not  a  thickness  parameter  but  a  cross- 
section  area  value).  This  is  illustrated  in 
Figure  2.3.1. 
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2.8.7  Spar  Definition 


(required) 


CARD 

COLUMN 

DATA 

DESCRIPTION 

NOTES 

1 

1-5 

SPARS 

Literal  'SPARS' 

6-10 

NSPAR 

Number  of  spars  to  be 
defined 

1 

2 

1-5 

MAT 

Material  property  code  for 
spar 

2 

6-10 

MATC 

Material  property  code  for 
spar  cap 

2 

16-25 

XI 

Inboard  x-coordinate  of 
si  r 

3 

26-35 

YI 

Inboard  y-coordinate  of 
spar 

3 

36-45 

XO 

Outboard  x-coordinate  of 
spar 

3 

46-55 

YO 

Outboard  y-coordinate  of 
spar 

3 

3 

1-10 

SPROP (1) 

Spar  thickness 

- 

11-20 

SPROP (2 ) 

Spar  cap  cross-sectional 

area 

4 

NOTES : 

(1)  Repeat  cards  2  and  3  until  NSPAR  spars  have  been 
defined.  If  NSPAR  =  0  or  blank  card  1  completes 
this  item. 


(2)  The  program  predefines  two  material  property 
codes  as  follows: 


LINEAR  AND  NONLINEAR  ANALYSIS 

Material  1 

Material  2 

Material 

aluminum 

steel 

Elastic  Modulus 
lb/in2 

.10*108 

.30*10® 

Poisson’s  Ratio 


.30 


.33 


2  /  •  4 

Mass  Density  lb  sec  /in 

. 259*10-3 

.725*10"3 

EQ  Stress  at  1st  Yield 
Ih/i-Vi 

.io*io5 

. 30*105 

NONLINEAR  ANALYSIS  ONLY 

Strain  Hardening  Curve 

1 

1 

Uniaxial  Stress/Strain 

Data  Curve 

2 

2 

The  user  may  use  these  or  may  require  additional 
ones.  If  codes  other  than  1  or  2  are  used  the 
user  must  take  care  to  add  these  material  properties 
to  the  WINGEN  output  file  (TAPEll)  prior  to 
executing  an  analysis  run. 

If  MAT  =  0  or  blank  no  spar  webs  will  be  defined 
for  this  spar.  If  MATC  =  0  or  blank  no  spar  caps 
will  be  defined  for  this  spar. 

(3)  Define  each  spar  in  terms  of  inboard  and  outboard 
x-  and  y-coordinates;  the  z-component  is  supplied 
by  the  program  based  on  wing  depth  distribution. 

This  is  illustrated  in  Figure  2.8.12. 

(4)  Each  spar  web  will  have  an  associated  spar  cap  if 
MATC  >  0.  The  user  must  define  the  cross-sectional 
area  of  these  caps.  (Note:  this  is  not  a  thickness 
parameter  but  a  cross-section  area  value) .  This  is 
illustrated  in  Figure  2.3.1. 
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Figure  2.8.12. 


WINGEN  will  define  nine  spars  for  this  model 
in  the  sequence  indicated. 
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CARD 

COLUMN 

DATA 

DESCRIPTION 

NOTES 

1 

1-6 

MODIFY 

Literal  'MODIFY' 

_ 

7-10 

ND 

Number  of  modify  data 
specifications 

1 

A.  Class  1  -  Removal  of  Selected  Elements  j 

j 

'j 


CARD 

COLUMN 

DATA 

— 

DESCRIPTION 

NOTES 

2 

1-7 

ID 

Modify  type  literal 

2 

8-10 

(blank) 

11-15 

INDXI 

Spanwise  bay  of  element  to 
be  removed 

3 

16-20 

INDXJ 

Chordwise  bay  of  element  to 
be  removed 

3 

B.  Class  2  -  Removal  of  Elements  in  a  Specified  Region 


CARD 

COLUMN 

DATA 

DESCRIPTION 

NOTES 

2 

1-7 

ID 

Literal  'LOCATE'  (modify  type) 

__ 

8-20 

(blank) 

21-30 

DX 

x-coordinate  of  centroid  of 

area  to  be  removed 

4 

31-40 

DY 

y-coordinate  of  centroid  of 

area  to  be  removed 

4 

41-50 

RAD 

Spherical  radius  of  area  to  be 

removed 

4 

NOTES : 


(1)  Repeat  card  2  ND  times  to  define  all  modifications 
to  model.  Classes  1  and  2  may  be  mixed  freely. 

If  ND  =  0  or  blank  then  card  1  completes  this  input 
item. 


(2)  There  are  eight  designations  for  Class  1  modifications. 
These  literals  should  be  entered  left  justified.  All 
Class  1  and  Class  2  modifications  should  equal  ND 
modifications.  The  designations  for  element  removal 
are: 


ID  =  SKINL 
SKINU 
RIB 
RCAPL 
RCAPU 
SPAR 
SCAPL 
SCAPU 


for  lower  skin  element 
for  upper  skin  element 
for  rib  panel 
for  lower  rib  cap 
for  upper  rib  cap 
for  spar  panel 
for  lower  spar  cap 
for  upper  spar  cap 


These  may  be  entered  in  any  combination  or  order 
to  alter  the  model  elements. 


(3)  Spanwise  bays  are  numbered  from  root  to  tip;  chordwise 
bays  are  numbered  leading  to  trailing  edges.  For 
spars  and  spar  caps,  INDXI  is  the  spanwise  nodal 

line  affected;  for  ribs  and  rib  caps,  INDXJ  denotes 
the  chordwise  nodal  line  affected. 

(4)  DX  and  DY  specify  the  center  of  a  sphere  of  radius 
RAD.  All  elements  whose  centroids  lie  within  this 
sphere  will  be  removed.  All  values  are  given  in 
model  coordinates  supplied  by  the  user.  It  is 
assumed  the  user  will  generate  an  undamaged  model 
first,  then  will  generate  a  damaged  model.  DX,  DY 
and  RAD  should  be  easily  obtainable  from  the  nodal 
coordinates  data  of  the  model. 
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2.8.9  Post  Definition 


(required) 


CARD 

COLUMN 

DATA 

DESCRIPTION 

NOTES 

1 

1-5 

POSTS 

-  ' 

Literal  'POSTS' 

6-10 

NPST 

Number  of  posts  to  be 
defined  ; 

1 

2 

1-5 

IC 

Chordwise  station 

2 

6-10 

IS 

Spanwise  station 

2 

11-15 

MAT 

Material  property  code  for 
post 

3 

3 

1-10 

BPROP (1) 

(blank) 

11-20 

BPROP (2 ) 

Post  cross-sectional  area 

4 

NOTES : 

(1)  If  NPST  >  0,  repeat  cards  2  and  3  for  each  post  to 
be  defined  in  the  model.  If  NPST  <  0,  insert  only 
one  of  each  card  2  and  3  leaving  columns  1-10  blank 
on  card  2 .  The  program  will  generate  a  post  at 
each  rib-spar  intersection  with  the  given  material 
properties.  A  note  of  caution:  defining  a  model 
without  posts  will  yield  erroneous  analysis  results 
due  to  collapsing  of  the  model  under  loading. 

(2)  The  chordwise  stations  are  defined  from  the  wing 
coordinate  origin  located  at  the  root  chord 
(inboard)  and  leading  edge  (fore) .  To  determine 
the  appropriate  chordwise  station  move  aft  (to 
trailing  edge)  along  the  root  chord  one  station  per 
defined  chordwise  depth  station  as  defined  in  2.8.4. 
Wing  Depth  Distribution.  This  is  illustrated  in 
Figure  2.8.13.  It  is  best  to  define  a  post  at  the 
point  of  intersection  of  rib  and  spar  and  have  it 
extend  from  the  lower  skin  to  the  upper  skin.  The 


spanwise  stations  are  defined  beginning  at  the 
wing  coordinate  origin  and  moving  outboard  (to  the 
wing  tip) .  Each  spanwise  depth  station  represents 
a  spanwise  station  as  illustrated  below  and  defined 
in  2.8.4  Wing  Depth  Distribution.  The  parameters 
IC  and  IS  effectively  become  a  local  coordinate 
system  for  defining  where  a  post  should  be  placed 
where  IC  is  in  the  chordwise  direction  along  the 
root  chord  and  IS  is  in  the  spanwise  direction 
along  the  leading  edge. 

(3)  The  program  predefines  two  material  property  codes 
to  be  used  for  model  definition: 


LINEAR  AND  NONLINEAR  ANALYSIS 


Material  1 


Material  2 


Material 


aluminum 


steel 


.  2 


Elastic  Modulus  lb/in 


.10*10 


8 


.30*10 


8 


2  4 

Mass  Density  lb  sec  /in 


.259*10 


-3 


.725*10 


-3 


EQ  Stress  at  1st  Yield 
lb/in2 


.10*10' 


30*10' 


NONLINEAR  ANALYSIS  ONLY 


Strain  Hardening  Curve 


Uniaxial  Stress/Strain 
Data  Curve 


Additional  codes  may  be  entered  for  preprocessor 
generation  but  the  user  must  be  careful  to  define 
any  additional  codes  he  utilizes  by  adding  the 
appropriate  data  to  the  WINGEN  output  file 
(TAPEll ) . 
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(4)  The  user  must  define  a  cross-sectional  area  for 
the  posts.  This  area  should  be  neither  too 
small  nor  too  large  or  abnormally  large  stresses 
and  strains  will  develop.  (An  area  of  .001  -  .1 
should  be  considered) . 


Illustration  of  chordwise  and  spanwise  nodal  stations 
designation  to  place  a  post  at  location  P  would  be 


2.8.10  Refinement  Directives 


(required) 

Any  refinements  made  to  the  skins  will  result  in 
meaningless  analysis  results  unless  the  shell  elements  option 
is  requested  while  executing  WINGEN.  See  Section  2.3  WINGEN  - 
MODEL  DEFINITION  for  more  information  about  this. 


CARD 

COLUMN 

DATA 

DESCRIPTION 

NOTES 

1 

REFINE 

Literal  'REFINE' 

__ 

7-10 

NREF 

Number  of  refinement  data 
specifications 

1 

2 

1-7 

TYPE 

Refinement  type  literal: 

2 

'SPAN'  =  spanwise 
refinement 
'CHORD'  =  chordwise 
refinement 

'THICK'  =  thickness  (depth) 

refinement 

8-10 

(blank) 

■s 

11-15 

IDIV 

Number  of  divisions 

16-20 

_ 

IBAY 

Bay  number 

H 

NOTES : 

(1)  Repeat  card  2  NREF  times  to  define  all  refinement 
specifications.  If  NREF  =  0  or  is  blank,  no 
refinement  is  to  be  specified  and  card  1  completes 
this  input  item. 

(2)  Three  refinement  types  are  available: 

'SPAN'  results  in  a  spanwise  refinement 
'CHORD'  results  in  a  chordwise  refinement 
'THICK'  results  in  a  depthwise  refinement. 

The  repeated  use  of  any  one  type  in  a  single  execution 
of  WINGEN  will  yield  additive  results  unless  the 
IBAY  option  is  utilized.  Any  number  of  each  or 
all  types  may  be  input  in  any  order.  Figure  2.8.14 
illustrates  the  three  refinement  types. 
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(3)  IDIV  represents  the  number  of  divisions  each 
specified  bay  will  be  divided  by.  If  IDIV  =  1 
each  bay  affected  will  be  divided  in  half  (one 
additional  bay) .  If  IDIV  =  2  each  bay  affected 
will  be  divided  into  thirds  (two  additional  bays) . 

(4)  For  I BAY  =  0  or  blank,  all  bays  in  the  refinement 
type  will  be  increased  by  IDIV  divisions.  If 
IBAY  >  0  only  the  IBAY  in  the  refinement  type  will 
be  increased  by  the  number  of  divisions. 
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Ik 


'-3  8  wing  model  with  THICK  refinement  of  1 
livision  and  nodes  labeled. 
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Figure  2.8.14f.  T-38  wing  model  root  chord  elements  before 

and  after  THICK  refinement  with  nodes 
labeled. 


2.8.11  Damage  Specifications 


(required) 


CARD 

COLUMN 

DATA 

DESCRIPTION 

NOTES 

1 

1-6 

DAMAGE 

Literal  'DAMAGE* 

7-10 

ND 

Number  of  damage  data 
specifications 

1 

A.  Class  1  -  Deletion  of  Selected  Elements 


CARD 

COLUMN 

DATA 

DESCRIPTION 

NOTES 

2 

1-7 

ID 

Damage  type  literal 

2 

8-10 

(blank) 

- 

11-15 

INDXI 

Spanwise  bay  of  element  to 
be  deleted 

3 

16-20 

INDXJ 

Chordwise  bay  of  element  to 
be  deleted 

3 

B.  Class  2  -  Deletion  of  Elements  in  a  Specified  Region 


CARD 

COLUMN 

DATA 

DESCRIPTION 

NOTES 

2 

1-7 

ID 

Literal  'LOCATE'  (damage 

type) 

- 

8-20 

(blank) 

21-30 

DX 

x-coordinate  of  centroid  of 

area  to  be  deleted 

4 

31-40 

DY 

y-coordinate  of  centroid  of 

area  to  be  deleted 

4 

41-50 

RAD 

Spherical  radius  of  area  to 

be  deleted 

4 

I 


i 

i 
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NOTES : 


(1)  Repeat  card  2  ND  times  to  define  all  damages  to 
model.  Classes  1  and  2  may  be  mixed  freely.  If 

ND  =  0  or  blank,  then  card  1  completes  this  input  item 

(2)  There  are  eight  designations  for  Class  1  damages. 

These  literals  should  be  entered  left  justified. 

All  Class  1  plus  Class  2  deletions  should  equal  ND 
damage  specifications.  The  designations  for  element 
deletions  are: 

ID  =  SKINL  -  for  lower  skin  element 
SKINU  -  for  upper  skin  element 
RIB  -  for  rib  panel 
RCAPU  -  for  upper  rib  cap 
RCAPL  -  for  lower  rib  cap 
SPAR  -  for  spar  panel 
SCAPU  -  for  upper  spar  cap 
SCAPL  -  for  lower  spar  cap 

These  may  be  entered  in  any  combination  or  order 
to  define  the  model  damage. 

(3)  Spanwise  bays  are  numbered  from  root  to  tip;  chordwise 
bays  are  numbered  leading  to  trailing  edges.  For 
spars  and  spar  caps,  INDXI  is  the  spanwise  nodal 

line  affected.  For  ribs  and  rib  caps,  INDXJ  is  the 
chordwise  nodal  line  affected. 

(4)  DX  and  DY  specify  the  center  of  a  sphere  of  radius 
RAD.  All  elements  whose  centroids  lie  within  this 
sphere  will  be  deleted.  All  values  are  given  in  model 
coordinates  supplied  by  the  user.  It  is  assumed  the 
user  will  generate  an  undamaged  model  first,  then  will 
wish  to  generate  a  damaged  model.  DX,  DY  and  RAD 
should  be  easily  obtainable  from  the  nodal  coordinates 
data  of  the  model. 
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2.8.12  Loads  Definition 


1 


(required) 


CARD 

COLUMN 

DATA 

DESCRIPTION 

NOTES 

1 

| 

LOADS 

Literal  'LOADS' 

_ 

1 

NREF 

Number  of  load  data 
specifications 

Literal  'CENTER' 

1 

2 

CENTER 

- 

7-10 

MREF 

Change  center  of  load  plane 
flag 

2 

3 

1-10 

XC 

x-coordinate  of  new  center 
of  load  plane 

3 

11-20 

YC 

y-coordinate  of  new  center 
of  load  plane 

3 

21-30 

ZC 

z-coordinate  of  new  center 
of  load  plane 

3 

i  A.  Class  1  -  Fixture  Loads  -  Specifying  point 

!  loads  to  be  applied. 


CARD 

COLUMN 

DATA 

- 

DESCRIPTION 

NOTES 

4 

1-7 

TYPE 

Literal  'FIXTURE* 

8-10 

(blank) 

11-15 

NLODS 

Number  of  loads 

4 

5 

1-10 

FI 

x-direction  load  value 

5 

11-20 

F2 

y-direction  load  value 

5 

21-30 

F3 

z-direction  load  value 

5 

31-40 

R1 

x  distance  from  structure 
centroid  to  applied  load 

6 

41-50 

R2 

y  distance  from  structure 
centroid  to  applied  load 

6 

51-60 

R3 

z  distance  from  structure 
centroid  to  applied  load 

6 

l/l  ui  ** 


B.  Class  2  -  Test  Loads  -  Specifying  loads  to  be 

distributed  across  all  nodes  on  the  loaded  plane. 


CARD 

COLUMN 

DATA 

DESCRIPTION 

NOTES 

4 

1-4 

TYPE 

Literal  'TEST' 

__ 

5-10 

(blank) 

11-15 

ALPHA 

Angle  (degrees)  of  sweep  of 

wing 

7 

5 

1-10 

FI 

x-direction  load  on  loaded 

plane 

8 

11-20 

F2 

y-direction  load  on  loaded 

plane 

8 

21-30 

F3 

z-direction  load  on  loaded 

plane 

8 

31-40 

F4 

moment  about  x-direction  of 

loaded  plane 

8 

41-50 

F5 

moment  about  y-direction  of 

loaded  plane 

8 

51-60 

F6 

moment  about  z-direction  of 

loaded  plane 

8 

NOTES : 

(1)  Repeat  cards  4  and  5  until  NREF  sets  of  load  data 
have  been  input.  If  NREF  =  0  or  blank,  no  load 

is  to  be  specified  and  card  1  completes  this  item. 

(2)  To  change  the  center  of  the  loaded  plane  from  that 
which  is  calculated  by  the  program  at  MREF  >  0.  If 
MREF  >  0  the  program  will  input  card  3  and  establish 
the  center  of  the  loaded  plane  as  XC,  YC,  ZC.  If 
MREF  =  0  or  blank,  no  change  is  desired  and  card 

3  is  omitted. 

(3)  The  coordinates  XC,  YC,  ZC  specify  the  new  center 
of  the  load  plane.  It  is  advisable  to  have  the 
new  center  lie  within  the  plane  of  loaded  nodes. 

The  effect  of  this  command  will  be  to  alter  the 
load  axis  of  the  wing.  This  is  further  illustrated 

in  Figure  2.8.15.  Default  is  to  use  the  geometric  center 
of  the  wing  tip  load  plane. 
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Figure  2.8.15.  Illustration  of  the  loaded  plane  of  wing  tip. 

12  nodes  (n.-n^-)  are  shown  and  points  A  and 
B.  Point  Aiisxat  the  geometric  center  of  the 
loaded  plane.  Point  B  is  at  the  intersection 
of  the  load  axis  and  the  load  plane.  See 
Figure  2.8.16. 
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(4)  NLODS  represents  the  number  of  points  at  which  a 
load  will  be  applied.  Repeat  card  5  until  NLODS 
load  specifications  have  been  input. 

(5)  Actual  loads  to  be  applied  must  be  resolved  into 
x-,  y-  and  z-components  as  requested  by  parameters 
FI,  F2  and  F3  of  card  5. 

(6)  The  point  at  which  the  load  is  to  be  applied  is 
defined  as  an  x-,  y-  and  z-displacement  from  the 
centroid  of  the  structure.  Card  5  must  be  repeated 
NLODS  times  to  define  all  points  and  their  load 
values. 

(7)  ALPHA  is  the  sweep  angle  of  the  wing  (in  degrees) 
from  the  y-axis  for  purposes  of  identifying  the 
proper  x-,  y-  and  z-components  of  the  moment  loads 
that  can  be  applied.  The  sweep  angle  is  the  angle 
between  the  load  axis  of  the  wing  and  the  y-horizontal 
as  illustrated  in  Figure  2.8.16. 

(8)  Two  methods  exist  for  introducing  loads  in  the 

test  case:  1)  x-,  y-  or  z-direction  loads  components 
and  2)  x- ,  y-  or  z-direction  moments  of  a 

loaded  plane. 

The  loaded  plane  is  defined  as  being  all  the  wing 
tip  nodes  included  in  elements.  The  program  will 
determine  the  center  of  this  loaded  plane  unless 
the  user  overrides  it  with  a  'CENTER'  card.  All 
moment  loads  will  be  resolved  to  the  appropriate 
x-,  y-  and  z-components  and  will  be  distributed 
across  the  load  plane  at  each  node  based  on  the 
area  that  node  encompasses  such  that  there  is  a 
constant  load  distribution  throughout  the  wing  tip. 
This  is  further  explained  in  Section  2.11.1. 


Figure  2.8.16. 


and  Load  Plane 
x 


Alpha  is  the  angle  of  intersection  of  the  load 
axis  with  the  y  axis.  This  angle  is  used  to 
determine  correct  moments  of  loads  applied 
about  the  center  of  the  loaded  plane.  Also 
see  Figure  2.8.15. 
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2.9  MAGNA  LOAD  DECK  CREATION 

The  primary  purpose  of  the  preprocessor  WINGEN  is  the 
generation  of  nodes,  element  connectivities  and  all  other 
pertinent  data  for  the  creation  of  a  card  deck  (or  card-image 
deck)  which  will  be  acceptable  input  to  the  finite  element 
model  analysis  program  MAGNA.  This  section  discusses  the  third 
phase  of  the  WINGEN  operation  process:  the  creation  of  the 
actual  load  deck  data  file  once  the  model  has  been  fully  defined 
and  all  damages  included.  The  user  may  request  that  thin-shell 
elements  be  utilized  on  the  upper  and  lower  skins  to  allow  for 
analysis  of  a  refined  or  damaged  model.  WINGEN  reformats  the 
data  contained  on  the  load  deck  data  file  to  allow  for  the  new 
elements  and  nodes.  If  the  user  responds  to  the  program 
control  directives  that  he  wishes  to  use  shell  elements,  the 
program  will  automatically  generate  a  load  deck.  All  the 
information  necessary  for  a  MAGNA  analysis  execution  is 
contained  in  the  load  deck  which  is  on  file  TAPE11  at  the 
conclusion  of  the  WINGEN  execution.  This  file  may  be  made 
permanent  and  attached  later  to  be  batched  to  input  for  a  MAGNA 
analysis.  Refer  to  Chapter  3  for  details  on  performing  a  MAGNA 
analysis. 

2.9.1  Load  Deck  Creation 

Upon  completion  of  the  model  generation  WINGEN 
prompts  the  user  for  answers  to  several  questions  concerning 
load  deck  creation.  A  sample  of  these  questions  is  included 
in  Figure  2.9.1.  The  user  should  familiarize  himself  with 
the  parameters  requested  concerning  the  actual  analysis  by 
MAGNA.  Reference  1  contains  full  explanations  of  all  MAGNA 
input  data.  Each  question  prompted  by  WINGEN  will  be  explained 
below. 

ENTER  JOB  CARD  ID . : 

This  question  requests  the  user's  computer  ID  symbol,  a  four 
character  alphanumeric  symbol  for  identifying  computer  generated 
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ENTER  JOR  CAM  ID 


Figure  2.9 


_ Ot  CAM  COMMENT 

(PMQDLEN  NO..MANE,SYNROL,ETC. ) . 

MATERIAL  NONLINEARITIES  FLAG  l 
1  -CLASTIC  (LINEAR)  ANALYSIS 
t  -ELASTIC  (NONLINEAR)  ANALYSIS 

ENTER  MATERIAL  HOHLDCARITIES  FLAG . 1 

SEOMETRIC  NONLINEARITIES  FLAG  l 

I  SMALL  DISPLACEMENT  (LINEAR)  ANALYSIS _ 

8  -URGE  DISPLACEMENT  (NONLINEAR)  ANALYSIS 
ENTER  GEOMETRIC  NONLItCARITIES  FLAG....'  1 
ENTER  TME  MURDER  OF  LOAD  STEPS  TO  DE 

PCSrOMMED  IN  TW  SOLUTION . •  1 

ENTER  LOAD  INCREMENT  STEP  SIZE .  I 


.All 


1.  An  example  of  Load  Deck  creation  questions 
prompted  by  WINGEN  upon  generation  of  the 
model  if  the  load  deck  option  was  selected. 


output.  This  symbol  is  required  because  MAGNA  must  be  run  as 
a  batch  job. 

ENTER  JOB  CARD  COMMENT 

(PROBLEM  NO.  NAME,  SYMBOL,  ETC.) . 

The  problem  number  is  the  only  item  required 
here.  A  problem  number  is  a  six  digit  number  preceded  by  a 
letter.  This  number  is  issued  by  the  computing  center  for  cost 
accounting  purposes.  Other  information  may  also.be  of  value 
including  the  individual ' s  name  and  phone  number  in  the  event 
the  output  listing  is  sent  to  the  wrong  location. 

MATERIAL  NONLINEARITIES  FLAG: 

1  -  Elastic  Analysis  (linear) 

2  -  Elastic-Plastic  Analysis  (nonlinear) 

ENTER  MATERIAL  NONLINEARITIES  FLAG . : 

The  material  nonlinearities  flag  selects  whether  the  material  is 
linear  (Elastic  Analysis)  or  nonlinear  (Elastic-Plastic  Analysis) . 
Should  the  user  select  a  nonlinear  material  flag,  the  large 
displacement  analysis  flag  (geometric  nonlinearities  flag) 
will  automatically  be  invoked  by  MAGNA. 

GEOMETRIC  NONLINEARITIES  FLAG: 

1  -  Small  Displacement  Analysis  (linear) 

2  -  Large  Displacement  Analysis  (nonlinear) 

ENTER  GEOMETRIC  NONLINEARITIES  FLAG....: 

The  geometric  nonlinearities  flag  should  be  set  to  small 
displacements  analysis  (1)  for  linear  static  analysis  of  elastic 
materials  or  set  to  large  displacement  analysis  (2)  for 
nonlinear  static  analysis  of  either  elastic  or  plastic  materials. 

ENTER  NUMBER  OF  LOAD  STEPS  TO 

BE  PERFORMED  IN  THE  SOLUTION . : 

ENTER  LOAD  PARAMETER  STEP  SIZE...: 

These  two  input  items  instruct  the  MAGNA  program  how  tc  apply 
the  designated  loads  to  the  model  during  analysis.  The  total 
loads  specified  will  be  divided  by  the  number  of  load  steps 
specified.  The  analysis  will  then  occur  adding  a  load  value 
calculated  to  be  one  load  step  load  (=  total  load/#  steps)  times 
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the  load  parameter  step  size  until  the  full  load  has  been 
applied.  The  user  will  generally  want  both  parameters  equal 
to  1  for  a  linear  analysis.  A  value  of  20  load  steps  and  a 
load  parameter  step  size  of  1  is  usually  adequate  for  a 
nonlinear  analysis. 

WINGEN  will  complete  the  load  deck  data  file 
once  the  above  questions  have  been  answered.  The  load  deck 
will  be  placed  on  local  file'  TAPEll  at  the  conclusion  of  the 
WINGEN  program  execution.  It  is  the  user's  responsibility  to 
make  this  a  permanent  file  for  later  access  or  to  batch  the 
file  for  punching  out  a  card  deck.  Section  3.2  explains  the 
load  deck  data  file  created  and  illustrates  how  alterations 
can  be  effected. 


2.10  WINGEN  PLOTTING  CAPABILITIES 


The  ability  to  represent  finite  element  models  graphically 
has  proved  to  be  invaluable  in  model  development  and  subsequent 
interpretation  of  model  analysis  results.  WINGEN  provides  the 
model  developer  graphics  capability  for  visual  verification 
of  the  model  being  generated.  The  user  has  options  for 
selecting  the  viewing  position  of  the  model,  labeling  nodes 
and/or  elements,  a  limited  zoom  capability,  orthogonal  or 
perspective  viewing,  selected  plotting  of  either  type  3  elements 
or  type  4  elements  and  plotting  of  x-,  y-  and  z-axes.  These 
features  are  necessary  for  adequate  model  verification  by  the 
user  with  the  least  amount  of  knowledge  required  to  execute. 
Should  the  user  require  more  comprehensive  plotting  capabilities 
(e.g.  finer  zoom,  selected  elements,  expand  elements,  etc.) 
the  PLOTBOB  postprocessor  plotting  program  can  be  utilized 
to  view  the  input  data  to  MAGNA  which  WINGEN  places  on  TAPEll. 
Instructions  for  executing  PLOTBOB  are  contained  in  Section  4.2. 

To  enable  the  graphics  mode  the  user  must  respond  '  Y'  to 
the  question  in  the  program  control  directives : 

GENERATE  GRAPHS?  (Y,N) . : 

which  is  requested  at  the  initiation  of  the  program  execution. 
Following  the  node  and  element  connectivity  generation  will  be 
questions  concerning  the  load  deck  generation,  if  one  was 
requested.  The  program  will  then  enter  the  plotting  phase 
and  begins  with  the  question: 

ENTER  THE  CPS  RATE . : 

to  which  the  user  responds  with  the  characters  per  second 
transmission  rate  at  which  he  is  operating  (usually  either  30 
cps  or  120  cps) .  The  program  will  request  responses  to  several 
questions  concerning  eye  position,  minimum  and  maximum  values 
and  orthogonal  or  perspective  viewing.  The  entire  model  will 
be  displayed  on  the  first  plot  but  the  user  may  select  to  plot 
only  selected  element  types  with  or  without  labels  on  subsequent 


plots.  Figures  2.10.1,  2.10.2  and  2.10.3  illustrate  the  various 
plotting  commands  available.  These  commands  are  described 
below. 

2.10.1  Eye  Position 

WINGEN  generates  3-D  finite  element  models.  To 
view  a  3-D  model  the  user  must  imagine  the  model  in  a  3-D 
space  with  the  origin  of  the  space  at  the  model  coordinates 
origin  (see  Figure  2.8.1).  The  user  then  specifies  x-,  y-  and 
z-coordinates  relative  to  the  model  origin  to  place  his  'eye'. 

The  program  then  scales  all  the  data  to  display  the  model  as 
though  the  viewer  was  actually  at  that  relative  distance. 

A  useful  initial  viewing  distance  is  (100,100,100)  for  most 
models . 

2.10.2  Model  Viewing 

A  model  is  represented  as  a  3-D  object  on  a  2-D 
screen.  Occasionally  the  user  may  wish  to  view  only  part  of  that 
object.  The  program  will  display  the  maximum  and  minimum  x,  y 
and  z  values  of  the  model.  The  user  may  alter  these  and 
effectively  move  closer  or  farther  away  from  the  model.  Utilizing 
the  minima  and  maxima  option  the  user  can  isolate  any  portion  of 
the  model  he  wishes  and  then  alter  the  viewing  angle  for  improved 
perspective.  When  altering  the  maxima  and  minima  the  user 
should  input  the  values  as  MINX,  MAXX,  MINY ,  MAXY,  MINZ,  MAXZ. 

Only  those  nodes  which  lie  within  the  selected  range  will  be 
labelled.  (See  Figure  2.10.3). 

2.10.3  Projection  Type 

Two  projection  types  of  view  are  available: 
orthogonal  and  perspective.  An  orthogonal  view  is  one  where  the 
lines  of  the  model  are  projected  perpendicular  to  the  viewing 
plane  or  screen.  A  perspective  view  occurs  when  the  lines  of  the 
model  are  drawn  so  they  intersect  the  viewing  plane  or  screen 
while  being  drawn  to  the  eye  position.  Figure  2.10.4  illustrates 
the  two  viewing  types.  The  perspective  view  tends  to  distort 
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the  structure  such  that  the  elements  closer  to  the  eye  appear 
larger  than  those  farther  away. 

2.10.4  Axes 

The  user  may  wish  to  have  x-,  y-  and  z-axes 
drawn  to  help  with  the  orientation  of  the  model  in  space.  The 
axes  are  drawn  from  the  model  coordinate  origin  to  extend  10% 
beyond  the  maximum  coordinate  distance  in  each  respective 
direction . 

2.10.5  Node  and  Element  Labels 

Finite  elements  become  very  difficult  to 
discriminate  as  models  increase  in  size  and  complexity.  It 
is  useful,  therefore,  to  be  able  to  label  the  nodal  points  and 
elements  to  better  ensure  that  the  model  has  been  properly 
generated.  If  the  user  selects  to  label  the  nodes,  a  node 
number  will  be  written  beside  all  nodes  not  outside  the  minima 
and  maxima  specified  for  the  model.  Plate  (element  type  3)  and 
beam  (element  type  4)  elements  will  also  be  labeled  when 
requested  if  they  do  not  occur  outside  the  minima  and  maxima. 
For  easier  identification  all  plate  element  numbers  are 
prefixed  with  a  '  P'  and  all  beam  element  numbers  are  prefixed 
with  a  'B'.  There  is  no  prefix  for  the  node  numbers. 

2.10.6  Selective  Element  Plotting 

The  user  will  generally  desire  to  view  only  some 
part  of  the  model  in  determining  if  it  has  been  generated 
correctly.  The  program  will  allow  the  user  the  option  of 
plotting  all  of  the  elements,  just  the  plate  elements  (element 
type  3)  or  just  the  beam  elements  (element  type  4) .  Should 
the  user  wish  better  element  plotting  selection  it  is  suggested 
that  he  use  PLOTBOB  plotting  program  with  the  load  deck  file 
created  on  TAPEll.  PLOTBOB  program  utilization  is  discussed 
in  Section  4.2. 

The  user  should  be  careful  to  enter  the  proper 
data  for  each  of  the  plotting  options  above  the  first  time 
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as  there  is  no  means  to  reset  the  values  before  a  plot  is  drawn. 
If  wrong  values  are  entered  the  user  will  have  to  allow  the 
plot  to  be  drawn  then  request  another  graph  and  set  the 
parameters  correctly  the  second  time.  Should  the  user  request 
shell  elements  for  the  upper  and  lower  skin  he  will  need  to 
utilize  PLOTBOB  to  obtain  a  pre-analysis  plot  of  the  structure. 
WINGEN  will  generate  plots  of  the  structure  prior  to  conversion 
to  the  shell  elements  so  the  user  may  utilize  WINGEN  plotting 
to  verify  that  the  model  geometry  is  correct.  The  model  will 
be  converted  following  the  plotting  phase  of  WINGEN. 


Figure  2.10.2.  Sample  plot  from  WINGEN  plotting  option.  All  beam  elements  are 

plotted  and  labeled  (Type  4  elements  or  truss  elements).  All 
options  selected  for  this  plot  are  listed  in  the  upper  right 
quadrant  of  this  plot. 
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Figure  2.10.4a.  Sample  plot  of  replica  wing  structure  with  the  orthogonal  viewing 

option  selected. 


ttNEMTt  ANOTHER  CR*PH?<Y.N)  . 


2.11  SPECIAL  PROGRAM  CONSIDERATIONS 


Several  areas  of  WINGEN  processing  provides  the  program 
user  with  time  savings  convenience  for  such  areas  of  finite 
element  modeling  as  refinement  of  the  element  mesh,  generation 
of  loads  at  the  wing  tip  nodes  and  specification  of  damage 
conditions  of  the  wing.  Each  of  these  program  tools  can  be 
very  useful  to  the  user  if  their  capabilities  and  limitations 
are  fully  understood.  The  remaining  parts  of  this  section  will 
discuss  each  of  these  components  in  detail  to  allow  the  user 
full  use  of  these  features. 

2.11.1  Loads  Application 

WINGEN  provides  the  user  with  two  methods  of 
inputting  loads  to  be  distributed  over  the  wing  tip  nodes. 

The  nature  of  the  loading  situation  in  the  experimental 
condition  being  simulated  should  determine  the  loading  option 
selected  by  the  user  in  the  preprocessor.  The  'TEST'  option 
provided  assumes  the  test  specimen  is  attached  to  the  test 
frame  such  that  the  only  areas  carrying  the  load  are  the  skin 
thicknesses  across  the  top  and  bottom  wing  skins.  This  option 
computes  the  area  of  the  skin  thickness  associated  with  each 
node  and  distributes  the  load  at  each  node  as  a  proportion  of 
the  area  encompassed  by  that  node.  The  'FIXTURE'  option  is 
provided  for  the  situation  where  the  load  is  distributed  across 
the  entire  end  area  of  the  test  specimen.  This  whole  area  is 
used  to  calculate  the  proportion  of  the  load  to  be  distributed 
at  each  node.  Each  of  these  options  is  discussed  further  below. 

The  process  of  applying  loads  to  a  simulated  model 
must  make  assumptions  to  achieve  some  degree  of  realism,  while 
the  model  being  analyzed  experimentally  will  experience  a 
continuous  distribution  of  loading  across  its  loaded  end  a 
simulated  model  will  only  have  loads  applied  at  discrete  points 
on  a  load  plane.  This  situation  is  illustrated  in  Figure 
2.11.1a  and  b  where  a  is  an  experimental  test  frame  generating 
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a  distributed  load  on  a  test  specimen  and  b  is  a  diagram  of  a 

wing  model  with  nodes  at  which  loads  will  be  applied.  The 

survivability/vulnerability  engineer  will  be  primarily  concerned 

with  five  types  of  forces  applied  to  the  wing  section: 

chordwise  shear  (Vc) ,  spanwise  shear  (Vs) ,  torque  (Mt) , 

chordwise  moment  (Me)  and  spanwise  moment  (Ms)  as  illustrated 

in  Figure  2.11.2.  WINGEN  provides  for  all  five  of  these  types 

of  loadings  with  both  of  the  input  options  'TEST'  and  ‘FIXTURE' 

described  in  Section  2.8.12.  In  either  case  all  the  loads  are 

resolved  into  x,  y  and  z  load  magnitudes  which  are  then 

distributed  among  the  nodes  based  on  a  calculated  effective 

nodal  area.  This  effective  nodal  area  is  different  for  the  two 

options  but  the  effect  is  identical  in  that  each  node  on  the 

load  plane  is  assigned  a  fraction  of  the  total  load  based  on  its 

total  effective  area.  The  following  equations  will  help  to 

exolain  this  where  A  .  is  the  total  effective  area  of  the  i*"*1 
*■  ei 

node,  A,p  is  the  total  effective  area  of  all  the  nodes  and  L 
is  the  total  load  applied  to  the  simulated  model  system. 
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i=l 
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ei 


ei  L  =  L . 
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Z  L.  =  L 
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(2.11.1) 

(2.11.2) 

(2.11.3) 


Equation  2.11.2  yields  the  value  L.  which  represents  the  total 
load  applied  to  the  i  node.  The  sum  of  all  the  nodal  loads 
will,  of  course,  equal  the  total  load  applied  as  illustrated 
in  Equation  2.11.3.  Since  the  total  load  is  translated  into 
x,  y  and  z  direction  components.  Equations  2.11.1-3  are  applied 
for  each  component  of  the  load  and  the  sum  of  all  three  components 
equals  the  total  load  applied  to  the  system. 
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FULL  WING 


2.11.1.1  Load  Directions 

All  loads  designations  in  this  manual 

are  such  that  loads  will  be  positive  or  all  resultants  will  be 

positive.  The  chordwise  and  spanwise  shear  forces  (V  ,  V  ) 

c  s 

are  positive  in  the  positive  x  and  z  directions  respectively. 

A  torque  moment  (Mfc)  is  such  that  a  positive  torque  force  will 

yield  a  resultant  force  in  the  positive  y  direction  (the  right- 

hand  rule  will  have  the  thumb  pointing  in  the  positive  y 

direction).  Positive  chordwise  and  spanwise  moments  (M  ,  M  ) 

c  s 

will  likewise  yield  resultant  forces  in  the  positive  z  and 
positive  x  directions  respectively  (e.g.,  the  right-hand  rule 
will  yield  the  thumb  pointing  in  the  positive  z  and  x  directions) . 

2.11.1.2  Test  Loads 

The  'TEST'  loads  input  option  is 

designed  to  input  any  or  all  of  five  types  of  loads  pertinent 

to  the  survivability/vulnerability  engineer.  These  five  loads 

are  chordwise  and  spanwise  shear  loads  (V  ,  V  )  torque  load 

c  s 

(M  )  and  spanwise  and  chordwise  moment  loads  (M  ,  M  ).  These 

u  SO 

are  illustrated  in  Figure  2.11.2.  WINGEN  will  input  these 
loads  and  convert  them  to  x,  y  and  z  components  and  distribute 
them  across  the  load  plane  at  the  nodes  based  on  the 
effective  load  area  of  each  node.  The  resolution  of  the 
input  loads  into  the  various  components  is  a  multi-step  process. 

In  the  first  step  the  loads  are  multiplied  by  the  appropriate 
equations  to  yield  the  valid  forces  components  for  swept 
wings.  Equations  2.11.4  -  2.11.9  illustrate  the  data  input 
procedure.  The  value  ALPHA  is  the  sweep  angle  of  the  wing. 

Fx  -  Vc  (2.11.4) 

Fy  =  Vy  (2.11.5) 

F,  =  Vo  (2.11.6) 

Mx  =  Mx*  COS  (ALPHA) +  My*SIN (ALPHA)  (2.11.7) 
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M  =  M  *  COS  (ALPHA)  -  M  *SIN (ALPHA)  (2.11.8) 

y  y  x 

Mz  =  Mfc  (2.11.9) 

The  program  next  calculates  the  effective  load  area  for  each 
node  and  computes  the  total  effective  load  area.  To  do  this 
the  program  assumes  certain  experimental  conditions  are  being 
simulated.  The  experimental  test  fixture  illustrated  in 
Figure  2.11.1a  is  an  example  of  the  clamp  and  bolt  test 
apparatus  this  program  was  designed  to  simulate.  This  type  of 
apparatus  is  characterized  by  attaching  the  main  clamp  of  the 
loading  end  to  the  skin  panels  of  the  test  specimen.  The  net 
effect  of  this  procedure  is  that  the  skin  area  alone  carries 
the  load  distribution  from  the  actuators  applying  the  load  to 
the  remaining  structure.  In  reproducing  the  experimental 
situation  a  model  simulation  must  distribute  the  loads  to  the 
individual  nodes  in  a  fashion  similar  to  the  load  distribution 
in  the  experimental  situation.  To  accomplish  this  the  ’TEST* 
option  determines  the  total  area  of  the  skin  coincident  with 
the  load  plane  and  apportions  it  to  each  individual  node.  This 
is  accomplished  by  dividing  each  bay  in  half  and  combining  the 
two  halves  to  either  side  of  a  node  to  yield  a  length  1^  for 
each  node.  This  length,  1,  is  then  multiplied  by  the  thickness 
of  the  skin  to  generate  the  effective  load  area  of  the  node. 
This  process  is  illustrated  in  Figure  2 . 11 . 3 .  All  the 
effective  nodal  load  areas  are  added  to  yield  the  total 
effective  load  area.  The  loads  are  then  distributed  over  the 
nodes.  Each  node  is  taken  in  turn  through  equations  2.11.10  - 
2.11.12  until  all  nodes  in  the  load  plane  have  had  the  full 
load  distributed  over  them. 


IX 


F  M  *A  .  *R 
x  ,  y  ei  z 

n  A. 


F  M  *A  .  *R  M  *A  . *R 
L  _JL  x  ei  z  y  ei  x 

iy  n  A.  A 


(2.11.10) 

(2.11.11) 
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3.  Test  option  loads  distribution.  Each  bay  is 
divided  in  half  and  each  half  bay  is  associated 
with  the  node  adjacent  to  it.  The  length  of  this 
upper  skin  associated  with  each  node  is  then 
multiplied  by  the  skin  thickness  to  obtain  an 
effective  area  for  the  node,  Ag.  Ae  is  divided 
by  the  total  effective  area,  Am,  to  find  the 
fraction  of  the  load  to  be  applied  at  that  node. 
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(2.11.12) 


F  M  *A  . *R 

l  =  —  -  y  el  x 

iz  n  A. 


In  these  equations,  L.  represents  the 
load  applied  to  the  i1"  node  in  the  x,  y  or  z  direction;  F 
represents  the  total  shear  force  in  the  x,  y  or  z  direction;  n 
is  the  number  of  nodes  to  which  loads  will  be  applied;  M  is  the 
moment  in  the  x  or  y  directions;  Afc  is  the  total  effective  node 
load  area  for  loads  application;  A  .  is  the  effective  node  load 
area  for  the  i  node;  and  R  is  the  radius  from  the  center  of 
the  load  plane  to  the  node  in  the  x  or  z  direction.  The  load 
plane  is  defined  by  WINGEN  to  be  perpendicular  to  the  y  axis 
and  all  values  affected  are  adjusted  accordingly  by  the 
following  equations: 


x±  =  xi*COS (THETA)  -  Y^*SIN (THETA)  (2.11.13) 


y.  =  y. *COS (THETA)  +  x. *SIN (THETA)  (2.11.14) 

11  1 

where  x^  and  y.^  represent  the  x  and  y  coordinates  of  the  ith 
node  and  theta  is  the  angle  between  the  perpendicular  load 
plane  and  the  wing  tip  that  contains  all  the  nodes  to  which 
the  loads  are  to  be  applied.  WINGEN  calculates  the  angle 
theta  and  does  all  transformations. 

2.11.1.3  Fixture  Loads 

An  alternate  method  for  detailing 
load  input  to  the  nodes  is  with  the  FIXTURE  option.  This 
method  allows  the  user  to  specify  the  loads  in  terms  of  a 
magnitude  of  the  load  and  the  x,  y  and  z  distances  from  the 
center  of  the  loaded  plane  to  where  the  load  should  be  applied. 
WINGEN  will  take  this  load  and  resolve  it  into  its  components 
and  distribute  it  across  the  nodes  of  the  load  plane.  As  with 
the  'TEST'  case  the  center  of  the  load  plane  may  be  adjusted 
to  be  wherever  the  user  wishes.  The  method  of  determining  the 
area  over  which  the  loads  should  be  applied  is  slightly 
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different  than  discussed  above  for  the  'TEST'  option.  In  the 
'TEST*  option  the  effective  node  area  for  loads  application 
was  calculated  using  the  thickness  of  the  skin  and  the  length 
of  the  skin  'affected*  by  each  node.  The  'FIXTURE'  option  is 
based  on  a  different  experimental  test  concept  which  has  the 
entire  load  plane  active  in  distributing  the  load  through  the 
structure  via  some  type  of  end  plate  or  other  similar  device. 

To  compute  the  effective  node  area  the  program  uses  the  full 
depth  of  the  wing  and  computes  a  real  area  bounded  by  the  node 
as  illustrated  in  Figure  2.11.4.  To  obtain  this  area  each  bay 
is  divided  into  halves  and  each  half  is  associated  with  the 
nearest  node.  The  total  of  each  of  these  half  areas  equals 
the  total  effective  area  of  that  node.  The  sum  of  all  of  these 
area  yields  the  total  effective  area  for  the  load  plane. 
Equations  2.11.15  -  2.11.20  show  how  the  input  loads  are 
resolved  to  standard  forces : 


F 

V 

(2.11.15) 

X 

X 

F 

= 

V 

(2.11.16) 

y 

y 

F 

= 

V 

(2.11.17) 

z 

z 

M 

X 

= 

V«y 

-  V  *R 

y  z 

(2.11.18) 

M 

— 

V  *R 

-  V  *R 

(2.11.19) 

y 

X  z 

Z  X 

M 

= 

V  *R 

-  V  *R 

(2.11.20) 

z 

y  * 

X  y 

when  R  ,  R  and  R  are  the  distances  from  the  center  of  the 
x  y  z 

loaded  plane  to  the  point  where  the  load  is  applied.  WINGEN 
then  applies  Equations  2.11.10  -  2.11.14  to  the  forces 
calculated  and  yields  a  distribution  of  loads  across  the  load 
plane  nodes. 

2.11.2  Damage  and  Modify  Directives 

WINGEN  provides  the  user  with  a  set  of  options  to 
specify  alterations  to  the  basic  wing  model  geometry.  These 
options,  described  in  Sections  2.8.8  and  2.8.11,  adhere  to  basic 


fundamental  criteria  for  finite  element  modeling  such  that  no 
elements  or  nodes  are  in  violation  of  the  basic  modeling 
criteria.  The  first  of  these  criteria  is  that  no  element 
should  intersect  more  than  one  element  on  a  side  without  special 
restrictions.  An  example  of  this  is  given  in  Figure  2.11.5 
where  an  additional  rib  has  been  added  to  the  T-38  wing  model. 

In  order  to  accommodate  this  violation  certain  linear  constraints 
must  be  defined  in  the  input  data  file  (load  deck)  before  an 
analysis  can  be  performed  by  MAGNA.  These  constraint  equations 
are  disected  to  the  various  components  and  entered  in  the 
data  file  as  described  in  Sections  8.6  and  8.7  of  Reference  1. 
Table  2.11.1  lists  the  two  sets  of  constraint  data  entered  in 
the  program  for  the  above  problem. 

TABLE  2.11.1 

LINEAR  CONSTRAINT  EQUATION  TERMS 


EQUATION 


lower 

skin 

constraint 

equation 

terms 


upper 

skin 

constraint 

equation 

terms 


NODE  # 

DIRECTION 

MULTIPLIER 

I 

53 

1 

1. 

39 

1 

-.4102 

67 

1 

-.5898 

53 

2 

1. 

39  1 

2 

-.4102 

67 

2 

-.5898 

53 

3 

1. 

39 

3 

-.4102 

67 

3 

-.5898 

54 

1 

1. 

40 

1 

-.4102 

68 

1 

-.5898 

54 

2 

1. 

40 

2 

-.4102 

68 

2 

-.5898 

54 

3 

1. 

40 

3 

-.4102 

68 

3 

-.5898 
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These  equations  serve  to  effectively  prevent  the  elements  a  and 
b  of  Figure  2.11.6  from  displacing  in  a  manner  different  than 
element  c  under  loading  conditions. 


Figure  2.11.6. 


An  example  of  an  incompatible  element  alignment 
between  elements  a,  b  and  c,  above.  (This  is 
presented  in  Figure  2.11.5.) 


This  is  accomplished  by  essentially  specifying  that  node  n,. 
will  not  be  displaced  any  more  than  a  linear  displacement  of 
nodes  n 2  and  The  terms  listed  in  Table  2.11.1  would  relate 

node  53  with  ncf  node  39  with  n_  and  node  67  with  n_  where 
node  n^  is  the  constrained  node.  To  obtain  the  multipliers 
for  this  case  (remember  this  is  a  skin  element  constraint)  the 
length  of  the  bay  is  taken  and  a  proportion  is  established  for 
the  relative  constraining  requirements  placed  on  each  node. 

The  multipliers  of  nodes  n 2  and  n^  are  given  by  equations 
2.11.21  and  2.11.22.  The  multiplier  for  node  n5  is  1.0  as 
demonstrated  by  equation  2.11.23. 


11+12 


h+12 


(2.11.21) 


(2.11.22) 


11+12  11+12 


(2.11.23) 
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For  these  equations  1^  represents  the  distance  between  nodes 
n7  and  n^  while  1^  represents  the  distance  between  nodes  n ^  and 
n^  as  illustrated  in  Figure  2.11.6. 

Any  alterations  made  to  the  model  element  to  yield 

a  structure  which  is  perhaps  more  accurate  but  that  contains 

elements  which  are  not  uniform  throughout  the  structure  must 

have  linear  constraints  described  for  the  non-uniform  union 

of  the  incompatible  element  types.  To  obtain  a  structure  such 

as  the  T-38  wing  illustrated  in  Figure  2.11.5  the  user  had  to 

alter  the  load  deck  file  by  appending  the  additional  four  nodes 

to  the  node  coordinate  list  and  add  the  element  connectivity 

for  the  new  rib  section.  In  addition  new  element  connectivities 

were  given  for  the  plate  elements  (upper  and  lower  skins)  and 

two  new  plate  elements  were  generated  outboard  of  the  new  spar. 

Linear  constraints  were  also  added  as  detailed  above  to  ensure 

continuity  of  the  stresses  and  displacements  between  elements. 

An  alternate  method  to  accomplish  the  addition  of  this  short  rib 

would  have  been  to  specify  the  definition  of  a  normal  rib  when 

the  model  was  created.  Once  the  load  deck  file  was  generated 

the  user  could  essentially  remove  the  unneeded  rib  elements  by 

defining  a  new  material  property  code  with  identical  values 

for  all  components  as  the  original  property  code  but  with  a 

20 

considerably  higher  yield  stress  value  (1.1*10  ).  To  prevent 

the  elements  from  being  utilized  in  the  analysis  a  very  small 

-20 

thickness  for  the  element  should  be  specified  (1.0*10  ) 

along  with  the  change  to  the  new  material  property  code.  These 
two  changes  to  an  element  will  remove  it  from  any  effective 
role  in  the  analysis.  This  latter  method  is  very  convenient 
if  a  large  file  is  being  generated  by  the  preprocessor  that 
would  take  a  considerable  period  of  time  to  alter  or  a  large 
number  of  short  ribs  need  to  be  defined.  This  second  method 
also  eliminates  the  requirements  of  inputting  the  linear 
constraints  although  model  plots  will  show  the  additional 
rib  which  may  or  may  not  be  important  to  the  project  engineer. 


A  second  important  criteria  for  finite 
element  modeling  is  that  all  nodes  defined  for  a  model  must 
either  be  utilized  in  element  connectivities  or  must  be 
constrained  in  all  directions.  Any  superfluous  nodes  generated 
by  the  program  will  be  so  constrained  but  if  the  user  should 
add  any  additional  nodes  or  delete  any  elements  which  would 
result  in  nodes  not  being  utilized  in  elements,  then  the  user 
must  be  careful  to  add  these  to  the  list  of  nodes  to  be 
constrained.  Figure  2.11.7  shows  a  portion  of  the  constraint 
data  contained  on  the  WINGEN  Load  Deck  after  user  modification 
to  constrain  nodes  for  an  additional  rib  section  (Element  number 
80  in  Figure  2.11.5)  inserted  into  the  T-38  wing  model.  The 
first  three  lines  of  this  figure  contain  information  about  the 
last  element  type  defined  and  end-of-element  processing  cards 
(2  blank  cards) .  The  fourth  card  contains  information  defining 
the  boundary  conditions  for  the  model.  The  first  three  values 
of  the  fourth  card  specify  the  number  of  types  1,  2  and  3 
constraints  to  be  defined  while  the  fourth  and  fifth  variables 
define  the  number  of  linear  constraint  equations  and  the 
maximum  number  of  terms  per  constraint  equation,  respectively. 

Two  lines  are  required  for  each  type  1  or  2  constraint 
definition.  Since  three  type  1  and  2  constraints  are  specified, 
the  next  six  lines  tell  MAGNA  that  nodes  1  through  16,  25  through 
28  and  43  through  52  are  all  constrained  to  zero  displacements 
in  all  directions  (for  designation  of  directions  MAGNA  interprets 
1  as  x,  2  as  y  and  3  as  z  based  on  a  Cartesian  coordinate 
system) .  Following  the  boundary  constraints  are  the  linear 
constraints.  The  eleventh  line  tells  how  many  terms  are  in  the 
next  equation  and  the  program  then  reads  those  terms.  Two 
constraint  equations  are  defined  here,  with  nine  terms  per 
equation.  The  last  line  indicates  the  number  of  load  cases  to 
follow.  The  linear  constraints  illustrated  here  are  derived 
from  the  element  situation  illustrated  in  Figure  2.11.6.  Further 
information  may  be  obtained  from  the  MAGNA  User's  Manual 


(Reference  1).  Nodes  that  are  not  utilized  in  element 
generation  must  be  constrained  using  type  1  or  2  constraints. 

The  Damage  and  Modify  options  are 
deletion  directives  in  that  they  specify  the  removal  of  various 
elements.  With  practice  the  user  can  become  adept  at  defining 
structures  and  specifying  damage  or  modification  to  very 
particular  areas  by  combining  several  of  the  damage/modification 
directives  with  the  refinement  options.  Additional  model 
editing  as  described  above  should  allow  the  definition  of 
any  wing  structure  likely  to  require  modeling. 
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SECTION  3 

FINITE  ELEMENT  ANALYSIS  PROGRAM  MAGNA 

3.1  INTRODUCTION 

The  backbone  of  any  finite  element  modeling  and  analysis 
is  the  program  that  actually  performs  the  analysis.  This 
application  makes  use  of  the  MAGNA  finite  element  analysis 
program.  MAGNA  was  developed  as  a  finite  element  program  for 
the  materially  and  geometrically  nonlinear  analysis  of  3-D 
structures  subjected  to  static  and  transient  loading  conditions. 
The  requirements  of  this  application  call  for  a  program  with 
restart  capabilities,  elastic-plastic  (material  nonlinear 
solution)  capability,  large  deflection  (geometric  nonlinear 
solution)  capability  and  an  appropriate  finite  element  library. 
MAGNA  meets  these  needs  easily.  The  requirements  of  wing  flight 
load  simulation  do  not  utilize  a  number  of  features  incorporated 
into  MAGNA  thus  these  features  will  not  be  discussed  in  this 
report.  Reference  1  provides  a  detailed  discussion  of  MAGNA 
and  how  to  utilize  it  for  finite  element  analysis. 

Illustrated  below  is  the  basic  operating  scheme  for 

MAGNA. 


This  section  of  the  report  will  deal  with  these  three  components 
of  the  MAGNA  environment.  The  user  must  first  define  an  input 
file  of  data  which  defines  the  nature  of  the  problem  to  be 
solved  and  provides  all  of  the  information  in  terms  of  model 
definition,  boundary  conditions  and  loads  to  be  applied.  The 
generation  of  such  a  file  can  be  very  time  consuming  and 


3.1 


laborious.  As  a  result  the  preprocessor  program  WINGEN  was 
developed  to  alleviate  the  burden  of  generating  much  of  the 
file  required  by  MAGNA.  The  file  discussed  in  Section  3.2  is 
the  file  output  by  WINGEN  for  MAGNA  analysis.  Once  the 
analysis  is  complete  MAGNA  will  produce  two  output  files.  The 
first  will  be  a  printed  tab  of  the  analysis  results  in  a 
legible  form  for  the  engineer  to  consult;  the  second  output  will 
be  a  disk  file  or  WINGMPOST  file  for  use  by  the  postprocessor 
programs  CONTOUR  and  PLOTBOB.  The  WINGMPOST  file  will  contain 
essentially  the  same  data  as  is  printed  on  the  tab  but  in 
different  formats. 


3.2 


3.2  PREPROCESSOR  OUTPUT  FILE/ANALYSIS  INPUT  DATA  FILE 

The  sole  function  of  the  preprocessor  program  WINGEN 
is  to  produce  a  complete  data  file  (load  deck)  of  a  finite 
element  model  for  analysis  by  MAGNA  from  simplified  input 
directives.  Upon  input  of  the  load  deck  parameters  WINGEN 
will  place  all  necessary  data  on  local  file  TAPE11  for  the  user's 
access  following  the  preprocessor  program  execution.  Included 
on  TAPE11  are  necessary  items  such  as  job  control  cards,  load 
data,  nodal  coordinates,  element  connectivities,  MAGNA  control 
parameters,  etc.  Figure  3.2.1  is  a  complete  listing  of  a  sample 
load  deck  created  by  WINGEN  (for  input  to  MAGNA) .  Should  the 
user  require  more  detailed  information  than  is  provided  below, 
he  should  consult  Reference  1. 

Following  the  execution  of  WINGEN,  the  user  will  find  all 
the  load  deck  information  described  in  this  section  on  TAPEll. 
This  is  a  local  file  and  must  be  saved  by  the  user  if  he  desires 
to  use  it  again  in  the  future.  The  file  may  be  saved  with 
standard  Request,  Copy  and  Catalog  commands  as  described  in  the 
CDC  NOS/BE  operating  manuals  (Reference  2). 

3.2.1  Job  Control  Cards 

The  data  stored  on  TAPEll  consists  of  3  data 

records : 

(1)  the  job  control  cards; 

(2)  program  directive  and  model  definitions;  and 

(3)  space  requirements  for  the  analysis  program 

to  run. 

The  job  control  cards  are  standard  and  are  fully  explained  in 
the  CDC  reference  manual  (Reference  2).  Figure  3.2.2 
illustrates  the  first  data  record  in  the  load  deck. 

This  record  contains  the  MAGNA  job  control  cards. 
The  cards  are  numbered  from  100  to  170  for  discussion  below. 

Line  100  is  the  job  card  and  contains  information  for  the 
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Figure  3.2.2.  Load  deck  job  control  language.  Illustrated  here  are  the  job 
control  cards  provided  by  WINGEN  to  execute  the  MAGNA  analysis 


computer  to  know  who  to  charge  and  where  to  send  the  output 
results  as  well  as  information  on  total  cpu  time,  input/output 
limits  and  central  memory  use  limits.  These  defaults  for 
time,  I/O  and  central  memory  should  be  sufficient  for  any  model 
being  investigated.  Should  you  obtain  a  diagnostic  in  any  of 
these  three  areas  indicating  the  value  (time)  has  been  exceeded, 
simply  increase  the  value  (time)  and  rerun  the  program.  (In 
most  cases  it  is  possible  to  restart  the  MAGNA  program.  Consult 
Reference  1  on  the  use  of  MAGNA  if  this  becomes  a  frequent  or 
costly  problem.)  Lines  110  -  120  set  initial  execution  conditions 
for  MAGNA  and  place  the  model  analysis  input  data  for  MAGNA  on 
local  input  file  TAPE5.  Line  130  provides  for  the  capability 
to  make  the  output  file  from  the  analysis  (the  MPOST  file)  a 
permanent  file;  the  command  for  which  is  executed  at  line  160. 
Lines  140  and  150  actually  execute  MAGNA.  All  postprocessing 
output  from  the  program  will  be  made  a  permanent  file  with  the 
name  "WINGMPOST".  There  is  a  limit  of  five  such  files  with  the 
same  name  (cycles)  so  the  user  must  take  caution  to  make  sure 
there  are  no  more  than  four  cycles  present  under  the  name 
"WINGMPOST"  prior  to  batching  the  job  to  the  input  queue. 

3.2.2  Finite  Element  Model  Data 

The  second  record  on  the  load  deck  file  contains 
all  the  information  about  the  model  needed  by  MAGNA  to  perform 
the  indicated  analysis.  Figure  3.2.1  contains  parts  1-10 
which  are  discussed  briefly  below.  More  detailed  information 
about  these  sections  may  be  obtained  from  Reference  1. 

Part  one  is  the  problem  title.  There  are  three 
lines  available  for  the  user.  The  first  line  contains  the  title 
card  entered  by  the  user  to  the  preprocessor;  the  subsequent 
two  cards  are  blank  and  may  be  altered  to  suit  particular  user 
needs.  All  three  title  lines  are  placed  on  CONTOUR 
postprocessor  plots  when  the  LABE  command  is  utilized. 
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Part  two  defines  the  MAGNA  program  controlling 
directives.  These  cards  may  vary  considerably  depending  on  the 
type  of  analysis.  Some  of  these  values  are  the  response  made 
by  the  user  to  the  load  deck  concerning  the  nature  of  the 
analysis  which  are  asked  by  WINGEN.  The  additional  parameters 
and  controlling  cards  are  supplied  by  the  program  as  default 
values  for  the  user-selected  options. 

Part  three  contains  the  nodal  coordinate  data. 
This  section  is  identified  by  the  header  keyword  "COORDINATES" 
and  an  associated  number  indicating  the  number  of  nodes 
defined  for  this  particular  model.  Following  the  header  card 
are  enough  cards  to  define  all  the  model  nodes.  Each  card 
contains  a  node  number  and  the  x,  y,  and  z  coordinates  for 
that  node  number . 

Part  four  contains  header  card  information  and 
material  property  data  for  an  element  type  (element  type  3  in 
this  case) .  The  first  card  of  this  section  defines  the 
following  section  of  data  to  be  element  type  3;  that  there  are 
two  material  property  codes  and  the  total  number  of  type  3 
elements  to  be  defined.  The  next  two  lines  are  the  material 
property  codes  discussed  earlier  in  the  report.  If  the  user 
is  interested  in  utilizing  additional  property  codes,  the 
total  number  of  material  property  codes  listed  on  the  element 
header  card  (parameter  two)  must  be  increased  by  1  for  each 
such  addition  and  the  necessary  values  (discussed  in  Section 
2.3  -  MODEL  DEFINITION)  to  define  the  material  properties  must 
be  placed  after  the  first  two  material  properties  and  before 
the  first  element  connectivity  data. 

Part  five  lists  all  the  element  numbers,  element 
subtype,  material  property  codes,  nodal  connectivities  and 
element  thicknesses,  respectively,  for  the  element  type 
discussed  in  part  four  above.  The  user  should  be  aware  that 
element  subtype  1  of  element  type  3  will  be  converted  to  shell 
elements  if  the  shell  elements  option  is  implemented.  The 


shell  elements  are  3-D  elements  of  which  two  dimensions  are 
defined  by  the  original  membrane  element  and  the  third 
dimension  added  is  composed  of  the  element  thickness.  Four 
new  nodes  are  generated  for  each  membrane  element  to  convert 
it  to  a  shell  element.  Each  of  these  new  nodes  is  given 
coordinates  (see  part  three)  to  place  it  the  thickness  above 
the  upper  skin  elements  or  the  thickness  below  the  lower  skin 
elements.  The  user  may  alter  these  thicknesses  or  newly 
generated  coordinates  (z  components)  to  more  accurately  reflect 
the  skin  depth  variations. 

Part  six  contains  information  on  element  type  4 . 
As  with  part  five,  there  exists  a  header  card  with  the 
element  type  (type  4  in  this  case) ,  the  number  of  material 
property  codes  and  the  number  of  type  4  elements  to  be  defined. 
This  header  card  is  followed  by  the  material  property  codes 
definitions  cards  discussed  earlier  in  this  report  (Section 
2.3  -  MODEL  DEFINITION).  Alterations  may  be  made  to  the 
material  property  codes  as  discussed  in  part  four. 

Part  seven  contains  the  data  defining  the  type 
4  elements:  element  number,  material  property  code,  nodal 
connectivities  and  element  cross-sectional  area,  respectively. 

Part  eight  consists  of  boundary  conditions  and 
other  constraint  information.  The  preprocessor  program  will 
automatically  constrain  all  the  nodes  on  the  root  chord  of 
the  wing  model  in  the  x,  y,  and  z  directions  as  this  end  will 
generally  be  fixed  by  the  constraints  of  the  experiment.  Other 
constraints  will  be  generated  by  WINGEN  if  the  model  is  refined 
with  a  "THICK"  directive  generating  a  new  depthwise  level  of 
nodes.  The  user  may  wish  to  add  further  constraints  depending 
on  experimental  designs.  Should  this  occur,  he  is  referred  to 
Section  8.6  -  BOUNDARY  CONDITIONS  in  Reference  1  concerning 
MAGNA  input  data  and  Section  2.11  of  this  report. 

Part  nine  is  present  only  for  nonlinear  analysis 
runs  and  contains  information  for  loading  curves  for  nonlinear 


static  analysis  and  uniaxial  stress-strain  data  for  elastic- 
plastic  materials.  These  curves  are  provided  bv  WINGEN  and 
should  not  require  any  changes  or  modifications.  Figure  3.2.1 
illustrates  a  linear  analysis  and,  therefore,  does  not  contain 
this  part. 

Part  ten  is  the  definition  of  the  load  data.  This 
information  is  supplied  as  the  number  of  load  cases  for  this 
analysis  run  followed  by  as  many  cards  as  necessary  to  define 
all  loads  to  input.  Only  one  load  case  is  permitted  for  each 
nonlinear  analysis;  however,  multiple  load  cases  are  permitted 
for  linear  analyses.  The  data  defined  is  the  load  case  number, 
the  node  number  to  which  the  load  is  applied,  the  x,  y,  or  z 
direction  that  the  load  is  to  affect,  and  the  magnitude  of  the  load. 

3.2.3  MAGNA  Space  Parameters 

This  completes  the  model  data  required  for  an 
analysis  run.  One  record  trails  the  model  data.  This  record 
defines  the  amount  of  space  required  by  MAGNA  to  actually  execute 
the  analysis.  The  user  should  not  have  to  alter  these  WINGEN 
defined  values. 

All  this  data  is  supplied  by  WINGEN  at  the 
conclusion  of  the  run  on  TAPE11.  The  user  may  wish  to  make 
alterations  to  this  data  prior  to  making  it  a  permanent  file. 

Nodes  and  elements  may  be  added  or  deleted  but  the  user  must  be 
careful  to  resequence  all  nodes  or  elements  and  must  make  all 
references  to  any  newly  labeled  nodes  correspond  to  the  new  node 
number.  Following  this  section  is  an  additional  sample  load 
deck  created  by  WINGEN  (Figure  3.2.3)  illustrating  a  nonlinear 
analysis  input. 

3.2.4  MAGNA  Program  Execution 

Due  to  the  large  space  requirements  of  analyzing 
finite  element  models,  MAGNA  has  been  designed  as  a  Batch 
execution  job.  This  means  the  user  has  no  control  over  the 
program  execution  once  the  job  is  submitted.  To  initiate  a 
MAGNA  analysis  job  do  one  of  the  following; 
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a)  Submit  a  punched  copy  of  the  load  deck  created 
on  TAPE11  by  WINGEN  to  a  card  reader. 

or  b)  if  TAPE11  was  made  a  permanent  file,  attach 
it  as  a  local  file  and  execute  a  BATCH 
command  as  follows: 

ATTACH, A, WINGLOADECK, . . . 

BATCH, A, INPUT 

These  two  commands  will  result  in  an  analysis 
being  performed . 

Regardless  of  the  method  selected,  the  results 
will  be  returned  as  is  customary  for  remote  job  processing.  The 
following  method  may  be  employed  to  make  certain  a  load  deck 
is  made  a  permanent  file  following  the  WINGEN  program 
execution: 

REQUEST, P,*PF 

REWIND, TAPE11,P 

COPYBF , TAPE 11 , P 

CATALOG, P, WINGLOADECK, RP=300, . . . 

This  will  then  allow  the  user  to  execute  the 
Batch  command  for  a  MAGNA  job  run. 
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Figure  3.2.3.  (continued) 


3.3  MAGNA  ANALYSIS 


The  survivability/vulnerability  engineer  is  primarily 
concerned  with  analysis  results  of  structures  and  determining 
whether  a  structure  is  capable  of  supporting  the  desired  load 
given  certain  extenuating  conditions  such  as  inflicted  damage 
on  the  structure  or  unusual  load  conditions.  Other  information 
such  as  alternate  load  paths  present  on  the  structure  following 
damage,  failure  loads  under  various  conditions  and  confirmation 
of  wing  structure  response  to  experimental  loading  conditions 
may  also  be  of  interest.  The  finite  element  analysis  program 
MAGNA  will  provide  the  user  with  the  type  of  information  when 
the  appropriate  finite  element  data  is  supplied.  MAGNA 
(Materially  And  Geometrically  Nonlinear  Analysis)  is  a  large 
scale  computer  program  for  the  static  analysis  of  complex 
three-dimensional  engineering  structures.  Isoparametric 
modeling  techniques  are  coupled  with  state-of-the-art  numerical 
analysis  and  programming  methods  to  provide  accurate  and 
efficient  solutions  for  large  problems  involving  highly 
nonlinear  response. 

The  modeling  capabilities  for  MAGNA  include  structural 
elements  for  truss  members  (bars) ,  plane  stress  and  plane 
strain  sections  (plates),  "shear  panels",  general  three- 
dimensional  solid  and  thin  plates  and  shells.  All  the  finite 
elements  are  arbitrarily  oriented  and  are  fully  compatible 
in  three-dimensional  space.  Each  of  the  finite  elements  in 
MAGNA  includes  the  effects  of  full  geometrical  nonlinearities 
(large  displacements,  large  strains),  using  a  Lagrangian 
(fixed  reference)  description  of  motion.  Material  nonlinearities, 
in  the  form  of  elastic-plastic  behavior;  are  analyzed  using  a 
subincremental  strategy  which  minimizes  the  error  in  following 
the  material  stress-strain  curve.  Other  features  and 
theoretical  developments  are  contained  in  the  MAGNA  User's 
Manual,  Reference  1.  The  user  is  referred  to  that  manual  for  a 
detailed  discussion  of  MAGNA. 
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3.4  MAGNA  OUTPUT  FILE 


No  structural  analysis  program  is  useful  if  the  engineer 
cannot  locate  or  interpret  the  data  generated.  The  nature  of 
finite  element  analysis  is  that  considerable  amounts  of  printed 
output  is  generated  making  it  difficult  to  isolate  what  is 
important  in  determining  answers  to  problems.  Plots  are  used 
frequently  in  this  area  to  isolate  that  data  which  is  most 
relevant  and  present  it  in  a  form  readily  comprehendable  by  the 
engineer.  Printed  tabs  of  the  output  data  are  also  convenient 
for  a  more  thorough  evaluation  of  the  structural  analysis. 

Figure  3.4.1  is  a  listing  of  the  major  components  of  a  nonlinear 
analysis  output.  A  linear  static  analysis  will  yield  a  similar 
output  listing.  The  differences  between  the  two  are  discussed  in 
each  section  below.  The  reader  will  note  that  the  sample 
analysis  run  in  Figure  3.4.1  is  divided  into  fourteen  parts. 

Each  of  these  parts  will  be  discussed  below. 

3.4.1  Header  for  MAGNA  Program  (Figure  3.4.1a) 

This  banner  identifies  a  MAGNA  structural  analysis 
run  and  provides  the  user  with  important  current  information 
on  who  to  contact  if  problems  arise  and  when  the  last  system 
updates  were  effected.  A  page  of  system  notes  is  also 
included  containing  information  that  may  or  may  not  be 
important.  If  there  are  any  problems  with  MAGNA,  contact 
the  individual  named  on  these  pages. 

3.4.2  Program  Space  Requirements  (Figure  3.4.1b) 

Part  two  is  a  summary  of  the  working  areas  in  the 
program  required  for  the  analysis  of  the  structural  problem. 

The  user  should  not  have  to  be  concerned  with  this  unless 
significantly  complex  models  are  analyzed  or  models  with 
considerable  refinements  are  submitted  for  analysis  by  MAGNA. 

If  this  is  the  case  and  problems  arise  consult  the  individual 
listed  on  the  Header  Page. 
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3.4.3  Main  Program  Listing  (Figure  3.4.1c) 

This  part  contains  a  listing  of  the  main  program 
segment  for  MAGNA.  This  is  of  no  value  to  the  engineer  or 
user. 

3.4.4  Input  Data  File  (Figure  3. 4. Id) 

MAGNA  will  always  copy  the  input  data  file  to  the 
output  file  prior  to  executing  the  analysis.  Part  four  is  the 
complete  load  deck  generated  by  the  preprocessor  WINGEN  on 
TAPEll .  This  file  is  discussed  in  Section  3.2  and  in  detail  in 
Reference  1. 

3.4.5  MAGNA  Options  Selected  (Figure  3.4.1e) 

MAGNA  is  a  general  purpose  finite  element  program 
and  because  of  that  has  a  large  number  of  available  solution 
options  to  accommodate  a  wide  range  of  requirements.  Part 
five  lists  the  solution  options  selected  for  the  current 
analysis.  Only  a  very  few  of  these  options  will  be  pertinent 
for  simulated  wing  loads  analysis.  The  reader  may  refer  to 
Section  2.9  for  a  discussion  of  the  load  deck  creation  and 
options  selection  available. 

3.4.6  Nodal  Coordinates  (Figure  3. 4. If) 

Part  six  provides  a  listing  of  the  input  nodal 
coordinates  for  the  model  finite  elements  definition. 

3.4.7  Element  Connectivities  (Figure  3.4.1g) 

Each  model  is  composed  of  several  types  of 
elements  including  the  truss  (bar) ,  plane  stress  (membrane 
plate) ,  shear  panel  (shear  plate)  and  shell  elements.  This 
section  lists  all  the  input  values  and  the  material  property 
codes  defined  for  each  element  type. 

3.4.8  Constraint  Data  (Figure  3.4.1h) 

Part  eight  is  concerned  with  model  displacement 
restrictions  or  contraints.  The  program  provides  a  summary 
of  the  boundary  conditions  imposed  on  the  model  both  by  the 
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experiment  and  due  to  damage/modify  specifications.  Following 
the  boundary  conditions  is  a  nodal  variable  table  providing 
a  listing  of  the  model  degrees  of  freedom  at  each  node.  The 
linear  constraints  imposed  on  the  model  are  then  listed 
followed  by  a  summary  of  the  matrix  required  for  the  problem 
solution.  Finally)  a  matrix  profile  map  is  supplied  to  assist 
the  user  in  assuring  that  the  problem  was  set  up  correctly. 

3.4.9  Stress-Strain  Curves  (Figure  3.4.1i) 

Nonlinear  analysis  only,  requires  the  input  of 
stress-strain  data  curves  to  be  utilized  in  performing  the 
solution.  These  curves  are  referenced  by  the  materials 
definition  data  for  each  element  and  are  utilized  by  MAGNA 
for  material  nonlinearities. 

3.4.10  Loads  Data  (Figure  3.4. lj) 

Part  nine  repeats  the  loads  input  information 
that  will  be  utilized  in  the  problem  solution.  This  completes 
the  data  input  process.  The  next  four  sections  deal  with  the 
output  generated  by  the  analytical  solution. 

3.4.11  Node  Displacement  (Figure  3.4.1k) 

MAGNA  will  generate  a  listing  of  the  node 
displacement  and  element  stress-strain  data  for  each  load 
increment  performed  in  the  problem  solution.  A  linear  static 
analysis  will  have  only  one  increment  where  as  a  nonlinear 
static  analysis  may  have  any  practical  number  of  increments 
in  the  solution.  The  node  displacement  list  generated  for 
each  increment  gives  the  total  displacement  of  the  node  from 
the  original  position.  A  zero  displacement  followed  by  an 
asterisk  indicates  a  node  constrained  in  that  direction. 

3.4.12  Element  Stresses  (Figure  3.4.11) 

All  the  elements  of  each  element  type  will  be 
listed  with  the  stresses  and  strains  accumulated  up  to  the 
increment  listed.  A  linear  static  analysis  will  only  have  one 
increment  and  therefore  only  one  solution  listing  for  each 
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element.  A  nonlinear  static  analysis  will  have  any  reasonable 
number  of  solution  steps  and  a  listing  of  stresses  and  strains 
will  be  generated  for  each  solution  increment. 

3.4.13  Solution  Time  and  I/O  Summaries  (Figure  3.4.1m) 

Part  12  is  a  listing  of  solution  times  for  various 
parts  of  the  program  execution  and  of  the  input  and  output 
operations  for  program  execution.  The  user  will  generally  not 
be  too  concerned  about  these. 

3.4.14  Day  File  (Figure  3. 4. In) 

This  last  part  is  a  listing  of  the  CDC6600 
computer  job  control  cards  and  execution  diagnostics  for  the 
execution  of  MAGNA.  The  important  items  to  note  here  are  that 
the  analysis  ran  successfully  and  that  a  new  WINGMPOST  file  was 
successfully  cataloged.  These  conditions  are  bracketed  in 
Figure  3. 4.  In. 

This  is  all  the  information  that  is  present  on  the  output 
tab.  The  engineer  should  find  that  all  this  data  is  relevant 
to  the  structural  problem  solution.  This  output  listing  is 
designed  to  be  used  in  conjunction  with  plots  generated  by 
CONTOUR  and  PLOTBOB  which  help  to  describe  the  alternate  load 
paths  and  stress/strain  distributions  as  well  as  nodal 
deformation  under  the  specified  loads.  Any  questions  regarding 
the  output  tabs  or  MAGNA  should  be  directed  to  the  individual 
listed  on  the  program  header  as  illustrated  in  Figure  3.4.1a. 
Figure  3.4.2  is  a  listing  of  a  linear  static  analysis  output 
for  a  T-38  wing  model. 
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Figure  3.4.1c.  (continued) 
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Figure  3.4.1c.  (continued) 
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Figure  3.4.1c.  (continued) 


Figure  3. 4. Id.  Input  Data  File 
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Figure  3. 4. Id.  (continued) . 
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Figure  3. 4. If.  (continued) 
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Figure  3.4.2.  (continued). 
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Figure  3.4.2.  (continued). 
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Figure  3.4.2.  (continued). 
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Figure  3.4.2.  (continued). 


2  ■ 
Ui 


a. 

o 


4 


> 

Ui 

o 

.1 

a 


-J 

*  K 

2  C/» 

<  w 

• 

oc 

<  *- 

Ui  l/J 

2  tt 

«/> 

M  H 

ui 

-J  U. 

-J 

CO 

V  I 

*- 

o 

2  . 

u. 

M  III 

*  in 

a 

O' 

«« 

> 

.«!  ; 

«*  ■ 

a 

o 


N  3  c  a  j  L/  o  i 


)ooooeDao3nvOirN^eH'}aa«ioKQ*o{9|N^c«iiKoi<)\p9lNUvc*4-}K9(*)tC 

«4  w4  w4  tsi  m 


iood0aoooaNiAcwJNDK)Deoat0O>Mn«w«KDii)iA9lNUN0W4KonlDff>MU\ 
y-i  Wi  *4  r>J  i\J  NNntOn^^<fVMM^<rtOtOiOKKK*«CCO>(r 


«49  OAOfSOO 


Of3OO33ODO*iJkO»0v00'MO3Ooli'«^'t^OMvDirNir»HJNcWvDff'fyU'  «  y*  ** 
«4*4*4«4CSJ  NNnnn«^«J^tfSU\«(Di0NKKS««C^O< 


oeeor>eO&OeeeCSC0OKKN.KK^NKoOD9KK.KKKKKN.K.KK^I 


.KKKKNK 


o 

o 

!S 

* 


- - *‘  w  >  ~  *  "  ‘  ‘  *  - - 


HH«4HAiNNIM(Vi(MNnnn^<tt^lMMA(OO)l0NKUKtt«e(rir 


-4(Vm^U>«flK4D(T3^Nm^ir<£Ke(ri9^NnJU\v0K0ir<9'HfMn^ir«OSe(T»9«4IMnt^«OK«(rD*4(M 
«4«4«4«4«4«4*4«4*4*4CM*J<M<NJ<StNN  <4  t'4tUtO*>tO’0’OIOt«>»OV*>IO'*  .*.»'?**'*  *  Jf  Jt  *  V\  tf>  IT> 


3.00 


Figure  3.4.2.  (continued) 
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Figure  3.4.2.  (continued). 
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Figure  3.4.2.  (continued). 
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Figure  3.4.2.  (continued). 
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Figure  3.4.2.  (continued). 
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Figure  3.4.2.  (continued). 
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3.5  MAGNA  OUTPUT  WINGMPOST  POSTPROCESSOR  FILE 

In  addition  to  the  printed  output  tab  of  the  problem 
solution  MAGNA  also  produces  a  special  file  for  postprocessing 
of  the  data  by  the  programs  CONTOUR  and  PLOTBOB.  This  file  is 
listed  in  Figure  3.5.1.  The  file  contains  information 
essentially  the  same  as  on  the  printed  tab  but  in  a  condensed 
format  to  save  space  on  auxiliary  (disk)  storage.  If  the  user 
is  executing  very  large  nonlinear  problems  with  a  large  number 
of  increments  the  WINGMPOST  file  should  be  written  to  a  tape 
for  storage  rather  than  having  it  remain  on  disk.  This  will 
alleviate  potential  problems  with  losing  the  data  file  as  a 
result  of  the  ASD  permanent  file  management  system.  All  important 
files  should  be  placed  on  tape  at  any  rate,  for  insurance. 

The  user  is  cautioned  against  making  any  modifications 
to  this  output  fiJe  as  it  will  almost  certainly  result  in  errors 
while  trying  to  be  utilized  by  CONTOUR  or  PLOTBOB.  It  is 
important  for  the  user  to  be  certain  that  there  is  sufficient 
space  available  for  the  WINGMPOST  file  and  that  a  free  cycle 
number  is  available  for  the  WINGMPOST  permanent  file  name 
under  the  problem  number  to  be  used.  If  there  is  no  available 
cycle  position, the  day  file  (Section  3.4.14)  will  give  a  message 
to  that  effect.  If  there  is  insufficient  space  available  for 
your  office  the  file  will  be  automatically  purged. 
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Figure  3-5.1.  (continued). 
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Figure  3.5.1.  (continued) 
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Fiqure  3.5.1.  (continued) 
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Fiaure  3.5.1.  (continued). 


Figure  3.5.1.  (concluded) . 


SECTION  4 

POSTPROCESSING  PROGRAMS  PLOTBOB  AND  CONTOUR 

4 . 1  INTRODUCTION 

Finite  element  analysis  tends  to  utilize  a  large  portion 
of  computer  resources  to  accomplish  the  task  of  determining 
whether  or  not  a  structure  is  safe  and  reliable.  The  output 
from  such  a  program  can  be  voluminous  and  therefore  very 
difficult  to  interpret  or  comprehend.  Postprocessing  or 
post-analysis  manipulation  of  the  output  data  from  MAGNA  is 
the  primary  function  of  PLOTBOB  and  CONTOUR  plotting  programs. 
The  programs  address  themselves  to  the  requirements  of 
efficiently  and  understandably  presenting  the  MAGNA  analysis 
results  to  the  user  to  provide  a  more  comprehensible 
understanding  of  what  the  results  actually  indicate  about  the 
model  in  question  and  to  directly  tie  the  model  developed  with 
the  results  achieved. 

PLOTBOB  provides  the  user  with  the  capability  to  plot  a 
structure  prior  to  analysis  utilizing  the  MAGNA  input  file 
(load  deck  created  on  TAPE  11  by  WINGEN)  as  well  as  the  ability 
to  selectively  plot  deformed  structures  utilizing  an  -MPOST 
file.  The  preprocessor  output  file  may  be  input  to  PLOTBOB 
without  any  changes  after  it  is  created  by  WINGEN.  The  user 
may  wish  to  view  the  model  (especially  if  shell  elements  were 
selected  for  upper  and  lower  skins)  with  the  more  sophisticated 
capabilities  of  PLOTBOB,  as  compared  with  the  plotting 
capabilities  of  WINGEN,  to  better  ensure  the  model  was  properly 
defined  prior  to  MAGNA  analysis.  The  user  should  be  aware  that 
no  deformed  plotting  may  occur  with  a  MAGNA  input  file  because 
there  is  no  data  indicating  what  a  deformed  model  would  be  until 
after  the  analysis  run.  PLOTBOB  can  also  be  utilized  for  post¬ 
analysis  plotting  of  the  model.  This  is  limited  to  plotting 
deformed  structures  (contour  plotting  is  limited  to  CONTOUR) 
but  all  the  capabilities  for  selective  viewing  are  available 
for  use  with  the  deformed  structure  plots. 


CONTOUR  is  a  second  general-purpose  plotting  program  for 
use  with  MAGNA  analysis  output  data.  Unlike  PLOTBOB  it  cannot 
plot  any  MAGNA  input  files  but  it  has  the  capability  to  plot 
deformed  geometry  structures  and  superimpose  on  structures 
contour  lines  representing  different  stresses,  strains  or 
displacements.  This  can  be  very  informative  in  assisting  the 
user  with  determining  where  and  when  structural  failure  might 
occur.  With  capability  for  significant  labeling  of  the  model, 
contour  lines  and  the  entire  plot,  CONTOUR  provides  the  user 
with  excellent  documentation  of  output  results.  Figures  4.1.1 
and  4.1.2  illustrate  sample  plots  from  both  PLOTBOB  and  CONTOUR 
to  allow  the  user  a  comparison  of  their  capabilities.  The 
remainder  of  this  chapter  deals  with  how  to  use  these  two 
programs.  Reference  4  provides  some  additional  information 
on  the  postprocessing  programs.  The  program  names  in  Reference 
4  are  different,  however,  PLOTBOB  has  been  renamed  to  GPLOT  and 
CONTOUR  has  been  renamed  to  CPLOT. 


Figure  4.1.1.  Typical  plot  from  PLOTBOB  plotting  program.  This  illustrates  the  lower 
wing  skin  of  a  Replica  Wing  Test  Specimen  with  the  elements  shrunken  an< 
damage  to  the  skin  element  in  chord  bay  1  span  bay  3. 


4.2  PLOTBOB  -  INTRODUCTION 


PLOTBOB  is  a  versatile,  interactive  finite  element  model 
(FEM)  plotting  program.  The  program  can  plot  a  model  both 
before  (preprocessor  output  file  TAPEll)  and  after  (MAGNA 
generated  WINGMPOST  file)  MAGNA  program  execution.  Any  portion 
of  the  model  can  be  viewed  by  specifying  individual  elements  for 
plotting;  zooming  on  the  structure  utilizing  either  of  two 
methods  available;  or  clipping  away  part  of  the  structure  to 
view  inside.  An  exploded  view  option  is  also  available.  The 
model  can  be  rotated,  translated  and  reflected  about  axes  as 
necessary  for  optional  viewing.  It  can  be  viewed  to  scale  with 
an  axonometric  or  three  point  perspective  projection.  The 
user  can  label  the  elements  and  nodes  by  element,  type,  surface 
location  and  edge  location. 

4.2.1  Procedure  for  Executing 

PLOTBOB  may  be  executed  using  the  methods  outlined 
below  for  the  CDC6600  computer: 

LOGIN, . . . 

ATTACH , F , PLOTTINGPROCEDURES , ID=XXXX , SN=XXXX 

ATTACH, TAPE5 ,pfn  (desired  data  file) 

BEGIN, PLTBOB,F  ( ,H  for  HP  7221  plotter) 

(blank  for  Tektronix  terminal) 

PLOTTING  PROCEDURES  is  a  job  control  procedure 
which  performs  the  process  of  attaching  and  defining  the 
libraries  and  attaching  and  executing  the  program  PLOTBOB.  The 
parameters  on  the  BEGIN  command  above  specify  whether  the  user 
is  utilizing  a  Tektronix  graphics  terminal  or  a  Hewlett-Packard 
7221  bed  plotter.  No  additional  information  is  required  if 
using  a  Tektronix  terminal  only:  BEGIN, PLTBOB, F .  An  additional 
character  should  be  appended  to  the  command  if  one  is  using  the 
HP  plotter  as  follows:  BEGIN, PLTBOB, F,H. 

PLOTBOB  will  always  look  for  a  local  file  with  the 
name  'TAPES'  which  contains  the  input  data  defining  the  fem 
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model  to  be  displayed.  This  data  may  be  either  a  preprocessor 
output  file  (TAPE11 )  containing  all  the  data  necessary  for  a 
FEM  analysis  (type  1)  or  it  may  be  a  pem  analysis  output 
file  (-MPOST  file)  containing  all  the  analysis  data  along  with 
the  model  definition  data  (type  2).  The  program  will  ask  the 
user  which  file  type  he  is  using:  1  or  2 .  The  preprocessor 
output  file  is  detailed  in  section  3.2.  The  FEM  analysis  output 
file  is  explained  in  Section  3.3. 

4.2.2  Command  Structure 

Once  PLOTBOB  has  been  initiated  the  program  will 
prompt  the  user  with  several  questions: 

ENTER  THE  CHARACTERS  PER  SECOND . : 

TYPE  OF  PLOTTER  (T,H) . : 

ENTER  DATA  FILE  TYPES: 

1  =  PREPROCESSOR  DATA  FILE 

2  =  MPOST  FILE 

ENTER  DATA  FILE  TYPE  (1,2) . . 

These  questions  control  program  functioning.  The 
characters  transmitted  per  second  (cps)  are  a  function  of  the 
baud  rate:  on  a  CDC  1/10  *baud  rate  yields  the  cps.  The  type  of 
plotter  can  be  either  a  Tektronix  4010  series  (T) ,  a  Tektronix 
4010  series  emulator  (T) ,  a  Hewlett-Packard  7221  (H)  or  any 
terminal  or  plotter  operating  under  Tektronix  PLOT10  Terminal 
Control  system  (TCS)  routines  (T)  or  Hewlett-Packard  HPPLOT21 
graphics  library  routines  (H) .  The  PLOTBOB  program  will  handle 
two  types  of  data  files  as  discussed  above  for  TAPE5,  a 
preprocessor  output  file  or  a  MAGNA  generated  WINGMPOST. 

Once  satisfactory  responses  have  been  received 
by  the  program  it  will  prompt  the  user  with  an  asterisk  (*). 

The  user  may  then  enter  any  valid  command  (for  a  list  of 
commands  type  'HELP')  by  typing  the  first  four  letters  of  the 
command  and  pressing  the  <CR>  (carriage  return) .  Any  number 
of  commands  may  be  entered  and  reentered  until  the  time  limit 


4.6 


of  the  program  has  expired  ,  the  user  issues  the  'STOP*  command 
or  there  is  an  execution  error.  Section  4.2.3  discusses  the 
commands  available. 

4.2.3  Command  Summary 

This  section  lists  each  command  available  to  the 
PLOTBOB  program  user  in  alphabetical  order  with  a  brief 
description  of  each.  Only  the  first  four  letters  of  each 
command  are  required.  Ample  illustrations  have  been  provided 
to  allow  the  user  a  more  complete  understanding  of  each  command. 
Additional  information  may  be  provided  in  Section  4.3.4  for 
the  CONTOUR  program  commands  which  perform  similar  functions. 
Figure  4.2.1  lists  the  commands  covered  in  this  section. 

4. 2. 3.1  AXES 

This  command  allows  the  user  to  plot  a 
set  of  orthogonal  axes  labeled  "x",  "y" ,  and  "z" .  Each 
structure  (finite  element  model)  is  defined  using  a  3-D 
coordinate  system.  The  program  accepts  input  for  Cartesian, 
cylindrical,  and  spherical  coordinates  though  it  can,  by  program 
modification,  accept  other  user  coordinate  systems.  Those 
coordinates  which  are  linear,  rather  than  angular,  in  nature 
can  be  expressed  in  any  units  of  length  the  user  desires.  In 
general,  these  units  are  the  same  as  those  used  in  taking 
measurements  of  the  actual  structure.  All  locating  of  points 
on  the  structure  with  the  user  coordinate  system  and  unit  of 
length  is  done  from  an  arbitrary  (0,0,0)  point  called  the 
origin. 

The  axes  extend  from  the  origin  in  the 
positive  x-,  y-,  and  z-directions.  The  length  of  each  axis  is 
10  percent  greater  than  the  maximum  distance  of  the  structure 
from  the  origin  in  the  same  direction.  For  example,  if  the 
maximum  distance  of  the  structure  from  the  origin  in  the 
x-direction  is  10  mm  (i.e.,  the  maximum  x-coordinate  is  10) 
then  the  length  of  the  x-axis  if  11  mm.  This  is  true  regardless 
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Figure  4.2.1.  PLOTBOB  commands  summary.  This  list  may  be  obtained  by  issuing  the 
' HELP '  command . 


of  whether  or  not  the  structure  coincides  with  the  origin.  The 
x-axis  length  is  based  only  on  the  maximum  x-coordinate .  No 
axes  are  plotted  when  the  default  value  is  specified. 

4. 2. 3. 2  CLIP  (Figures  4.2.2  and  4.2.3) 

CLIP  gives  the  user  an  opportunity  to 
clip  off  a  portion  of  the  structure  to  be  plotted  in  order  to 
look  inside.  The  user  specifies  the  location  of  the  clip  plane 
by  entering  a  clip  factor,  which  is  defined  as 

CLP FAC  =  C/P 

where  C  is  the  distance  from  the  eye  to  the  clip  plane  and  P 
is  the  distance  from  the  eye  to  the  site  position.  The  site 
position  is  the  centroid  of  the  structure.  Any  portion  of  the 
structure  between  the  eye  and  the  clip  plane  will  not  be  mapped 
onto  the  virtual  plane  (see  SCALE) ,  and  thus  will  not  appear 
at  the  plotter.  The  default  value  for  CLPFAC  is  0.01. 


Figure  4.2.3.  Illustration  of  CLIP  command  definitions 

4.2. 3.3  CUBE  (Figure  4.2.4) 

CUBE  changes  the  viewing  box  size  in 
order  to  give  close-ups  of  particular  parts  of  a  structure. 


Figure  4.2.4.  Illustration  of  clamped  beam  elements  viewed  after 
utilization  of  the  CUBE  command.  This  is  a  zoom  on 
the  structure  illustrated  in  Figure  4.2.14. 


r 


For  default,  the  vertices  of  the  imaginary  viewing  box  coincide 
with  the  maximum  and  minimum  points  of  the  structure  in  each 
direction.  Using  CUBE,  the  user  can  change  these  maximum 
and  minimum  values  to  contain  only  that  portion  of  the  structure 
he  wishes  to  observe.  Once  these  values  are  chosen,  the  user 
can  then  view  this  portion  from  any  position  he  desires  without 
having  to  re-enter  the  CUBE  command  each  time  to  obtain  a  new 
virtual  plane  and  image  (which  the  user  must  do  with  the  ZOOM 
command) .  To  use  this  command  effectively,  the  user  should  be 
familiar  with  the  dimensions  of  the  structure  being  plotted  in 
order  to  select  suitable  maximum  and  minimum  values  for  the 
new  viewing  box. 

4 . 2 . 3 . 4  DEFAULT 

DEFAULT  sets  the  plotting  parameters 
controlled  by  user  commands  to  correspond  with  the  default 
values  given  in  the  descriptions  of  each  command. 

4. 2. 3. 5  DEFORM  (Figure  4.2.5) 

DEFORM  provides  the  user  with  the 
capability  of  one  of  three  plotting  options:  (1)  undeformed 
structure  plot;  (2)  deformed  structure  plor,  or  (3)  both 
deformed  and  undeformed  structure  plots.  The  undeformed  option 
gives  a  plot  of  the  structure  as  it  would  appear  prior  to 
loading  being  applied.  If  the  user  is  executing  this  program 
with  file  type  1  (input  load  deck  data  file)  this  is  the  only 
option  available.  A  postprocessor  file  (-MPOST  file)  will 
contain  information  allowing  for  deformed  structure  plots. 

A  deformed  structure  plot  represents  only 
the  deformed  structure  as  a  result  of  loads  applied  during 
analysis.  Only  -MPOST  post  analysis  files  may  have  deformed 
structure  plots. 

The  last  option  available  is  for  both 
deformed  and  undeformed  structure  plots.  The  undeformed 
structure  is  drawn  with  a  dashed  line  and  the  deformed  structure 
is  drawn  in  solid  lines,  superimposed  on  the  undeformed  structure. 
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Figure  4.2.5.  Deformed  (solid  line) 
unde formed  (dash  line) 
a  loaded  beam  structur 
illustrates  the  DEFO  c 


Once  the  plot  option  has  been  selected 
the  program  will  display  the  maximum  and  minimum  coordinate 
values  for  the  undeformed  case.  If  the  user  has  selected  a 
deformed  structure  plot  the  program  will  then  request  a  file 
increment  number.  A  file  increment  number  corresponds  to  one 
of  the  load  intervals  requested  during  the  analysis  run.  When 
the  load  deck  was  created  (Section  2.9)  the  user  input  a  load 
increment  value  and  the  number  of  steps  or  increments  over  which 
the  load  was  to  be  applied.  The  increment  requested  here 
reflects  which  step  of  the  loading  process  the  user  wishes  to 
view;  1  is  low  (lightest  load  or  only  load  in  most  linear 
analyses)  and  N  is  the  highest  load  increment  available  on  the 
file  (generally  20  for  nonlinear  static  analyses) .  Once  the 
user  selects  a  load  increment  the  program  will  return  the 
maximum  and  minimum  displacements  of  the  deformed  structure. 

The  user  can  compare  this  max/min  value  with  the  undeformed 
max/min  value  to  determine  an  appropriate  scale  factor  to 
multiply  the  deformed  structure  coordinates  by  in  order  to 
better  resolve  the  two  structures  as  well  as  to  expand  or 
exaggerate  the  model  deformities.  The  program  will  request 
the  scale  factor  as  the  last  parameter  for  this  command.  The 
default  for  DEFORM  is  undeformed  structure  plot  only. 

4 . 2 . 3 . 6  DRAW 

The  DRAW  command  plots  a  structure 
according  to  the  values  stored  in  the  plotting  parameters. 

4.2. 3. 7  ELEMENTS  (Figure  4.2.6) 

The  user,  employing  the  ELEMENTS 
command,  can  choose  which  elements  he  wishes  to  plot.  After 
entering  the  number  of  different  element  types  to  be  plotted, 
the  program  requests  the  first  element  type  from  which  to  select 
elements  for  plotting.  The  available  element  types  are 


his  figure  illustrates  a  new  vertical 
VERTICAL  command  see  Figure  4.2.5.  fo 


TYPE 


DESCRIPTION 


Variable  number  of  nodes 

1  solid  (up  to  27  nodes) 

2  Eight  node  brick 

3*  Four  node  plate  element 

4*  Two  node  bar  element 

5*  Eight  node  shell  element 

6  Twenty  node  solid  element 

7  Variable  number  of  nodes 
solid  (up  to  20  nodes) 

8  Sixteen  node  solid  element 
^denotes  element  types  used  for  wing  modeling 

Of  course,  not  every  structure  will  contain  all  of  these,  so 
the  user  must  be  familiar  with  the  data  he  wishes  to  plot.  After 
this  step,  the  user  enters  the  elements  of  his  choice  using  one 
of  three  methods. 

Method  1  -  Select  random  elements  by 
entering  the  number  of  elements  to  be  plotted  followed  by  the 
element  numbers  themselves. 

Example  -  entering  4,  1,  8,  10,  23  would 
plot  four  elements,  numbers  1,  8,  10,  and  23. 

Method  2  -  Select  a  range  of  elements 
by  entering  the  first  and  last  element  in  the  range,  followed 
by  an  increment  number. 

Example  -  entering  1,  13,  3  would  plot 
elements  1,  4,  7,  10,  and  13. 

Method  3  -  Plot  all  of  the  elements  for 
the  present  type.  Elements  are  numbered  consecutively  through 
each  element  type.  Thus,  if  a  structure  consists  of  15  type  1 
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elements  and  five  type  2  elements,  the  first  set  of  elements 
would  be  numbered  1  through  15,  and  the  second,  16  through  20. 
This  procedure  is  then  repeated  for  the  next  element  type  and 
so  on,  until  all  the  elements  the  user  desires  to  plot  have 
been  entered.  The  default  is  that  all  elements  for  all  types 
present  are  plotted. 

4. 2. 3. 8  EYE  (Figure  4.2.7) 

A  structure  may  be  viewed  from  any 
position  in  space.  This  position  is  defined  by  a  3-D  Cartesian 
coordinate  measured  relative  to  the  structure  origin  (see  AXIS). 
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Figure  4.2.7.  Default  Viewing  Position  for  EYE  Command. 

The  default  value  is  (100,100,100).  EYE  allows  the  user  to  move 
to  a  different  viewing  location.  There  is  no  r»«*-»-iction  on  this 
location;  however,  if  it  is  located  too  close  to  or  within  the 
structure,  portions  of  the  resulting  plot  will  be  clipped,  i.e., 
part  of  the  structure's  image  will  be  cut  off  the  screen. 
Subroutine  CLIP  contains  a  full  explanation  of  the  clipping 
process. 

4.2. 3. 9  HELP  (Figure  4.2.8) 

HELP  gives  a  list  of  all  available 
commands  with  a  brief  description  of  each. 
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Figure  4.2.8.  Illustration  of  command  listing  generated  by  HELP  command. 
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4.2.3.10  LABEL  (Figures  4.2.9  and  4.2.10) 

LABEL  permits  the  user  to  label  the 
elements  and  nodes  of  a  structure.  If  the  user  responds 
affirmatively  to  the  question  concerning  the  labeling  of  the 
elements,  all  elements  which  are  plotted  will  be  labeled. 

Labeling  for  each  element  occurs  at  its  centroid. 

To  label  the  nodes,  the  user  first 
specifies  the  number  of  element  types  which  will  receive  node 
labeling.  He  then  enters  an  element  type  followed  by  a 
number  corresponding  to  the  surface  he  wishes  to  label.  The 
choices  available  are  0  through  6  inclusive,  where  0  denotes 
the  labeling  of  all  the  nodes.  Choices  1  through  6  for  solid 
and  plane  elements  are  shown  in  Figures  4.2.9  and  4.2.10. 

Note  that  the  surface  and  edges  are 
defined  by  the  local  node  numbers  encompassing  them,  rather 
than  by  their  orientation  in  space.  Thus,  for  example,  surface 
5  is  always  the  surface  containing  local  nodes  1,  2,  5,  and  6. 

The  process  is  then  repeated  for  each  element  type.  Default  is 
no  labeling  of  elements  and  nodes. 

4.2.3.11  NEW 

With  this  command  the  user  can  request 
that  a  different  structure  be  plotted.  There  must  be  more  than 
one  structure  present  on  the  data  file  (TAPE 5) .  A  different 
structure  is  specified  using  the  number  corresponding  to  its 
location  on  TAPE5,  i.e.,  to  specify  the  third  structure,  the 
user  enters  a  3.  Each  set  of  data  has  a  header  card  with  the 
string  "COOR"  in  the  first  four  columns.  This  line  immediately 
precedes  all  the  information  concerning  the  nodal  coordinates 
and  element  connectivities  for  that  structure.  Generally, 
multiple  linear  analyses  may  be  executed  on  one  structure 
generating  one  structure  with  several  increments,  similar  to  a 
nonlinear  analysis  with  several  load  steps.  Several  such  output 
files  may  be  placed  on  a  single  file  following  the  analyses 
which  would  result  in  multiple  structures  on  the  file. 
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Figure  4.2.10.  Example  of  clamped  beam  structure  with  nodes  labeled  on  surface 
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4.2.3.12  PROJECTION  (Figures  4.2.11,  4.2.12  and  4.2.13) 

PROJECTION  gives  the  user  his  choice 
of  pictorial  projection.  The  two  choices  available  are 
orthogonal  (axonometric) ,  which  is  requested  by  the  letter  0, 
and  perspective  (three-point  perspective) ,  which  is  requested 
by  the  letter  P.  The  orthogonal  projection  is  formed  by  using 
projection  lines  which  are  perpendicular  to  the  virtual  plane 
described  in  SCALE.  The  result  is  an  image  in  which  lines 


Structure 


Projection  Lines 


Projection^ 
(Virtual)  Plane 


Final  Image 


Orthogonal 


Perspective 


Figure  4.2.11.  Illustration  of  the  derivation  of  the 
2  projection  types,  orthogonal  and 
perspective. 


parallel  to  the  vertical  axis  remain  vertical,  and  none  of  the 
other  lines  recede  away  from  the  viewer  as  in  a  perspective 
view.  The  perspective  view  causes  all  lines  to  recede,  as 
shown  in  Figure  4.2.12.  The  default  projection  is  perspective. 


4.2.3.13  REFLECT  (Figure  4.2.14) 

REFLECT  allows  the  user  to  reflect  a  plot 
about  the  yz-,  xz-,  or  xy-plane.  Both  the  original  image  and 
its  reflection  are  plotted.  Default  is  no  reflection. 

4.2.3.14  ROTATE  (Figure  4.2.14) 

ROTATE  gives  the  user  the  chance  to 
rotate  the  structure  about  the  x-,  y-,  and  z-axes.  The  user 
specifies  a  rotation  angle,  measured  in  degrees,  for  each  axis. 
The  angles  are  positive  if  they  are  counter-clockwise  when 
looking  down  the  positive  side  of  the  axis  toward  the  origin. 
Default  is  no  rotation. 

4.2.3.15  SCALE  (Figures  4.2.13  and  4.2.15) 

PLOTBOB  will  map  the  3-D  structure  onto 
a  2-D  "virtual  plane"  as  shown  in  Figure  4.2J.3.  The  dimensions 
A  and  B,  which  have  the  same  units  of  length  as  the  structure, 
define  a  plane  just  large  enough  to  contain  the  image.  This 


Figure  4.2.13.  Illustration  of  the  relationship  between  the 
virtual  plane,  the  structure  and  the  viewing 
position. 


plane  and  the  image  on  it  must  then  be  stretched  and  shrunk  as 
necessary  to  fill  the  output  device's  bed  or  screen  where  it  is 
to  be  viewed.  This  process  naturally  distorts  the  proportions 
of  the  final  structure  image,  making  it  out  of  scale. 
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Figure  4.2.14.  Plot  of  deformed  and  undeformed  clamped  beam  structure 
illustrating  TRANSLATE,  ROTATE  and  REFLECT  command.  T 
original  structure  is  illustrated  in  Figure  4.2.10. 


The  purpose  of  the  SCALE  command, 


therefore,  is  to  give  the  user  an  opportunity  to  view  the  final 
image  with  proper  proportions.  This  is  accomplished  as  follows: 

For  a  Tektronix  4014  screen, 

IF  A  >  1 . 33*B  then  B  =  0.75*A 
IF  A  <  1 . 33*B  then  A  =  1.33*B 


For  an  HP  7221  bed  plotter, 

IF  A  >  1 . 286*B  then  B  =  0.778*A 
IF  A  <  1.286*B  then  A  =  1.286*B 


This  command  allows  the  program  to  set  the  size  of  the  virtual 
plane  to  be  the  same  proportion  as  the  plotter  screen  without 
changing  the  size  of  the  image.  The  default  for  this  command 
plots  the  structure  with  the  proper  scale  for  either  PLOTIO  or 
HPPLOT21  library  users. 

4.2.3.16  SHRINK  {Figure  4.2.16) 

Three  basic  shapes  of  elements  are 
plotted  by  the  program  -  beams  (line  segments) ,  plates  (planes) , 
and  six-sided  solids.  When  a  beam  coincides  with  a  plate  or 
solid  edge,  or  a  plate  coincides  with  a  solid  surface  in 
multiple-element-type  plots,  the  user  cannot  distinguish  the 
various  elements  used  and  their  locations.  SHRINK  gives  the 
user  an  exploded  view  so  that  all  elements  may  be  seen  clearly. 

SHRINK  achieves  this  by  contracting  the 
elements  around  their  centroids,  thus  withdrawing  them  from  all 
adjacent  elements.  The  user  specifies  a  shrink  factor  from  0  to 
1  inclusive,  with  0  being  no  shrinkage  and  1  being  shrinkage 
down  to  the  centroid.  SHRINK  determines  the  centroid  for  each 
element  and  then  determines  the  new  location  of  each  node 
according  to  the  general  formula. 


P 

an 


C 

a 


F  +  (1  -  F  ) 
s  s 


p 

ao 
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where 


P  is  the  new  node  coordinate 
an 

P  is  the  old  node  coordinate 
ao 

C  is  the  centroid  coordinate 

d 

F  is  the  shrink  factor, 
s 

A  shrink  factor  of  0.1  works  well.  Default  is  no  shrinkage. 

4.2.3.17  STOP 


the  program. 


Issuing  a  STOP  terminates  execution  of 


4.2.3.18  SUMMARY  (Figure  4.2.17) 

SUMMARY  lists  the  current  options  chosen 
by  the  user  for  each  command. 


4.2.3.19  TIME 


When  the  user  gives  the  TIME  command  on 
a  CDC6600  computer,  the  computer  prints  the  CPU  time  used  in 
seconds,  since  the  LOGIN  was  issued.  Since  the  total  available 
time  for  a  session  is  typically  limited  to  an  installation-defined 
maximum,  and  since  the  user,  after  some  experience,  will  know 
the  approximate  time  for  a  plot,  he  can  make  a  series  of  plots 
without  exceeding  the  time  limit. 

4.2.3.20  TRANSLATE  (Figure  4.2.14) 

TRANSLATE  translates  a  structure  when 
the  user  specifies  a  translation  vector.  The  units  of  the  vector 
are  the  same  as  for  the  coordinates  of  the  structure.  Thus,  if 
the  user  wishes  to  move  the  structure  two  units  down  the  x-axis 
and  five  units  up  the  z-axis,  the  user  would  enter  2.,0.,5. 

Default  is  no  translation. 


4.2.3.21  VERTICAL  (Figure  4.2.6) 

Using  VERTICAL,  the  user  can  choose 
which  axis  (x,  y,  or  z)  should  appear  vertically  on  the  plot. 
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Figure  4.2.17.  An  example  of  the  SUMMARY  command. 


The  coordinates  are  adjusted  accordingly  when  the  user  selects 
an  axis  other  than  the  default  z-axis  for  the  vertical  axis. 

4.2.3.22  ZOOM  (Figure  4.2.18) 

ZOOM  is  the  second  command  the  user  can 
employ  in  getting  a  blow-up  of  a  portion  of  a  structure.  The 
structure  must  first  be  plotted  with  the  view  desired.  Then 
the  user  specifies  a  virtual  window  which  encloses  that  part 
of  the  plot  on  the  screen  the  user  desires  to  view  more  closely. 

This  is  done  with  a  Tektronix  terminal 
using  the  cursor  (two  fine  crossed  lines  which  appear  on  the 
screen  once  the  ZOOM  command  is  implemented) .  The  user  adjusts 
the  cursor  to  the  location  of  the  lower  left-hand  corner  of  the 
window  he  wishes  to  create  using  the  cursor  controls.  He  then 
enters  any  alphanumeric  character  from  the  terminal  and  depresses 
the  "RETURN"  key  if  the  terminal  is  configured  to  require  a 
carriage  return.  He  repeats  this  procedure  to  define  the  upper 
right-hand  corner  of  the  window  he  wishes  to  create.  The  user 
may  now  plot  the  'windowed'  area  at  full  scale  on  the  screen  or 
plotter . 

To  ZOOM  on  an  HP7221  bed  plotter  the 
user  must  wait  until  the  'ENTER'  light  on  the  plotter  control 
panel  begins  to  flash;  he  then  must  move  the  pen  to  the  lower 
left  position  of  the  ZOOM  area  with  the  pen  controllers  (buttons 
with  arrows  on  them) .  To  center  the  coordinate  of  the  lower  left 
corner  the  user  must  press  the  'ENTER'  button.  This  enters  one 
x,y  coordinate  pair.  The  user  must  then  move  the  pen  to  the 
upper  right  corner  of  the  area  to  be  expanded  with  the  pen 
controllers  and  again  press  the  'ENTER'  buttom.  The  program 
will  take  these  two  (x,y)  coordinate  pairs  and  create  a 
'window'  onto  which  the  zoomed  area  will  be  plotted.  The  user 
may  now  plot  the  structure. 

To  get  a  different  view  of  the  same  area, 
the  user  must  repeat  the  entire  process.  It  should  be  noted 
that  the  final  image  is  out  of  scale  (the  problem  is  similar  to 
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the  one  of  the  virtual  plane  being  mapped  onto  the  screen ,  as 
described  in  SCALE) .  Refer  to  the  CUBE  command  for  an 
alternative  method  of  enlarging  a  specific  area  of  a  model. 

It  is  important  to  make  sure  your  terminal  or  plotter  has 
digitizing  capability  prior  to  executing  this  command. 
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4.3  CONTOUR  -  INTRODUCTION 

CONTOUR  is  an  interactive  computer  program  for  the 
graphical  display  of  data  obtained  through  finite  element  model 
analysis.  CONTOUR  displays  data  either  in  the  form  of  contours, 
relief  maps  or  some  combination  of  the  two.  It  is  structured 
in  such  a  way  as  to  make  the  execution  of  the  program  as  simple 
as  possible.  CONTOUR'S  capabilities  include: 

•  plotting  of  stress,  strain  or  displacements  contours, 

•  labeling  of  contours, 

•  plotting  of  relief  maps  of  stress,  strain  or 
displacement , 

•  undeformed  and  deformed  structure  plotting, 

•  variable  contour  quantity  and  quality, 

•  variable  3-D  views, 

•  exploded  views  for  easier  model  comprehension, 

•  zooming, 

•  selective  element  plotting, 

•  descriptive  plot  labeling. 

CONTOUR  is  used  only  with  MAGNA  analysis  output  files 
(-MPOST  files)  and  proves  to  be  most  effective  for  high  quality 
'finished'  contour  and  displacement  plots  of  the  finite  element 
model.  Further  plotting  capabilities  are  detailed  in  Section 
4.2  -  PLOTBOB  plotting  program. 

4.3.1  Procedure  for  Executing 

CONTOUR  is  executed  by  utilizing  a  command 
procedure  as  illustrated  below: 

LOGIN . 

ATTACH , F , PLOTTINGPROCEDURES , ID=XXXX , SN=XXXX 

ATTACH, TAPE9 9, pfn  (desired  -MPOST  file) 

BEGIN, PLOT,F  (,H  is  added  for  HP  7221  plotter) 

(blank  for  Tektronix  terminal) 
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PLOTTINGPROCEDURES  is  a  job  control  procedure  which  performs 
all  tasks  associated  with  executing  CONTOUR.  The  parameters 
on  the  BEGIN  command  specify  whether  the  user  is  utilizing  a 
graphics  terminal  with  PLOTIO  software  (Tektronix)  or  a  bed 
plotter  with  HPPL0T21  software  (Hewlett-Packard) .  If  a  graphics 
terminal  is  being  used  only  the  command  BEGIN, PLOT, F.  is 
required.  For  use  of  the  program  with  a  bed  plotter  an  additional 
parameter  is  required  as  follows:  BEGIN, PLOT, F,H.  to  utilize 
the  proper  graphics  commands. 

4.3.2  Command  Structure 

The  execution  of  CONTOUR  plotting  program  is 
achieved  in  two  phases:  Part  I  -  Data  Initialization  and 
Part  II  -  Command  Directives.  Data  Initialization  identifies 
the  data  set  to  be  plotted  and  subsets  of  the  data  file  to  be 
utilized.  Command  Directives  control  the  actual  plotting 
including  selection  of  viewing  angles,  element  plotting,  contour 
and  relief  maps,  etc.  Use  of  the  NEST  and  NEWD  commands  will 
return  the  user  to  the  Data  Initialization  phase.  All  other 
commands  will  control  the  data  set  established  by  the  previous 
Data  Initialization  phase. 

4.3.3  Data  Initialization 

This  phase  of  CONTOUR  defines  the  data  set  to  be 
plotted.  The  program  will  prompt  the  user  with  a  series  of 
questions  to  direct  the  input  of  the  correct  data  for  the  plots 
desired.  These  questions  are: 

STRUCTURE  NUMBER? 

NUMBER  OF  ELEMENT  TYPES  TO  BE  AVERAGED? 

ELEMENT  TYPES  TO  BE  USED  IN  COMPUTING 

SMOOTHED  VALUES? 

ELEMENT  TYPE  3  SUBTYPE? 

VALUE  PLOT  OPTION? 

DO  YOU  WISH  TO  LIST  NODAL  VALUES? 
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STRESS/STRAIN  (if  stress  or  strain  chosen  above) 

INCREMENT  NUMBER? 

The  Structure  number  pertains  to  the  order  in  which  the  structure 
to  be  plotted  appears  on  the  -MPOST  file.  Most  cases  will 
involve  only  one  structure  per  nonlinear  analysis  and  possibly 
several  structures  for  a  linear  analysis. 

MAGNA  will  generate  stress,  strain  and  displacement 
values  at  selected  points  within  each  element.  In  order  to 
generate  approximate  values  for  stress  strain  or  displacement 
CONTOUR  extrapolates  these  interior  "integration"  values  to 
the  nodes.  Since  each  node  may  be  utilized  in  several  different 
elements  it  is  best  to  average  all  the  nodal  values  obtained 
from  extrapolating  from  integration  points  to  a  common  node 
from  the  same  element  type.  If  a  node  is  included  in  more 
than  one  element  type,  the  user  should  only  average  one  of  those 
element  types  at  a  time  to  get  a  true  representation  of  the 
stresses,  strains  and  displacements  on  that  node  due  to  the 
various  components.  For  Wing  model  plotting  the  user  will  wish 
to  average  either  element  types  3  or  element  type  5.  (if  shell 
elements  were  used)  but  only  one  at  a  time.  In  structures 
where  several  element  types  are  used  that  do  not  involve 
common  nodes  more  than  one  element  type  may  be  averaged  at  a 
time. 

Once  the  user  determines  how  many  element  types 
he  wishes  to  average  he  is  asked  to  input  the  actual  element  types 
he  wishes  to  average.  The  number  of  ELEMENT  TYPES  TO  BE  USED 
IN  COMPUTING  SMOOTHED  VALUES  must  correspond  to  the  value 
entered  for  the  question  NUMBER  OF  ELEMENT  TYPES  TO  BE  AVERAGED? 

If  element  type  3  is  used  above  the  program  will 
request  a  subtype:  ELEMENT  TYPE 3  SUBTYPE?  There  are  three 
subtypes  available: 

1  =  plane  stress 

2  =  plane  strain 

3  =  shear  panel. 
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The  subtype  requested  must  be  included  in  the  model  data.  The 
model  will  always  have  subtype  3  and  will  utilize  subtype  1  if 
no  shell  elements  are  generated. 

A  VALUE  PLOT  OPTION  is  requested  by  the  program 
to  determine  what  type  of  information  will  be  plotted.  The  value 
plot  option  may  be  one  of  four  choices: 

0  =  stress  or  strain  plot 

1  =  x  displacement  plot 

2  =  y  displacement  plot 

3  =  z  displacement  plot. 

Should  the  user  select  0,  stress  or  strain  plot,  a  question  will 
follow  later  requesting  the  stress  or  strain  code  (see  below) . 

Occasionally,  it  _ s  important  to  know  the  actual 
values  of  the  nodes  as  determined  by  utilizing  the  integration 
points  and  averaging  the  values.  DO  YOU  WISH  TO  LIST  NODAL 
VALUES?  provides  the  user  with  the  option  of  having  a  list  of 
these  values  printed  at  the  terminal  once  the  data  set  has  been 
defined  and  before  any  plotting  commands  may  be  initiated. 

The  STRESS/STRAIN  CODE  is  requested  if  the  user 
responded  that  he  wishes  to  plot  stresses  or  strains  (option  0 
to  question  of  VALUE  PLOT  OPTION) .  The  following  13  codes 
are  value  stress/strain  codes: 


1  =  EXX 

2  =  EYY 

3  =  EZZ 

4  =  EYZ 

5  =  EXZ 

6  =  EXY 

7  =  SXX 

8  =  SYY 

9  =  SZZ 

10  =  SYZ 

11  =  SXZ 

12  =  SXY 

13  =  equivalent  stress 

where  E  =  strain  components 
and  S  =  stress  components 

One  of  the  above  codes  must  be  selected.  The  program  cannot 
determine  whether  a  selected  code  will  generate  the  proper 
stresses  and  strains  for  the  element  type  selected  and  its 
orientation  in  space,  this  is  up  to  the  user  to  determine. 
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The  last  question  to  be  answered  in  this  part  is 
which  increment  of  the  data  set  the  user  wishes  to  plot. 
Increments  are  established  as  either  time  intervals  of  a  time 
dependent  analysis  or  as  a  load  interval  for  multi-step  loading 
analysis.  Each  time  interval  or  load  increment  is  considered 
a  data  increment  and  ranges  from  1  (for  most  linear  static 
analyses)  to  N  where  N  can  be  any  number  of  intervals  the  user 
selects.  Twenty  intervals  are  suggested  for  a  nonlinear  static 
analysis  of  wing  models. 

4.3.4  Command  Summary 

This  second  phase  of  CONTOUR  execution  allows  the 
user  the  use  of  various  plotting  options  to  most  effectively 
view  the  data  set  established  in  the  first  phase  above.  These 
commands  may  also  allow  the  user  to  alter  the  data  set  or 
select  a  new  increment  or  structure  for  plotting.  If  the  user 
desires  to  return  to  phase  one  described  above  he  can  enter  the 
commands  NEWD  or  NEWS  as  desired. 

Commands  are  entered  as  four  letter  keywords  which 
prompts  the  program  to  request  the  pertinent  information  to 
alter  the  plot.  Any  number  of  commands  may  be  executed  and 
any  command  may  be  executed  any  number  of  times  until  the  plots 
are  to  the  user's  satisfaction.  CONTOUR  has  been  designed  in 
such  a  way  that  no  commands  will  be  necessary  unless  the  user 
wishes  to  alter  a  default  parameter.  The  remainder  of  this 
section  lists  each  command  alphabetically  and  explains  its  use. 
Figure  4.3.1  provides  a  list  of  the  commands  discussed.  Should 
the  user  desire  additional  information,  he  may  refer  to  similar 
commands  discussed  in  Section  4.2  for  PLOTBOB. 

4. 3. 4.1  ALEL 

The  user,  employing  the  ALEL  command, 
can  choose  which  elements  he  wishes  to  plot.  After  entering 
the  number  of  different  element  types  to  be  plotted,  the  program 
requests  the  first  element  type  from  which  to  select  elements 
for  plotting.  The  available  element  types  are 


CONTROL  COMMANDS 
xxxxxxx  XXXXXXXX 


COMMAND 
NEST 
NEUD 
CUBE 
TIME 
EXIT 
PLOT 
HELP 
RESE 
SUMM 
UA  IT 


FUNCTION 

PREPARE  NEU  STRUCTURE  DATA 

PREPARE  NEU  DATA  SET  FOR  PLOTTING 

CHANGE  MAX/MIN  UALUES 

PRINT  CURRENT  CPU  TIME  USED 

EXIT  PROGRAM 

PLOT  UITH  PRESENT  DATA 

PRINT  COMMAND  OPTIONS  AND  FUNCTION 

RESET  TO  DEFAUL“  OPTIONS 

PRINT  THE  CURRENT  OPTION  CONDITION 

PUT  PROGRAM  INTO  PAUSE 


OPTIONS 

XttXXXt 


COMMAND 

FUNCTION 

DEFAULT 

LELE 

LABEL  ELEMENTS 

NO 

LAXS 

LABEL  AXIS 

NO 

LABE 

PLOT  TITLE  BLOCK 

NO 

ALEL 

PLOT  ONLY  SOME  ELEMENTS 

NO 

ENTE 

PLOT  ENTIRE  ELEMENT 

NO 

EXPL 

PLOT  EXPLODED  UIEU 

NO 

200M 

EXPAND  AREA  OF  INTEREST 

NO 

POSA 

PLOT  A  POSITIUE  ARROW 

NO 

SUBT 

SUBTITLE  IN  TITLE  BLOCK 

NO 

DEFO 

SELECT  GEOMETRY 

UNDEFORMED 

CORE 

SELECT  CONTOUR  OR  RELIEF 

CONTOUR 

CLIP 

DISTANCE  TO  CLIP  PLANE 

.01 

EYEP 

SELECT  EYE  POSTION 

100, 100, 100 

SITE 

CHANGE  SITE  POSTION 

CENTROID 

UERT 

SELECT  UERT ICAL  AXIS 

Z 

REFL 

SELECT  A  REFLECT  OPTION 

NONE 

PROJ 

SELECT  PROJECTION  TYPE 

PERSPECT. 

SURF 

SELECT  PLOTTING  SURFACE 

6 

CONL 

LABEL  CONTOURS 

NO 

STEP 

CHANGE  STEP  SIZE 

.  05 

NODE 

LABEL  NODES 

NO 

Figure 

4.3.1.  List  of  commands  available  for  CONTOUR. 
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TYPE 


DESCRIPTION 


2 

3* 

4* 

5* 

6 

7 

8 


Variable  number  of  nodes 
solid  (up  to  27  nodes) 

Eight  node  brick 

Four  node  plate  element 

Two  node  bar  element 

Eight  node  shell  element 

Twenty  node  solid  element 

Variable  number  of  nodes 
solid  (up  to  20  nodes) 

Sixteen  node  solid  element 


♦denotes  element  types  used  for  wing  modeling. 


Of  course,  not  every  structure  will  contain  all  of  these,  so 
the  user  must  be  familiar  with  the  data  he  wishes  to  plot.  After 
this  step,  the  user  enters  the  elements  of  his  choice  using  one 
of  three  methods. 

Method  1  -  Select  random  elements  by 
entering  the  number  of  elements  to  be  plotted  followed  by  the 
element  numbers  themselves. 

Example  -  entering  4,  1,  8,  10,  23  would 
plot  four  elements,  numbers  1,  8,  10,  and  23. 

Method  2  -  Select  a  range  of  elements 
by  entering  the  first  and  last  element  in  the  range,  followed 
by  an  increment  number. 

Example  -  entering  1,  13,  3  would  plot 
elements  1,  4,  7,  10,  and  13. 

Method  3  -  Plot  all  of  the  elements  for 
the  present  type.  Elements  are  numbered  consecutively  through 
each  element  type.  Thus,  if  a  structure  consists  of  fifteen 
type  1  elements  and  five  type  2  elements,  the  first  set  of 


elements  would  be  numbered  1  through  15,  and  the  second,  16 
through  20.  This  procedure  is  then  repeated  for  the  next 
element  type  and  so  on,  until  all  the  elements  the  user  desires 
to  plot  have  been  entered.  The  default  is  that  all  elements  for 
all  types  present  are  plotted. 

4. 3.4.2  CLIP  (Figure  4.3.2) 

The  CLIP  command  allows  the  user  to 
designate  the  distance  from  the  eye  to  the  clip  plane. 


Figure  4.3.2.  Illustration  of  how  to  define  a  clip  plane  for 
the  CLIP  command. 

The  distance  is  defined  by  the  ratio  c/s,  where  c  is  the 
distance  from  the  eye  to  the  clip  plane  and  s  is  the  distance 
from  the  eye  to  the  site.  Any  part  of  the  structure  that  is 
closer  to  the  eve  than  the  clip  plane  will  be  clipped.  This  is 
espectially  useful  for  viewing  inside  or  very  close  to  the 
structure.  The  default  value  is  .01.  The  distance  should  never 
be  greater  than  one,  and  depending  on  the  eye  position  should 
be  a  fairly  small  number. 
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4. 3. 4. 3  CONL  (Figures  4.3.3  and  4.3.9) 


The  CONL  command  lables  the  contours 
so  the  user  can  easily  identify  values  of  the  contours.  When 
both  CONL  and  LABE  options  are  active,  a  table  of  contour  values 
is  printed  in  the  lower  right  hand  side  of  the  plot. 

4. 3. 4. 4  CORE  (Figures  4.3.3  and  4.3.4) 

The  CORE  command  allows  the  user  the 
choice  of  plotting  either  relief  maps,  contour  plots  or  both. 

A  relief  map  consists  of  a  surface 
constructed  above  and  below  the  plotting  surface.  At  each 
poini:  the  surface  is  scaled  proportional  to  the  value  of  the 
component  being  plotted  at  that  point.  The  relief  map  is  scaled 
so  the  maximum  displacement  from  the  surface  is  the  relief  scale. 

A  contour  plot  consists  of  contours 
of  equal  value.  The  value  is  that  of  the  component  being 
plotted.  These  contours  are  drawn  in  incremental  steps  (see 
STEP  command) ,  the  values  for  which  are  designated  by  the  user. 
The  contours  can  also  be  labeled  (see  CONL  command) . 

If  the  user  selects  the  contour  option, 
the  user  will  be  given  the  maximum  and  the  minimum  value  of  the 
component  being  plotted.  The  user  will  then  be  asked  to  enter 
minimum,  maximum  and  increment  values  for  the  contours.  If 
the  increment  value  is  zero,  the  input  minimum  and  maximum 
values  are  ignored  and  the  component  max  and  min  are  used  with 
a  computed  increment  such  that  ten  contours  are  plotted. 

If  the  user  selects  the  relief  option, 
the  user  will  be  asked  to  input  the  number  of  lines  per  element 
(NLIN) .  NLIN  must  be  greater  than  1,  and  determines  how  fine 
the  relief  map  is.  The  number  of  lines  that  should  be  chosen 
depends  on  how  large  the  elements  appear  on  the  plots,  for 
larger  elements  more  lines  are  suitable. 


NLIN  =  5  NUN  =  3 


Figure  4.3.4.  The  CORE  command  allows  the  user  to  select  the 

number  of  lines  per  element  that  may  be  utilized 
in  highlighting  the  structure.  The  use  of  these 
lines  is  illustrated  in  Figure  4.3.5. 

The  last  parameter  requested  of  the  user 
is  the  relief  scale.  The  relief  scale  is  the  maximum  number  of 
units  above  or  below  the  element  surface  the  relief  map  can  be. 
If  the  contour  only  or  the  contour  and  relief  option  is  used  the 
contour  will  be  scaled  off  the  surface.  In  the  contour  and 
relief  option,  this  results  in  having  the  contours  drawn  on  the 
relief  map.  In  the  contour  option  this  results  in  having  the 
contours  drawn  scaled  off  the  surface. 

The  default  value  is  contour  only  with 
the  minimum  and  maximum  being  the  component  minimum  and 
maximum  and  the  increment  computed  so  ten  contours  are  plotted. 
The  default  for  the  relief  scale  is  zero. 

4. 3. 4. 5  CUBE 

The  CUBE  command  allows  the  user  to 
change  and  minimum  and  maximum  x,  y  and  z  coordinates.  This 
command  can  be  used  either  to  zoom  in  on  a  particular  section 
or  to  pan  out.  Using  CUBE  to  zoom  has  an  advantage  over  the 
ZOOM  command  of  having  the  plot  automatically  scaled.  Also 
included  in  the  CUBE  command  is  a  change  of  site;  for  more 
information  see  the  SITE  command. 


4. 3. 4. 6  DEFO  (Figures  4.3.5a-c  and  4.3.6a-b) 

The  DEFO  command  allows  the  user  to 
decide  to  plot  the  undeformed  structure,  the  deformed  structure 
or  both.  If  the  user  chooses  to  plot  a  deformed  structure, 
a  scale  factor  must  also  be  input.  The  displacements  will  be 
multiplied  by  the  scale  factor  before  the  displacement  is 
added  to  the  undeformed  structure.  The  outline  of  either  the 
deformed  or  undeformed  structure,  when  plotted  alone,  will  be 
dotted,  but  when  they  are  plotted  together,  the  undeformed  is 
in  solid  lines  and  the  deformed  is  dotted.  The  default  for 
this  option  is  undeformed. 

4. 3. 4. 7  ENTE  (Figure  4.3.6b) 

The  ENTE  command  enables  the  user  the 
choice  of  plotting  the  entire  element  outline  or  just  the 
designated  surface.  If  the  user  selects  to  plot  the  entire 
element  for  element  types  1,  2,  5,  6,  7  and  8,  a  rectangular 
prism  will  be  drawn.  See  SURF  command  for  further  information 
on  surfaces.  The  default  is  to  just  plot  the  surface  6. 


4. 3. 4. 8  EXIT 


program. 


The  EXIT  command  is  used  to  exit  the 

4.3.4. 9  EXPL  (Figure  4.3.7) 


The  EXPL  command  gives  the  user  a  method 
of  separating  the  elements  from  each  other  so  the  plot  can  be 
more  easily  understood.  If  the  user  chooses  to  plot  an  exploded 
view,  a  scale  factor  must  be  entered.  This  factor  is  defined 
as  being  1~B/A  where  A  is  the  original  distance  from  a  node  to 
the  centroid  of  the  element  and  B  is  the  desired  distance  from 
the  node  to  the  centroid.  It  is  suggested  that  the  factor  be 
between  .25  and  .1.  In  this  range  the  elements  are  usually 
distinguishable  and  yet  large  enough  to  be  seen. 
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Figure  4.3.7.  This  figure  defines  the  parameters  used  in 

shrinking  the  elements  for  exploded  views.  Refer 
to  the  EXPL  command. 

4.3.4.10  EYEP 

The  EYEP  command  allows  the  user  to 
change  his  eye  position.  The  eye  position  is  given  in  x,  y,  z 
global  coordinates,  that  is,  the  coordinates  used  to  define 
nodal  positions.  The  default  eye  position  is  100,100,100. 

The  user  should  be  careful  in  choosing  the  eye  position.  If 
the  eye  position  is  too  close  to  the  structure,  the  plot  will  be 
clipped  by  the  clip  plane. 

4.3.4.11  HELP  (Figure  4.3.8) 

The  HELP  command  gives  a  list  of  the 
commands,  and  a  short  description  of  their  function.  Figure 
4.3.8  lists  the  commands  and  their  default  values. 

4.3.4.12  LABE  (Figures  4.3.5,  4.3.6  and  4.3.9) 

The  LABE  command  enables  the  user  to 
label  the  plot  which  will  help  in  plot  identification.  The 
label,  as  seen  in  the  figure,  identifies  the  structure,  the 
type  of  plot,  the  stress,  strain  or  displacement  plotted,  the 
view  and  if  desired,  contour  values  and/or  a  subtitle.  The 
title  of  the  structure,  which  appears  at  the  top  of  the  plot, 
consists  of  the  first  three  lines  of  the  mpost  file.  These 
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CONTROL  COMMANDS 
*******  ******** 

COMMAND 

NEST 

NEUD 

CUBE 

TIME 

EXIT 

PLOT 

HELP 

RESE 

SUMM 

UAIT 

OPTIONS 

******* 


Figure  4.3.8 


FUNCTION 

PREPARE  NEU  STRUCTURE  DATA 

PREPARE  NEU  DATA  SET  FOR  PLOTTING 

CHANGE  MAX/MIN  UALUES 

PRINT  CURRENT  CPU  TIME  USED 

EXIT  PROGRAM 

PLOT  UITH  PRESENT  DATA 

PRINT  COMMAND  OPTIONS  AND  FUNCTION 

RESET  TO  DEFAUL”  OPTIONS 

PRINT  THE  CURRENT  OPTION  CONDITION 

PUT  PROGRAM  INTO  PAUSE 


LELE 

LABEL  ELEMENTS 

NO 

LAXS 

LABEL  AXIS 

NO 

LABE 

PLOT  TITLE  BLOCK 

NO 

ALEL 

PLOT  ONLY  SOME  ELEMENTS 

NO 

ENTE 

PLOT  ENTIRE  ELEMENT 

NO 

EXPL 

PLOT  EXPLODED  VIEU 

NO 

ZOOM 

EXPAND  AREA  OF  INTEREST 

NO 

POSA 

PLOT  A  POSITIVE  ARROW 

NO 

SUBT 

SUBTITLE  IN  TITLE  BLOCK 

NO 

DEFO 

SELECT  GEOMETRY 

UNDEFORMED 

CORE 

SELECT  CONTOUR  OR  RELIEF 

CONTOUR 

CLIP 

DISTANCE  TO  CLIP  PLANE 

.01 

EYEP 

SELECT  EYE  POSTION 

100,100, 100 

SITE 

CHANGE  SITE  POSTION 

CENTROID 

UERT 

SELECT  VERTICAL  AXIS 

Z 

REFL 

SELECT  A  REFLECT  OPTION 

NONE 

PROJ 

SELECT  PROJECTION  TYPE 

PERSPECT. 

SURF 

SELECT  PLOTTING  SURFACE 

6 

CONL 

LABEL  CONTOURS 

NO 

STEP 

CHANGE  STEP  SIZE 

.05 

NODE 

LABEL  NODES 

NO 

This  list  is  obtained  by  executing  the  HELP 
command. 
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Figure  4.3.9. 


This  illustration  demonstrates  the  LABE  command 
allowing  the  plot  to  be  well -documented  with 
necessary  plot  information.  The  CONL  command 
was  also  invoked  to  label  the  contours. 
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lines  originate  from  the  beginning  of  the  input  file  and  may  be 
up  to  eighty  characters.  The  type  of  plot  is  designated  by 
telling  if  the  plot  is  contour,  relief,  or  both.  The  components 
of  stress  or  strain  are  given  using  the  codes  described  in 
Section  4.3.3.  When  displacements  are  plotted,  DIS  is  used 
with  the  direction.  The  eye  position,  view,  minimum  and 
maximum  of  component  plotted  and  certain  plot  parameters  are 
also  given.  The  contour  label  values  and  subtitle  are  discussed 
in  the  CONL  and  SUBT  commands  respectively. 

4.3.4.13  LAXS  (Figure  4.3.5a) 

The  LAXS  command  allows  the  user  to  plot 
the  x,  y  and  z  axis.  This  can  be  useful  in  determining  the 
orientation  of  the  structure.  The  axis  will  be  drawn  from  the 
(0,0,0)  point  to  the  maximum  point  on  each  axis.  The  coordinates 

used  are  global  coordinates,  that  is,  the  coordinates  used  to 

define  elements. 

4.3.4.14  LELE  (Figure  4 . 3 . 9 ) 

The  LELE  command  allows  the  user  to 
number  the  elements.  The  elements  are  numbered  starting  with 
the  first  element  type  used  and  numbered  consecutively  through 
element  types.  When  a  reflect  option  is  used  the  reflected 
elements  are  numbered  from  where  the  unreflected  elements  ended. 
NOTE:  Element  numbers  remain  constant  from  preprocessor 

generation  through  analysis  and  into  postprocessor  plotting. 

4.3.4.15  NEST 

The  NEST  command  allows  the  user  to 
plot  more  than  one  structure  on  the  same  file.  The  NEST  command 
sets  all  the  parameters  to  their  default  values  and  returns 
to  the  beginning  of  the  program.  See  Section  4.3.2. 

4.3.4.16  NEWD 

The  NEWD  command  allows  the  user  to 
change  the  data  set  being  plotted.  The  NEWD  command  does  not 
change  the  plotting  parameters.  See  Section  4.3.2. 


4.3.4.17  NODE 


The  NODE  command  allows  the  user  to 
label  nodes.  Depending  on  whether  or  not  the  entire  element 
is  to  be  plotted,  either  only  the  surface  nodes  or  all  the  nodes 
may  be  plotted.  For  more  information  about  the  surface,  consult 
the  SURF  command. 

4.3.4.18  PLOT 

The  PLOT  command  generates  a  plot  with 
the  current  data  set  and  parameters. 

4.3.4.19  POSA  (Figure  4.3.10) 

The  POSA  command  allows  the  user  to 
determine  the  positive  direction  on  a  relief  plot.  When  an 


Figure  4.3.10.  A  positive  direction  arrow  will  be  drawn  by  the 
POSA  command  to  inform  the  user  in  which 
direction  positive  stresses,  strains  or 
displacements  occur. 
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arrow  is  desired,  it  will  appear  in  the  center  of  the  first 
element  plotted,  and  will  point  in  the  positive  direction.  For 
example,  this  option  should  be  used  when  there  is  a  question  as 
to  whether  a  stress  represents  compression  or  tension. 

4.3.4.20  PROJ  (Figure  4.3.11) 

The  PROJ  command  allows  the  user  to  chose 
between  an  orthogonal  or  a  perspective  view.  An  orthogonal  view 
occurs  when  the  lines  of  projection  of  an  image  are  drawn 
perpendicular  to  the  projection  plane  whereas  a  perspective 
view  occurs  from  drawing  projection  lines  from  an  image  straight 
to  the  eye  position  of  the  viewer  and  being  intercepted  by  a 
projection  plane.  The  perspective  view  tends  to  distort  the 


Projection  Lines 

* 


Projection  Lines 


Structure 


Structure 


Projection  Plane 


Projection  Plane 


Orthogonal  Plot 


Perspective  Plot 


Figure  4.3.11.  This  figure  illustrates  the  differences  between 
orthogonal  and  perspective  projections. 


structure,  so  the  elements  closer  to  the  eye  appear  larger  than 
the  elements  further  away  from  the  eye.  The  default  is 
perspective . 

4.3.4.21  REFL  (Figure  4.2.14) 

The  REFL  command  allows  the  user  to  plot 
the  structure  plus  a  reflected  image.  The  possible  reflect 
options  are  x,  y,  z  or  none.  If  x,  y  or  z  are  chosen  the  image 
is  reflected  along  the  indicated  axis. 

4.3.4.22  RESE 

The  RESE  command  resets  all  the  plotting 
parameters  to  their  default  values.  For  a  list  of  the  default 
values  see  the  HELP  command. 

4.3.4.23  SITE 

The  SITE  command  allows  the  user  to 
change  the  site  position  of  the  structure.  The  site  position  is 
the  position  in  the  structure  that  appears  in  the  center  of  the 
plot.  By  changing  the  site  position,  the  user  can  concentrate 
his  attention  on  a  specific  section  of  the  structure.  When 
using  the  CUBE  command,  the  site  position  is  adjusted  to  the 
centroid  of  the  new  maximum  and  minimum  unless  otherwise 
specified . 

4.3.4.24  STEP  (Figure  4.3.12) 

The  STEP  command  allows  the  user  to 
produce  more  accurate  contours.  The  step  size  is  the 
distance  from  one  point  on  a  contour  to  the  next.  The  distance 
is  measured  in  local  element  coordinates.  This  step  will  be 
in  either  the  positive  or  negative  x  or  y  direction.  The 
default  value  is  .05.  The  suggested  values  for  the  step  are 
between  .1  and  .01.  The  user  should  be  warned  that  by 
decreasing  the  step  size  the  CPU  time  increases. 


p 


yi 


Figure  4.3.12.  The  STEP  command  allows  the  user  to  set  the 

increment  value  for  x  and  y  steps  to  be  used 
when  plotting  contour  lines.  The  same  value 
is  used  for  both  x  and  y  direction  steps. 

4.3.4.25  SUMM  (Figure  4.3.13) 

The  SUMM  command  gives  the  user  a  list 
of  the  plotting  parameters.  It  is  used  to  see  if  the  parameters 
are  properly  set.  Refer  to  Figure  4.3.13  for  an  example. 

4.3.4.26  SURF  (Figure  4.3.6) 

The  SURF  command  allows  the  user  to 
select  which  surface  the  contour  or  relief  map  will  be  plotted 
on.  Only  one  surface  per  element  type  is  permitted,  but  each 
element  type  may  have  a  different  surface.  The  surfaces  are 
defined  by  the  order  in  which  the  connectivities  are  arranged. 
The  surfaces  are  numbered  one  through  six  as  shown  in  Figure 
4.2.9.  The  default  is  surface  six. 

NOTE:  The  figure  shows  element  type  one,  however,  the  surface 
of  the  other  element  types  are  defined  as  if  they  are 
degenerate  type  one  elements. 

4.3.4.27  TIME 

The  TIME  command  allows  the  user  to 
print  the  CPU  time  used  since  the  user  logged  in.  On  some 
systems  where  it  is  not  possible,  other  functions  such  as  time 
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Figure  4 


LABEL  ELEMENTS . NO 

PLOT  AND  LABEL  AXES . NO 

LABEL . NO 

PLOT  ALL  ELEMENTS . VES 

PLOT  ENTIRE  ELEMENT . NO 

EXPLODED  VIEW . NO 

POSITIVE  ARROW  .  NO 

PLOT  STRUCTURE . UNDEFORMED 

PLOT  OPTION . RELIEF 

MINIMUM  VALUE-  53.429  MAXIMUM  VALUE-  19361.131 

CLIP  PLANE  FACTOR .  9. 910 

EVE  POSITION! 

X-EVE  POSITION  .  25.000 

Y-EVE  POSITION  .  50.000 

Z-EVE  POSITION  .  10.000 

SITE  POSITION .  5.000  0.063  0.250 

AXIS  VERTICAL . Z-AXIS 

AXIS  OF  REFLECTION . NONE 

PROJECTION  TYPE . PERSPECTIVE 

MAXIMUM  X-  10.000  MINIMUM  X-  0.000 

MAXIMUM  Y-  0.125  MINIMUM  V-  0.000 

MAXIMUM  Z-  0.500  MINIMUM  Z-  0.000 

ZOOM  IN  . NO 

STEP  SIZE  IS  .  0.050000 

LABEL  NODES . NO 

? 


3.13.  This  figure  is  an  example  of  the  SUMM  command. 
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of  day  may  be  substituted.  This  function  is  especially  useful 
when  the  user  can  be  automatically  logged  out  when  a  certain 
CPU  time  limit  is  exceeded. 

4.3.4.28  VERT 

The  VERT  command  allows  the  user  to 
decide  which  axis  will  appear  vertical  on  the  plot.  The  default 
is  the  z  axis. 

4.3.4.29  WAIT 

The  WAIT  command  puts  the  program  into 
pause.  This  is  useful  when  the  user  must  leave  the  terminal  for 
a  period  of  time  and  does  not  wish  to  be  automatically  logged 
out  due  to  lack  of  action.  To  resume  the  program  execution,  the 
user  simply  types  GO  <CR>  (carriage  return) . 

4.3.4.30  ZOOM  (Figure  4.2.18) 

ZOOM  is  the  second  command  the  user  can 
employ  in  getting  a  blow-up  of  a  portion  of  a  structure.  The 
structure  must  first  be  plotted  with  the  view  desired.  Then  the 
user  specifies  a  virtual  window  which  encloses  that  part  of  the 
plot  on  the  screen  the  user  desires  to  view  more  closely. 

This  is  done  using  the  cursor  (two 
fine  crossed  lines  which  appear  on  the  screen  once  the  ZOOM 
command  is  implemented) .  The  user  adjusts  the  cursor  to  the 
location  of  the  lower  left-hand  corner  of  the  window  he  wishes 
to  create  using  the  cursor  controls.  He  then  enters  any 
alphanumeric  character  from  the  terminal  and  depresses  the 
"RETURN"  key  if  the  terminal  is  configured  to  require  a  carriage 
return.  He  repeats  this  procedure  to  define  the  upper  right- 
hand  corner  of  the  window  he  wishes  to  create.  The  user  may  now 
plot  the  'windowed'  area  at  full  size  on  the  screen  or  plotter. 

To  ZOOM  on  an  HP7221  bed  plotter  the 
user  must  wait  until  the  'ENTER'  light  on  the  plotter  control 
panel  begins  to  flash;  he  must  then  move  the  pen  to  the  lower 
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left  corner  of  the  zoom  area  with  the  pen  controllers  (buttons 
with  arrows  on  them) .  To  enter  the  coordinates  of  the  lower 
left  corner  the  user  must  press  the  'ENTER'  button.  This 
enters  one  x,  y  coordinate  pair.  The  user  must  then  move  the 
pen  to  the  upper  right  corner  of  the  area  to  be  expanded  with 
the  pen  controllers  and  again  press  the  'ENTER'  button.  The 
program  will  take  these  two  (x,y)  coordinate  pairs  and  create 
a  'window'  onto  which  the  zoomed  area  will  be  plotted.  The  user 
may  now  plot  the  structure. 

To  get  a  different  view  of  the  same  area, 
the  user  must  repeat  the  entire  process.  It  should  be  noted 
that  the  final  image  is  out  of  scale  (the  problem  is  similar  to 
the  one  of  the  virtual  plane  being  mapped  into  the  screen,  as 
described  in  SCALE) .  Refer  to  the  CUBE  command  for  an 
alternative  method  of  enlarging  a  specific  area  of  a  model.  It 
is  important  to  make  sure  your  terminal  or  plotter  has 
digitizing  capability  prior  to  executing  this  command. 
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GLOSSARY  OF  REPORT  TERMS 


aft  The  trailing  edge  of  the  wing.  This  edge  is 

closer  to  the  tail  of  the  plane  than  the  leading 
edge  which  is  closer  to  the  nose  of  the  plane. 
Section  References  2.8. 

analysis  The  process  of  conducting  an  investigation  of 

structural  flight  loading  simulation  of  a  finite 
element  model.  The  MAGNA  computer  program  is 
utilized  to  perform  a  flight  loads  simulation 
analysis  of  the  wing  specimens.  Section 
References  1.1,  1.3,  3.1. 

ATTACH  NOS/BE  command  to  gain  the  interactive  computer 

user  access  to  a  data  file  or  computer  program 
stored  on  auxiliary  storage  on  the  CDC  6600 
computer.  Refer  to  the  NOS/BE  INTERCOM 
manual  (Reference  2). 

AXES  Command  in  PLOTBOB  plotting  program.  This  allows 

x-,  y-,  z  axes  to  be  drawn  with  the  structure. 

axonometric  A  method  of  representing  three-dimensional  images 
on  a  two-dimensional  surface  (or  screen) .  The 
orthogonal  projection  in  the  plotting  programs 
utilizes  the  axonometric  method  where  lines  are 
drawn  from  a  figure  in  space  perpendicular  to 
the  surface  onto  which  they  are  to  be  displayed. 
Refer  to  PROJ  (CONTOUR  and  PLOTBOB)  commands  in 
Chapter  4.  Section  references  4.2.3.12,  4.3.4.20. 

bar  element  Rib  and  spar  caps  and  posts  are  represented  with 
one-dimensional  bar  elements  also  called  beam 
elements,  truss  elements,  or  MAGNA  type  4  elements. 
These  elements  are  defined  with  two  nodes  and  a 
thickness.  Section  References  2.3. 

basic  model  The  basic  model  will  be  generated  by  WINGEN 
utilizing  the  input  data  prior  to  damage, 
modifications,  refinements  or  conversion  to  shell 
elements.  The  basic  model  uses  type  3  and  4 
elements  (plate  and  bar)  and  is  composed  of  ribs, 
spars,  skins,  rib  and  spar  caps  and  posts. 

Section  References  1.1,  2.3. 

BATCH  mode  The  CDC6600  computer  allows  two  methods  for  user 

interaction:  a)  BATCH  mode  provides  the  user  with 

the  full  computing  resources  of  the  system  but  the 
user  cannot  control  any  aspects  of  what  the 
computer  does  except  by  input  data  cards.  All 
cards  are  assembled  and  read  into  the  computer 
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via  a  card  reader  or  a  "card-image"  file  is 
created  interactively  and  the  BATCH  command  is 
utilized  to  send  the  information  to  the  computer 
to  perform  a  task,  b)  Interactive  user  of  the 
CDC  system  may  be  accomplished  through  INTERCOM, 
a  facility  that  allows  the  user  to  share  a 
selected  portion  of  the  CDC  operating  system  to 
accomplish  a  range  of  tasks  while  exerting 
immediate  control  over  the  processes.  Interactive 
use  of  the  CDC  requires  an  interactive  computer 
terminal  while  batch  use  requires  a  card  reader 
or  interactive  computer  terminal.  Section 
Re f e rence s  1.2,  2.2,  2.5,  2.8. 

baud  rate  Computers  may  send  and  receive  data  over  telephone 

wires  or  "dedicated"  wires  at  various  rates 
depending  on  the  equipment  available.  The  baud 
rate  is  a  designation  of  bits  of  information 
transmitted  every  second  between  the  computer  and 
a  terminal  or  other  device.  The  most  common 
speeds  today  are  rates  of  300  and  1200  baud. 

The  CDC  generally  works  in  CPS  units  which  are 
characters  transmitted  per  second.  See  CPS  for 
further  information.  Section  References  2.4, 
4.2.1,  4.3.1. 

beam  element  See  bar  element. 

Once  a  computer  program  has  been  written  and 
tested,  it  can  be  advantageous  to  store  it  in 
a  binary  program  form  as  opposed  to  a  source 
program  form.  The  source  program  is  a  legible 
format  of  the  higher  level  commands.  After  the 
program  has  been  compiled  and  linked,  which 
converts  the  programming  language  statements 
to  machine  language  and  adds  to  the  program 
all  the  additional  modules  required  from 
program  libraries  the  user  will  have  a  binary 
program  file.  This  file,  while  generally 
larger  than  a  program  source  file  is  easier 
to  execute  because  no  compilation  is  needed. 

Large  programs  that  require  segmented  loading 
are  generally  handled  this  way.  Section 
References  2.4. 

boundary  Finite  element  analysis  requires  that  structures 

conditions  which  are  modelled  adhere  to  certain  conditions 
so  the  analysis  will  be  accurate.  One  condition 
is  that  the  structure  must  be  fixed  at  some 
point  or  there  will  be  no  distribution  of  a 
load  over  the  structure.  A  second  condition 
is  that  elements  must  be  compatible  such  that 
the  boundary  of  any  one  element  does  not  contain 


binary 

program 
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more  them  one  other  element  without  special 
constraints  being  specified.  Third  a  model 
should  not  have  any  nodes  defined  which  are  not 
used  in  element  connectivities  unless  the  node 
is  constrained.  Section  References  2.11.1. 

card  file  Batch  mode  utilization  of  the  CDC  may  be 

accomplished  through  the  use  of  card  decks  which 
are  submitted  to  the  CDC  control  room  or  read 
into  remote  card  readers.  A  card  deck  or  file 
will  have  all  the  necessary  job  control 
information  to  tell  the  computer  what  to  do 
with  it.  The  user  has  no  means  of  interfering 
with  the  execution  of  a  batch  job.  Section 
References  1.2,  2.2. 

Users  of  the  CDC  may  create  a  card-image  file 
of  a  card  file  and  store  it  on  auxiliary  disk 
storage  in  the  computer.  Once  stored  on  the 
computer  a  simple  ATTACH  command  will  retrieve 
the  file  and  make  all  the  data  or  program 
source  available  to  the  requestor.  Consult 
the  NOS/BE  operating  manuals  (Reference  2) . 

Section  References  1.2,  2.2. 

CATALOG  This  NOS/BE  command  when  used  with  the  REQUEST 

command  allows  the  user  to  make  a  file  "permanent" 
on  the  CDC  for  future  access  through  ATTACH.  A 
CATALOG' ed  file  may  be  removed  with  a  PURGE 
command.  Consult  the  NOS/BE  operating  manuals 
(Reference  2). 

CENTER  A  preprocessor  (WINGEN)  control  card  to  reset 

the  center  of  the  loaded  plane.  Consult 
section  2.8.12. 

centroid  The  geometric  center  of  a  two-  or  three- 

dimensional  structure.  The  centroid  is  used 
in  determining  if  an  element  should  be  included 
in  the  element  list  by  investigating  whether 
the  centroid  of  an  element  lies  within  the  region 
to  be  deleted  by  damage  or  modification  directives. 
Section  References  4.2.3.16,  4. 3.4. 9. 

The  transmission  rate  over  phone  lines  or  "hard¬ 
wired"  terminals  is  usually  given  in  characters 
per  second  (CPS) .  For  most  instances  the  CPS 
rate  of  interest  will  be  30  or  120  which  relate 
with  300  baud  and  1200  baud  respectively.  These 
rates  are  utilized  as  dummy  arguments  for 
initializing  the  plotting  software.  Section 
References  2.10,  4.2.2,  4.3.3. 
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A  chord  bay  is  the  space  encompassed  by  two  spars. 

A  chord  bay  is  numbered  by  counting  from  leading 
edge  to  trailing  edge  along  the  root  chord  (or  rib) . 
The  chord  bay  extends  from  the  root  chord  to 
the  wing  tip.  Section  Reference  2.8. 

The  chord  dimension  refers  to  the  root  chord 
length  in  the  x-axis  direction.  Section 
Reference  2.8.3,  2.8.4. 

The  chordwise  direction  is  always  from  leading 
edge  to  trailing  edge.  Section  Reference  2.8. 

A  chordwise  depth  station  is  a  point  along  any 
rib  where  the  wing  depth  changes  with  or 
without  a  concurrent  nodal  station  defined  for 
rib  and  spar  intersection.  Section  References 
2.8.3,  2.8.4. 

WINGEN  will  apply  a  chordwise  moment  load  in  the 
positive  x  direction  to  the  upper  skin  nodes  and 
in  the  negative  x  direction  to  the  lower  skin 
nodes.  Section  References  2.8.12,  2.11. 

Any  change  in  the  wing  depth  will  be  defined  by 
a  nodal  station  and  all  intersections  of  ribs 
and  spars  will  be  defined  as  a  nodal  station. 

All  the  nodal  stations  along  a  rib  collectively 
define  a  chordwise  nodal  line.  Section  References 
2.8.3,  2.8.4. 

WINGEN  will  apply  a  chordwise  shear  load  in  a 
positive  x  direction.  Section  References 
2.8.12,  2.11. 

At  every  point  in  a  model  along  a  rib  where  the 
wing  depth  changes  or  a  spar  and  rib  intersect 
there  is  defined  a  chordwise  station  at  which 
a  node  is  defined  (hence  a  chordwise  nodal 
station).  Section  References  2.8.3,  2.8.4. 

A  command  for  the  plotting  programs  to  allow  the 
user  to  cut  away  part  of  the  structure  for 
better  viewing. 

A  two-dimensional  plane  defined  in  a  three- 
dimensional  space.  All  elements  and  nodes 
between  the  clip  plane  and  the  viewing  position 
(EYE)  will  not  be  displayed,  e.g.,  they  will  be 
"clipped"  away.  Section  References  4. 2. 3. 2, 

4. 3. 4. 2. 


5.6 


I 


CONL 


CONTOUR 


control 

directives 


coordinate 
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coordinate 

transformation 


CORE 


CUBE 


cursor 
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A  CONTOUR  program  command  which  allows  the  user 
to  label  the  contours  of  a  contour  plot.  A 
listing  of  the  contour  line  values  will  also  be 
given  if  the  LABE  command  is  optioned. 

One  of  two  plotting  programs  for  the  graphical 
display  of  finite  element  model.  CONTOUR  is 
specially  defined  for  the  display  of  models  to 
illustrate  stress,  strain  or  displacements  in 
the  form  of  contours  drawn  on  the  models. 

Section  References  1.4,  4.3.  t 

The  preprocessor  program  WINGEN  requires  several 
pieces  of  information  to  allow  it  to  properly 
handle  the  model  generation.  This  information 
is  requested  at  the  beginning  of  the  program  and 
is  denoted  as  program  control  directives. 

Section  References  2.5,  4.3.3. 

All  finite  element  models  are  generated  as  a 
combination  of  three  coordinate  points  in 
space.  The  points  or  nodes  are  defined 
relative  to  some  arbitrary  coordinate  origin. 

For  the  model  generation  in  the  preprocessor 
program  the  coordinate  origin  is  suggested  to 
be  at  the  intersection  of  the  root  chord  and 
the  leading  edge.  Section  Reference  2.8. 

To  assist  the  user  in  visualizing  a  correct 
model  it  is  sometimes  necessary  to  alter  the 
coordinate  origin.  This  is  accomplished  as 
coordinate  transformation  where  all  coordinates 
will  be  shifted  according  to  the  user's  wishes 
resulting  in  the  model  being  oriented  in  space 
differently  but  the  model  itself  will  remain 
the  same.  Section  References  2.8.3,  4.2.3.20, 

4.3.4.23. 

A  CONTOUR  plotting  program  command  which  allows 
the  user  to  select  for  contour  or  relief 
plotting  of  a  structure. 

A  plotting  program  command  which  allows  the 
user  to  create  a  cube  around  a  model  or  section 
of  a  model.  All  elements  that  lie  outside 
the  viewing  cube  will  be  eliminated. 

A  cursor  is  a  positioner  for  locating  points  on  a 
viewing  screen  or  data  tablet.  It  consists  of 
fine  cross-lines  or  cross-hairs  to  aid  the  user 
in  locating  points  very  exactly;  cursor 
controlers,  if  the  cursor  is  generated  on  a 
screen  or  other  similar  device;  and  a  mechanism 
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cycle 


DAMAGE 


damage 


data  file 


i 


DEFAULT 


DEFO 

DEFORM 


to  transmit  the  data  point  selected,  a  key 
board  character  or  special  button.  Section 
References  4.2.3.22,  4.3.4.30. 

The  CDC  computing  system  allows  up  to  five 
permanent  files  to  be  stored  with  the  same 
name  for  each  problem  number/ID  number.  Each 
file  is  given  a  cycle  number  indicating  when 
it  was  cataloged  relative  to  all  othei  files 
with  that  name.  The  cycle  numbers  increase 
from  1  to  999.  Section  Reference  2.6. 

A  preprocessor  control  card  for  inflicting 
damage  to  a  wing  model.  Section  Reference 
2.8.11. 

A  modal  may  have  several  types  of  damage 
inflicted  to  it  when  the  model  is  being 
created.  These  damages  are  essentially 
removal  of  elements  to  simulate  structural 
damage  that  may  be  sustained  by  aircraft  in 
combat  situations.  Section  Reference  2.8.11. 

Most  computer  programs  require  the  user  to 
input  certain  pieces  of  information  which  are 
then  processed  and  the  desired  output  is 
returned  to  the  user.  In  most  cases  this  input 
information  will  be  typed  into  the  computer 
either  through  cards  or  by  an  interactive 
terminal  and  stored  as  a  group  of  information 
within  the  computer.  Such  a  group  of 
information  is  known  as  an  input  data  file. 

Once  a  program  has  finished  processing  the 
data  input  it  will  generally  create  an  output 
data  file  which  is  a  group  of  information 
awaiting  further  processing  via  another 
program  or  review  by  the  engineer.  Section 
References  1.1,  1.2,  1.3,  1.4,  2.2. 

A  PLOTBOB  plotting  program  command  which 
resets  all  plotting  parameters  to  their 
default  values.  Same  as  RESE  command  for 
CONTOUR. 

Plotting  program  commands  which  allow  the  user 
to  plot  a  deformed  or  undeformed  structure. 

A  preprocessor  program  control  card  which 
specifies  the  depth  of  the  wing  as  being  the 
distance  between  the  upper  and  lower  skins  of 
the  wing.  Section  Reference  2.8.4. 
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displacement 
of  structure 


DRAW 


edge  location 
edge  surface 
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A  user  may  input  to  a  requesting  program  certain 
information  concerning  an  image  displayed  on 
a  graphics  screen.  Digitizing  is  the  process 
of  utilizing  a  cursor  controller  to  position 
a  cursor  (or  locator)  composed  of  two  fine 
cross-hairs)  at  a  particular  position  then 
depressing  an  appropriate  control  to  send  to 
the  program  the  coordinates  of  the  cursor's 
position  in  digital  format.  The  program  will 
do  all  the  necessary  conversions  to  make  the 
digital  data  into  normal  cartesian  coordinates 
(x  and  y  values) .  Section  References 
4.2.3.22,  4. 3.4.30. 

When  a  structure  undergoes  analysis  due  to  the 
application  of  a  load  or  loads  then  certain 
structural  changes  or  deformations  are  most 
likely  to  occur.  The  displacement  of  a 
structure  refers  to  the  distance  a  structure 
or  parts  of  a  structure  are  moved  from  the 
initial  position  when  the  load  is  applied. 
Displacement  or  deformed  plots  will  illustrate 
the  deformation  of  the  loaded  structure. 

Section  References  3.3,  3.4,  4.2. 3.5,  4. 3. 4. 6. 

A  PLOTBOB  plotting  program  command  which 
causes  the  program  to  draw  the  structure 
utilizing  the  options  selected  by  the  user. 

Finite  elements  have  surfaces  and  edges  if 
they  have  more  than  one  dimension.  An  edge 
location  or  edge  surface  is  one  aspect  of  a 
finite  element.  Section  References  4. 2. 3. 7, 
4.3.4.26. 

Generally,  an  effective  area  is  that  region  of 
a  structure  or  model  over  which  a  localized  event 
will  have  an  influence.  An  effective  nodal  area 
is  that  region  over  which  a  load,  applied  to  the 
node,  will  be  distributed  with  the  greatest  effect. 
Section  Reference  2.11.1. 

Structures  are  composed  of  materials  which  may 
behave  differently  under  different  circumstances. 
Elasticity  and  plasticity  denote  two  different 
mechanisms  by  which  materials  will  respond 
under  differing  analytical  conditions.  Elastic 
behavior  is  typical  of  structures  undergoing 
linear  analysis  (small  displacements)  while 
elastic-plastic  behavior  is  representative  of 
material  undergoing  nonlinear  analysis  (large 
displacements).  Section  Reference  2.9. 
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element 


A  structure  may  be  represented  by  a 
mathematical  model  known  as  a  finite  element 
structure  or  model.  This  fern  is  composed  of 
elements  which  are  linked  together  via  common 
nodes  to  accurately  represent  the  model.  An 
element  is  composed  of  two  or  more  nodes  that 
define  a  one,  two  or  three-dimensional  space 
which  becomes  one  building  block  for  a  model. 

A  user  defines  as  many  of  these  building  blocks 
as  necessary  to  accurately  piece  together  a 
representative  model  of  the  structure  to  be 
analyzed.  Section  References  1.1,  2.5. 

element  types  All  finite  element  models  are  designed  for  a 

particular  analysis  program.  The  preprocessor 
WINGEN  was  designed  to  utilize  three  element 
types:  beam  (or  truss),  plate  (or  membrane)  and 

shell  for  analysis  by  MAGNA.  Each  element 
type  is  defined  in  Reference  2.  Section 
Reference  2.5. 

ELEMENTS  A  PLOTBOB  plotting  program  command  which  allows 

the  user  to  plot  selected  elements  and/or 
element  types.  Same  as  LELE  command  for 
CONTOUR. 

element  The  same  as  node  connectivities. 

connectivities 

ENTE  A  CONTOUR  plotting  program  command  that  allows 

the  user  to  plot  an  entire  element  or  only 
one  surface. 

EXIT  A  contour  plotting  program  command  to  stop 

the  program  execution. 

EXPL  A  CONTOUR  plotting  program  command  to  provide 

the  user  with  an  exploded  view  of  the  finite 
elements  of  a  model.  Each  element  is  isolated 
from  all  the  surrounding  elements. 

A  plotting  program  command  that  allows  the  user 
to  alter  the  angle  at  which  the  structure 
being  plotted  is  displayed.  Changing  the  eye 
position  is  analgous  to  holding  a  model 
airplane  in  one's  hand  and  rotating  it  or 
moving  it  closer  or  farther  away. 

failure  The  analysis  of  a  structure  is  generally 

concerned  with  whether  or  not  the  structure 
can  support  the  test  load.  If  it  cannot  then 
structural  failure  results  when  some  part  of 
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the  structure  no  longer  carries  the  load. 

Parts  of  structures  may  fail  without  failure 
of  the  entire  structure. 

Abbreviation  for  finite  element  model. 

A  finite  element  is  a  mathematical  concept  of 
linking  together  points  in  space  to  represent 
a  physical  object.  A  number  of  finite  elements 
connected  together  define  a  whole  structure. 
Section  References  1.1,  2.5. 

The  analysis  by  computer  simulation  of  a  finite 
element  model  to  determine  response  of  system 
to  some  type  of  loading.  Section  Reference 
1.1. 

A  structure  represented  by  a  group  of 
three-dimensional  coordinates  in  space 
which  are  connected  together  to  form 
finite  elements.  These  finite  elements 
represent  the  structure.  Section  References 
1.1,  2.5. 

A  program  that  conducts  a  finite  element 
analysis.  Section  References  1.1,  1.3,  3.1. 

A  preprocessor  option  for  applying  the  loads 
of  an  aircraft  wing  section  so  that  it  is 
distributed  as  a  function  of  the  nodal 
effective  area  consisting  of  wing  skins  and 
spars  of  the  clamped  load  end.  Section 
References  2.8.12,  2.11. 

The  tip  or  nose  of  the  airplane  is  always  the 
fore  (forward)  part  of  the  structure.  The 
fore  part  of  a  wing  is  that  edge  which  lies 
closest  to  the  nose  of  the  aircraft.  Section 
References  2.8. 

A  FORTRAN  programming  language  specification 
mechanism  for  inputting  data  into  a  program. 
Consult  a  FORTRAN  language  manual  for  further 
details.  Section  Reference  2.8. 

Specifies  that  an  analysis  will  presumably  go 
beyond  the  linear  response  range  for  the 
structure  being  analyzed.  Different  mathematical 
techniques  are  required  for  nonlinear  response 
of  a  structure  under  analysis.  Section 
Reference  2.9. 
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A  plot  or  graphic  picture  illustrating  the  basic 
components  of  a  structure  which  has  been 
represented  by  a  finite  element  model. 

Section  References  2.10,  4.1. 

A  plotting  program  command  which  gives  the  user 
a  list  of  valid  program  commands. 

A  library  of  FORTRAN  and  COMPASS  subroutines 
which  enable  the  pre-  and  post-processors 
to  plot  structures  on  a  Hewlett-Packard 
HP7221  bed  plotter.  Section  References 

2.4,  4.2.1,  4.3.1. 

From  any  reference  point,  inboard  will  refer  to 
a  direction  which  points  to  the  centerline  of 
the  fuselage  of  an  aircraft.  Section 
Reference  2.8. 

An  analysis  must  be  performed  as  a  function  of 
several  specific  points  of  reference  in  a  finite 
element  model.  MAGNA  will  select  a  minimum 
number  of  integration  points  within  an  element 
to  compute  the  various  stresses,  strains  and 
displacements  the  structure  will  undergo 
during  analysis.  Section  Reference  3.3. 

The  CDC  computing  system  has  facilities  for 
both  Batch  (remote  job  processing)  and 
INTERCOM  (interactive  job  processing). 
Interactive  mode  allows  the  user  to  get  an 
immediate  response  to  inquires  or  commands 
at  an  appropriate  terminal  from  the  computer 
system.  Section  References  1.2,  2.1,  2.2, 

2.5,  2.6.  Also  see  BATCH . 

A  CONTOUR  plotting  program  command  to  provide 
the  user  with  a  title  block  and  contour  label 
values  if  requested. 

A  PLOTBOB  plotting  program  command  to  allow 
the  user  the  option  of  labelling  the  nodes 
and/or  elements  of  a  model.  Same  as  LELE 
command  for  CONTOUR. 

A  CONTOUR  plotting  program  command  which  permits 
the  drawing  and  labelling  of  a  set  of  axes  for 
a  model.  Same  as  AXES  command  for  PLOTBOB. 

The  edge  of  the  wing  which  lies  closest  to  the 
nose  of  the  aircraft  is  the  fore  edge  or  leading 
edge  of  the  aircraft  wing.  Section  References 

2.5,  2.8. 
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A  CONTOUR  plotting  program  command  which 
allows  the  user  the  option  of  labelling  the 
nodes  and/or  elements.  Same  as  LABELS  command 
for  PLOTBOB. 

A  group  of  auxiliary  subroutines  which  have  a 
common  purpose  for  a  number  of  unique  programs 
will  generally  be  placed  in  a  special  file 
and  utilized  only  when  required  by  the 
programs.  This  file  is  known  as  a  library 
and  will  be  searched  by  the  system  linker 
for  subroutines  not  otherwise  provided  by 
the  programmer.  Section  References  2.4, 

4.2.1,  4.3.1. 

Under  very  slight  loading,  a  structure 
behaves  linearly  -  that  is,  the  strains  are 
linear  functions  of  the  displacements,  and  the 
stress-strain  law  is  of  linear  elastic  type. 

A  linear  analysis  is  performed  under  these 
assumptions,  and,  therefore,  represents  a 
"straight  line"  approximation  to  the  initial 
load-deflection  response.  Section  Reference 
2.9. 

Restrictions  on  how  nodes  behave  under 
loading  conditions  may  be  imposed  by  limiting 
node  displacements  to  be  a  linear  combination 
of  displacements  of  surrounding  nodes.  Section 
Reference  2.11.2. 

Any  group  of  letters  and/or  numbers  to  be 
used  exactly  as  they  appear.  Section 
Reference  2.8. 

A  load  is  any  force  applied  to  the  structure 
either  directly  or  indirectly.  Section 
Reference  2.8.12,  2.11. 

A  preprocessor  control  card  for  designating 
the  amount  of  loading  (breakdown  of  forces) 
to  be  applied  to  the  structure.  Section 
Reference  2.8.12. 

Several  items  must  be  determined  in  order  to 
distribute  the  loads  on  the  finite  element 
model  properly  such  as:  the  center  of  the 
load  plane,  the  magnitude  and  directions  of  the 
loads,  the  load  axis  of  the  wing,  etc.  These 
all  comprise  the  load  conditions.  Section 
Reference  2.11. 
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The  main  function  of  the  preprocessor  is  the 
creation  of  a  file  of  finite  elements  and 
loading  specifications  for  input  to  the  MAGNA 
finite  element  program.  Such  a  file  is  called 
a  load  deck.  Section  References  1.1,  1.2, 

1.3,  2.9,  3.2. 

The  MAGNA  analysis  program  allows  the  user 
to  specify  the  size  of  a  step  in  load  for 
each  iteration  of  a  solution  as  a  number  of 
total  iterations  to  be  used  in  performing  the 
solution.  A  load  step  is  that  fraction  of 
the  load  to  be  applied  during  each  solution 
increment.  Section  Reference  2.9. 

The  CDC  has  three  file  designations:  permanent, 
temporary  and  local.  A  local  file  can  be 
either  permanent  or  temporary  but  exists  only 
for  the  duration  of  an  INTERCOM  terminal 
session  or  BATCH  job  execution.  If  the 
local  file  was  permanent  it  is  returned  to 
its  storage  after  the  user  is  finished:  if 
it  was  temporary,  it  will  be  disposed  of. 

No  files  are  ever  automatically  made 
permanent  by  the  CDC.  Section  References 

2.4,  2.6,  2.7,  2.8,  2.9. 

A  CDC  command  requesting  the  computer  to 
acknowledge  your  desire  to  utilize  the 
INTERCOM  interactive  features  from  a  computer 
terminal.  Section  References  2.4,  4.2.1, 

4.3.1. 


A  computer  program  for  the  Materially  And 
Geometrically  Nonlinear  AnaTysis  of  finite 
element  models.  See  Reference  1.  Section 
References  1.3,  3.1. 

A  production  of  the  preprocessor,  the  load 
deck  contains  all  the  information  to  cause  a 
BATCH  execution  of  MAGNA.  Section  References 
1.2,  3.2. 

Same  concept  as  for  geometric  nonlinearities 
except  that  material  forces  are  the  subject 
instead  of  geometric  displacement.  Section 
Reference  2.9. 

A  preprocessor  code  which  indicates  what 
type  of  material  a  particular  part  of  a 
wing  structure  was  made  of.  Aluminum  and 
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steel  property  codes  (1  and  2)  are  generally 
used.  Section  References  2.8.5,  2.8.6,  2.8.7, 
2.8.9. 

Chordwise  moment  loading  applied  to  wing. 
Section  References  2.8.12,  2.11. 

Two-dimensional  elements  which  include  only 
in-plane  stress  and  strain  effects.  Membrane 
elements  utilized  in  this  program  are  of  two 
subtypes:  shear  panel  (subtype  3)  and  plane 

stress  (subtype  1) .  They  are  also  referred 
to  as  shear  panel  (for  ribs  and  spars)  and 
plane  stress  (for  skins)  elements.  Section 
Reference  2.5. 

Refers  to  the  size  of  the  finite  element  with 
respect  to  the  structure  being  represented. 

A  coarse  mesh  is  one  where  very  few  finite 
elements  will  define  a  section  or  entire 
structure.  A  fine  mesh  may  contain  as  many 
as  two  dozen  smaller  elements  where  only  one 
coarse  element  was  defined  previously. 

Section  References  1.1. 

A  finite  element  mesh  may  be  made  finer  by 
specifying  a  refinement.  Generally  a  two  or 
three  fold  increase  in  the  number  of  elements 
will  be  desired  in  a  refinement  but  more  are 
possible  as  well  as  the  feasibility  of  only 
increasing  the  density  (reducing  the  size) 
of  elements  in  selected  regions  of  the  model. 
Section  References  1.1,  2.8.10. 

A  group  of  finite  elements  which  define  a 
physical  structure.  Section  References  1.1, 
2.5. 

The  point  at  which  all  nodal  coordinates 
are  referenced  from.  Generally,  the  origin 
is  selected  to  be  at  the  intersection  of 
the  leading  edge  and  root  chord  of  the 
wing.  Section  Reference  2.8. 

A  graphic  illustration  of  a  finite  element 
model.  Section  References  2.10,  4.1. 

A  preprocessor  control  card  specifying  that 
certain  areas  of  a  generated  model  should  have 
certain  structural  components  omitted.  Section 
References  2.8.8. 
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An  output  data  file  generated  by  MAGNA  as  a 
result  of  the  structural  analysis.  This  file 
is  used  by  CONTOUR  and  PLOTBOB  for  postprocessing 
of  the  data  into  various  kinds  of  plots. 

Section  References  1.3,  1.4. 

Spanwise  moment  load  to  wing.  Section 
References  2.8.12,  2.11. 

Torque  moment  load  applied  to  wing.  Section 
References  2.8.12,  2.11. 

A  CONTOUR  plotting  program  command  permitting 
the  user  to  utilize  the  data  of  another 
structure  on  a  multiple-structure  data  file. 

A  PLOTBOB  plotting  program  command  that  permits 
the  user  to  utilize  the  data  of  another  structure 
on  a  multiple-structure  data  file. 

A  CONTOUR  plotting  program  command  that  allows 
the  user  to  plot  a  different  solution  increment 
for  a  structure  with  more  than  one  load 
increment  on  the  file. 

Three-dimensional  coordinates  located  in 
reference  to  an  origin  (the  model  coordinate 
origin)  which  are  associated  together  to  form 
building  blocks  of  regular  geometric  shapes 
known  as  finite  elements.  These  finite 
elements  are  defined  so  they  accurately 
represent  a  structure  to  be  analyzed.  Section 
Reference  1.1. 

A  CONTOUR  plotting  program  command  to  allow 
the  user  the  option  of  labelling  the  nodes 
of  a  structure  plot. 

Nodes  defined  as  three  coordinate  points  in 
space  relative  to  a  model  coordinate  origin 
are  associated  in  particular  sequences  and 
numbers  to  form  geometric  shapes  called 
elements.  These  elements  may  be  one,  two  or 
three-dimensional  in  nature.  By  incorporating 
nodes  into  several  elements,  element 
connectivities  are  formed  associating  one 
element  with  another.  Enough  associated 
finite  elements  are  defined  to  represent  a 
physical  structure.  The  association  of  nodes 
to  form  elements  is  called  nodal  connectivity. 
Section  References  2.8.5,  3.2. 
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A  nonlinear  structural  analysis  in  general  takes 
into  consideration  the  nonlinear  relations 
between  strain  and  displacements  (which  is 
necessary  when  the  displacements  are  finite) 
and  between  stress  and  strain  (due,  for  example, 
to  plastic  deformation) .  Since  the  governing 
equations  are  nonlinear,  a  solution  is  typically 
obtained  in  an  iterative  fashion  at  any  level 
of  loading. 

From  any  reference  point  the  outboard  direction 
is  from  that  point  away  from  the  center  line  of 
the  fuselage  of  the  aircraft.  Generally, 
outboard  refers  to  a  direction  from  the  wing 
root  to  the  wing  tip.  Section  Reference  2.8. 

The  CDC  has  three  file  designations:  permanent, 
local  and  temporary.  A  permanent  file  is  one 
the  user  has  made  a  request  to  be  permanent 
and  then  cataloged.  A  temporary  file  will 
be  created  only  for  the  duration  of  an  INTERCOM 
session  or  a  BATCH  job  execution  and  is  a 
local  file.  Any  permanent  files  used  during 
an  INTERCOM  or  BATCH  session  will  be  made  a 
local  file  until  the  session  is  over  and 
returned  to  its  storage  location.  Section 
Reference  2.5,  2.6,  4.2.1,  4.3.1. 

Same  as  membrane  element  subtype  1.  Section 
Reference  2.5. 

Two-dimensional  finite  elements  in  which  both 
inplane  and  transverse  bending  effects  are  taken 
into  account.  Out-of-plane  bending  deformations 
may  complicate  the  finite  element  by  necessitating 
the  use  of  rotations  or  similar  variables  as 
nodal  degrees  of  freedom. 

A  planform  description  indicates  the  general 
pattern  of  how  the  wing  is  situated  with 
respect  to  the  fuselage  of  the  aircraft: 
either  swept  fully,  partially  swept  or 
rectangular  (no  sweep).  Section  Reference  2.8.3. 

A  preprocessor  card  which  specifies  the  planform 
code.  Section  Reference  2.8.3. 

A  CONTOUR  plotting  program  command  that 
provides  the  user  with  a  graphic  presentation 
of  the  finite  element  model  currently  being 
postprocessed  utilizing  the  plotting  options 
selected  prior  to  the  command. 


5.17 


PLOTBOB 


PLOT 3D 


POSA 


postprocessor 


posts 


POSTS 


preprocessor 


preprocessor 
data  file 


profile  class 


PROFILE 


A  pre-  and  post-processing  graphics  program 
which  will  plot  a  MAGNA  load  deck  file  or  a 
MAGNA  MPOST  file  (WINGMPOST) .  Section 
References  1.4,  4.2. 

A  library  of  subroutines  designed  to  allow  the 
3-D  plotting  of  structures  on  a  2-D  surface. 

Section  References  2.4,  4.2.1,  4.3.1. 

A  CONTOUR  plotting  program  command  allowing  the 
user  to  have  an  arrow  plotted  with  the  structure 
indicating  the  positive  direction  of  what  force 
or  displacement  option  has  been  selected. 

A  program  which  accepts  data  from  a  finite  element 
analysis  program  and  further  interprets  or  reduces 
that  data  for  easier  interpretation,  generally,  in 
the  form  of  plots.  Section  References  1.1,  1.4,  4.1. 

Addition  of  vertical  one-dimensional  elements 
between  the  upper  and  lower  skins,  known  as  posts, 
are  to  give  the  model  structural  integrity  for 
analysis  purposes.  These  are  essentially  theoret¬ 
ical  constructs  to  allow  the  model  to  be  analyzed 
properly.  Section  References  2.5,  2.8,  2.8.9. 

A  preprocessor  control  card  specifying  the 
generating  options  for  the  posts.  Section 
Reference  2.8.9. 

A  program  which  utilizes  simplified  input  to  perform 
the  tedious  chores  of  generating  a  detailed  finite 
element  model  ready  for  analysis  by  a  finite  element 
program.  Section  References  1.2,  1.2,  2.1. 

The  preprocessor  WINGEN  has  an  optional  input 
data  file  for  wing  model  generation  which  contains 
abbreviated  wing  descriptions.  This  information 
will  be  used  to  generate  an  optional  output  data 
file  called  the  load  deck  which  is  a  detailed 
listing  of  nodes,  elements  and  all  other 
information  for  a  finite  element  model.  Section 
Reference  1.1. 

The  preprocessor  requires  the  user  to  input  a 
profile  class  for  the  wing  model  to  be  generated. 

This  profile  is  a  combination  of  the  planform  code 
and  the  chordwise  and  spanwise  depth  distributions. 
Section  Reference  2.8.2. 

A  preprocessor  control  card  which  defines  the 
profile  class.  Section  Reference  2.8.2. 
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A  plotting  program  command  which  allows  the 
user  to  select  a  perspective  or  an  orthogonal 
projection  of  an  image  onto  the  screen. 

An  image  is  defined  in  a  three-dimensional  space 
but  must  be  represented  on  a  two-dimensional 
surface.  The  projection  plane  is  defined  as  the 
surface  on  which  all  points  in  the  3-D  model  are 
imposed  or  mapped  to  yield  a  2-D  representation 
of  the  model.  Section  References  4.2.3.12,  4.3.4.20. 

A  CDC  BATCH  or  INTERCOM  command  used  to  eliminate 
unneeded  permanent  files.  Consult  Reference  2. 

A  model  where  the  original  mesh  generated  by 
connecting  the  basic  nodes  to  form  elements 
yields  a  normal  and  usually  a  coarse  mesh  model. 

A  refined  model  is  one  where  the  mesh  has  been 
made  finer  by  the  addition  of  more  nodes  and  the 
elements  rearranged  so  they  are  smaller.  Section 
References  1.1,  2.8.10. 

A  plotting  program  control  card  which  specifies 
which  areas  of  a  model  require  a  finer  mesh 
generation.  Section  Reference  2.8.10. 

A  plotting  program  command  which  enables  the 
user  to  draw  a  plot  of  the  model  plus  its 
reflection  about  an  axis.  This  is  helpful  if 
half  of  a  symmetrical  structure  is  modeled. 

A  model  may  have  its  mirror  image  plotted  along 
with  the  basic  model.  This  is  reflection,  which 
must  be  done  about  one  of  the  axis.  Section 
References  4.2.3.14,  4.3.4.21. 

A  CDC  BATCH  and  INTERACTIVE  command  which  allows 
a  file  to  be  initialized  so  that  it  may  be  made 
permanent  by  a  CATALOG  command.  Further 
information  may  be  found  in  Reference  2 . 

A  CONTOUR  plotting  program  command  which  resets 
all  options  for  plotting  to  their  default  values. 
Similar  to  DEFAULT  command  for  PLOTBOB. 

A  structural  component  of  a  wing  which  runs  from 
the  leading  edge  to  the  trailing  edge  of  a  wing. 
Several  ribs  will  be  present  in  any  given  wing. 
Section  References  1.2,  2.8,  2.8.6. 

A  preprocessor  control  card  designating  required 
information  for  rib  generation.  Section 
Reference  2.8.6. 


5.19 


rib  cap 

root  chord 


ROTATE 

rotation  of 
model 

SCALE 

semi- 

interactive 

mode 

shear  panel 
elements 

shell  element 

shell  element 
model 

SHRINK 


A  structural  component  of  ribs  is  a  stringer 
that  is  placed  in  conjunction  with 
the  rib  to  lend  structural  support  and  allow 
attachment  of  the  skins  at  the  rib  junctions. 

Section  References  1.1,  2.8,  2.8.6. 

The  root  of  a  wing  is  that  part  of  the  wing  model 
which  lies  closest  to  the  centerline  of  the 
aircraft  being  modeled.  The  root  chord  is 
the  first  nodal  station  defined  for  the  wing. 

The  root  chordwise  chord  runs  from  leading 
edge  to  trailing  edge  in  an  almost  parallel 
line  with  the  centerline  of  the  aircraft. 

Section  Reference  2.8. 

A  PLOTBOB  plotting  program  command  which 
allows  the  user  to  rotate  the  model  about  an 
axis . 

A  rotated  model  is  a  model  that  has  simply 
been  turned  so  that  a  different  part  of  the 
model  faces  the  viewer.  No  other  changes 
are  incorporated.  Section  References 
4.2.3.14,  4.3.4.28. 

A  PLOTBOB  plotting  program  command  that  allows 
the  user  to  alter  the  scale  of  a  plot, 
generally,  to  make  the  horizontal  and  vertical 
distances  proportional. 

The  WINGEN  preprocessor  has  three  operating 
modes:  batch;  interactive  and  semi-interactive. 

The  semi-interactive  mode  utilizes  the  batch 
mode  data  file  (Section  2.8)  but  executes 
under  user  interaction  for  plot  generation  and 
monitoring  of  model  generation  characteristics. 
Section  References  1.2,  2.2,  2.5,  2.7. 

Same  as  membrane  elements  subtype  3.  Section 
Reference  2.5. 

A  finite  elememt  type  with  three-dimensional 
characteristics  requiring  eight  nodes  to  define. 
Section  References  2.5,  2.8.10. 

A  finite  element  model  utilizing  shell  elements.  In 
this  situation,  shell  elements  compose  the 
top  and  bottom  skins.  Section  Reference  2.5. 

A  PLOTBOB  plotting  program  command  which  allows 
the  user  to  shrink  all  elements  about  their 
centroids  yielding  an  "exploded"  view  of  the 
model. 
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A  CONTOUR  plotting  program  command  which  allows 
the  user  to  select  a  different  point  within  the 
model  to  be  the  viewing  "center"  of  the  model. 

The  material  covering  the  frame  of  the  aircraft; 
most  notably,  the  material  covering  the  wing 
frame  composed  of  ribs,  spars  and  stringers. 
Section  Reference  2.8.5. 

The  chordwise  section  between  two  adjacent  ribs 
forms  a  span  bay.  Section  Reference  2.8. 

The  perpendicular  distance  from  the  root  chord 
to  the  intersection  with  a  line  passing  through 
the  wing  tip  chord  (or  rib) .  Section  References 
2.8,  2.8.3,  2.8.4. 

The  direction  from  the  root  chord  to  the  wing 
tip  (outboard).  Section  Reference  2.8. 

A  point  along  a  spar  between  the  root  chord 
and  wing  tip  where  there  is  a  change  in  the 
depth  between  the  upper  and  lower  wing  skins. 
Section  Reference  2.8.3,  2.8.4. 

A  point  along  a  spar  between  the  root  chord 
and  wing  tip,  inclusive,  where  a  node  has  been 
defined.  Section  Reference  2.8.3,  2.8.4. 

A  torque  force  which  acts  about  the  root  chord 
of  the  wing  with  the  wing  tip  -  MS.  Section 
Reference  2.11.1. 

A  vertical  displacing  force  upon  the  wing, 
positive  in  the  positive  z  direction  -  VS. 
Section  Reference  2.11.1. 

A  wing  structural  component  which  runs  from 
the  center  line  of  the  fuselage  to  the  wing 
tip.  Section  References  1.1,  2.8,  2.8.7. 

A  preprocessor  control  card  specifying  required 
information  for  spar  generation  in  the  model. 
Section  Reference  2.8.7. 

A  wing  structural  component  assocated  with  a 
spar.  The  cap  runs  adjacent  to  the  spar, 
lends  structural  stability  and  allows  a  surface 
for  wing  skin  attachment.  Section  References 
1.1,  2.8,  2.8.7. 
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A  CONTOUR  plotting  program  command  which  allows 
the  user  to  alter  the  computed  step  size  while 
calculating  contours. 

In  wing  design  it  sometimes  is  desirable  to  add 
lightweight  structural  supports  to  the  skins  to 
keep  them  from  "buckling".  Stiffeners  are 
strips  of  metal  that  run  along  the  skin 
generally  between  rib  and  spar  webs  to  lend 
additional  support  to  the  wing  skin.  Section 
Reference  1.1. 

MAGNA  analysis  of  a  fem  will  yield  results  which 
are  broken  down  into  various  classes  of 
stresses  and  strains  such  as  xx,  xy,  yz,  etc., 
in  terms  of  a  local  coordinate  system  for  the 
type  3  element  types.  The  CONTOUR  plotting 
program  will  request  the  user  to  select  one 
type  of  stress  or  strain  from  a  list  of  codes 
available.  Section  Reference  4.3.3. 

Any  material  constructed  in  such  a  way  as  to 
carry  a  load.  Section  Reference  1.1. 

A  method  by  which  a  structure  is  represented 
by  a  system  of  finite  elements  and  analyzed 
by  a  finite  element  program.  Section 
Reference  1.1. 

A  CONTOUR  plotting  program  command  which 
allows  the  user  to  specify  a  subtitle  to  be 
placed  on  a  plot  if  the  LABE  option  is  also 
selected. 

A  plotting  programs  command  which  provides  the 
user  with  a  list  of  the  plotting  options 
currently  selected. 

A  CONTOUR  plotting  program  command  which 
allows  the  user  to  specify  which  surface  of 
an  element  type  or  types  to  plot. 

A  preprocessor  program  local  file  containing 
input  data  for  the  generation  of  a  wing  finite 
element  model.  Section  References  2.7,  2.8. 

A  PLOTBOB  plotting  program  local  file 
containing  a  data  file  of  which  plots  will  be 
drawn.  Section  Reference  4.2.1. 

A  preprocessor  program  local  output  file 
containing  the  data  necessary  for  the 
regeneration  of  a  model  by  WINGEN.  This  is 


created  by  an  interactive  terminal  session  and 
is  made  a  permanent  file  by  WINGEN  prior  to 
program  execution  provided  cycle  space  is 
available.  Section  Reference  2.6. 

TAPE11  A  preprocessor  program  local  output  file 

containing  the  MAGNA  analysis  load  deck.  The 
user  must  catalog  this  file  if  he  wishes  to 
keep  it.  Section  Reference  2.9. 

TAPE99  A  CONTOUR  plotting  program  local  input  file 

containing  an  MPOST  data  file  for  postprocessing 
into  plots.  Section  Reference  4.3.1. 

TEKLIB  A  library  of  subroutines  utilized  by  the  pre- 

and  postprocessors  for  creating  plots  on  a 
Tektronix  graphics  device.  Section  References 
2.4,  4.2.1,  4.3.1. 

test  load  WINGEN  provides  the  user  with  two  means  of 

specifying  a  load  upon  a  wing  model  for 
analysis  by  MAGNA.  The  test  load  applies  the 
loads  over  the  wing  tip  based  on  an  associated 
end  wing  skin  area  for  each  node.  Section 
References  2.8.12,  2.11.1. 

TIME  A  plotting  programs  command  that  provides  the 

user  with  the  total  CPU  time  that  has  elapsed 
since  the  start  of  the  program  execution. 

title  The  preprocessor  requests  a  problem  title  to 

be  placed  on  the  data  files  for  identification 
purposes.  Section  Reference  2.8.1. 

torque  One  of  the  loading  forces  that  can  be  applied 

to  the  wing.  Section  References  2.8.12, 

2.11. 

TRANSLATE  A  PLOTBOB  plotting  program  command  which  permits 

the  user  to  move  the  model  from  the  origin  a 
designated  distance  in  any  direction. 

translation  The  process  of  altering  the  model  coordinate 

of  model  origin.  Section  References  4.2.3.20,  4.3.4.23. 

truss  element  A  bar  element. 

USRDATAFILE  WINGEN  will  provide  the  user  with  a  permanent 

file  named  USRDATAFILE  at  the  conclusion  of 
an  interactive  program  execution.  This 
permanent  file  will  be  present  as  a  local  file 
TAPE7  until  the  user  executes  a  LOGOUT.  Section 
Reference  2.6. 
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VERT  A  plotting  programs  command  which  allows  the 

VERTICAL  user  to  alter  which  axis  is  to  be  vertical  for 

plots. 

viewing  box  The  plotting  programs  display  of  3-D  image  on 

a  2-D  surface  utilizing  a  range  of  coordinate 
values  (from  min  x  to  max  x  and  min  y  to  max  y) 
to  create  an  imaginary  box  around  the  structure. 
Only  those  elements  that  lie  within  the  box  are 
plotted.  Section  References  4. 2. 3. 3,  4. 3. 4. 5. 

virtual  plane  Prior  to  plotting  a  structure  it  will  be  mapped 

from  a  three-dimensional  set  of  coordinates  to  a 
two-dimensional  plane  surface  which  is  then 
drawn  on  the  screen.  The  object  may  be  mapped 
using  one  of  several  projection  options  but  the 
result  is  always  a  two-dimensional  virtual 
plane  containing  the  final  image.  Section 
References  4.2.3.12,  4.3.4.20. 

WAIT  A  CONTOUR  plotting  program  command  that  places 

the  program  in  a  pause  state  until  the  user 
types  'GO'  or  the  computer  idle  times  (time  of 
last  terminal  input/output  signal)  is  exceeded. 
Currently,  idle  time  before  a  forced  logout  is 
five  minutes. 

window  Once  an  object  has  been  mapped  onto  a  virtual 

plane  (three-dimensional  structure  converted  to 
two-dimensional  structure)  it  is  ready  to  be 
displayed.  A  window  is  established  at  the 
graphics  device  which  allows  the  user  to  decide 
how  much  of  the  image  on  the  virtual  plane  will 
be  displayed  or  more  properly  what  range  of  x 
and  y  coordinates  in  a  two-dimensional  system 
will  be  actually  drawn.  Section  References 
4.2.3.22,  4.3.4.30. 

wing  depth  The  distance  (vertical  or  z  direction,  usually) 

between  the  upper  skin  and  the  lower  skin. 

Section  References  2.8.3,  2.8.4. 

wing  depth  A  preprocessor  requirement  for  wings  with  changes 

distribution  in  wing  depth  is  that  the  user  specify  in  a 

regular  pattern  how  the  wing  depth  changes  over 
the  entire  wing.  Section  Reference  2.8.4. 

wing  planform  Several  basic  wing  types  are  accommodated  for 

fern  generation  by  the  preprocessor  including 
the  classes  of  full  swept  (leading  edge  is 
straight  from  fuselage  to  wing  tip) ,  partially 
swept  (wing  is  swept  in  sections,  e.g.,  the 
leading  edge  is  not  straight) ,  and  rectangular 
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wing  profile 


wing  skin 
WINGEN 

yielding 


1 


(The  leading  edge  is  perpendicular  to  the 
fuselage) .  When  wing  sections  are  being 
analyzed  the  whole  section  as  it  rests  in  the 
test  frame  must  be  regarded  as  a  wing  where  the 
test  frame  is  the  fuselage.  Section  References 
2.8.3. 

The  preprocessor  requires  the  user  to  specify 
the  combinations  of  planform  and  depth 
distributions  over  span  and  chord  sections 
of  the  wing.  These  three  wing  specifications 
form  the  wing  profile.  Section  Reference 
2.8.2. 

See  skin. 

Preprocessor  program  which  inputs  simplified 
wing  geometries  and  generates  a  complete  finite 
element  model  file  ready  for  analysis  by  MAGNA. 
Section  References  1.2,  2.1. 

For  stresses  greater  than  a  certain  level 
(called  the  yield  stress) ,  permanent  deformations 
of  a  material  may  occur,  and  the  relationship 
between  strain  and  stress  becomes  nonlinear. 

The  onset  of  permanent,  or  "plastic",  deformation 
is  called  yielding.  In  two  or  three  dimensions, 
yielding  is  determined  in  most  metals  by  the  von 
Mises  criterion,  which  equates  the  distortional 
energy  at  a  point  to  the  corresponding  energy  in 
a  uniaxial  specimen  at  the  point  of  first  yielding. 


ZOOM 


A  plotting  program  command  which  allows  the  user 
to  plot  a  closer  view  of  a  section  of  a  model. 
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Section  F  -  login  procedure  for  executing  CONTOUR 
plotting  program  on  a  Hewlett-Packard 
bed  plotter  at  300  baud  (30  cps) 
utilizing  an  MPOST  file. 
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